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HIGH FREQUENCY PROPERTIES OF PLASMA 


PLASMA PHYSICS AND THE PROBLEMS OF CONTROLLED THERMONUCLEAR REACTION 


ABSTRACT 


These articles present the results derived from 
theoretical and experimental investigations of high 
frequency plasma properties: the methods of high fre- 
quency plasma heating, propagation of electromagnetic 
waves in a magnetoactive plasma, thermal radiation of 
a plasma, and development of instabilities when employing 
high frequency methods of plasma heating. A description 
is provided of experimental equipment developed for high 
frequency heating and containment of a plasma. 


This collection is designed for scientific researchers 
and engineers dealing with the problems of a plasma and 
its technical application, as well as for students and 
graduate students in the physics departments of universi- 
ties and physical-technical institutes. 


SECTION I 


HIGH FREQUENCY PLASMA HEATING 


INVESTIGATION OF THE ENERGY OF CHARGED PARTICLES EMANATING 
FROM A MAGNETIC TRAP DURING HIGH FREQUENCY HEATING 


N. I. Nazarov, A. I. Yermakov, V. T. Tolok 


The methods of resonance excitation of the eigen fluctuations of a [5* 
plasma cylinder by outer electromagnetic fields have been extensively em- 
ployed in plasma heating with high frequency fields. Spatially periodic 
electromagnetic fields may be employed to excite the eigen fluctuations in 
a plasma located in a constant magnetic field to frequencies which are close 


* Note: Numbers in the margin indicate pagination in the original foreign 
text. 


to the gyrofrequency of ions w 7% One (an ion cyclotron wave) or to frequencies 
Wyy < OK Ope (rapid magnetosound wave), where Wye is the electron gyrofre- 


quency and wyi is the ion gyrofrequency. In the first case, the energy of 


electromagnetic fluctuations is transmitted directly to the ions, and in the 
second case the energy is transmitted to the plasma electrons. 


The articles (Ref. 1, 2) have investigated the conditions for the reso- 
nance excitation of these fluctuations, their propagation, and damping. It 
has been found that even at high electron temperatures the damping is signi- 
ficant and is caused by a collisionless mechanism. It was also found that, /6_ 
when this excitation method is employed, high frequency power is transmitted 
to the plasma very effectively. 


This study presents the results derived from measuring the energy of 
ions and electrons passing along the magnetic field, when ion cyclotron and 
rapid magnetohydrodynamic waves are employed for heating the plasma. The 
experiment was performed on a "Sneg" apparatus, which has been described in 
great detail previously (Ref. 1). The eigen fluctuations were excited in the 
plasma by spatially periodic electromagnetic fields at a frequency of 10 Mc 
with the appropriate selection of the magnetic field strength Hj. In contrast 


to preceding experiments, the power of the high frequency (hf) generator was 
increased to 300 kw. 


In order to increase the power introduced into the plasma, the pulse of 
the hf generator was programmed so that the strong loading on the 
circuit at the moment of its resonance loading by the plasma was compensated 
by a corresponding voltage increase in the pulse from the hf generator (Figure 
1). Thus, the necessity of a special electrical strengthening of the hf cir- 
cuit was avoided, even when a power greater than 100 kw was introduced into 
the plasma. 


The energy of the charged particles was measured by a transit time electro- 
static analyzer (Ref. 3) and a multigrid probe (Ref. 4). The first method 
made it possible to study the energy spectrum and the mass composition of 
plasma ions; the second method made it possible to measure the energy of 
ions and electrons. The plasma electron temperature was determined by a 
spectral method. 


The input slit of the analyzer, which was located 25 cm behind a mag- 
netic mirror, cut out a narrow plasma flux, from which an ion bundle was 
separated after passing a separation device. The energy of the bundle 
was analyzed by the electric field of the flat condenser. Ion fluxes were 
recorded by an ion-electron converter, which changed the ion bundle into a 
bundle of electrons accelerated up to 20 kev, after deflection in the 
analyzing condenser. These electrons were detected by a plastic scintilla~ 
tor with a photoelectric multiplier. 


In order to study the mass composition of the plasma ions, the flight 


Figure 1 


time was measured by ions having a drift section 56 cm long. When the analy- 

zer modulator was supplied with voltage having a rectangular form, it was 
possible to obtain short pulses of the ion current (t = 1; 0.5; 0.2 microseconds). 
Due to a difference in the velocities of different ions, the ions were separ- 
ated by mass in the drift space. The flight time of ions in the drift sec- /7_ 
tion was used to determine the ion velocities and masses, respectively. By 
measuring the amplitude of the current signals at different periods of time 

and by changing the voltage on the analyzing condenser, it was possible to 

record the energy spectra of ions having different masses and to observe 

the change in ion energy during heating, by means of this analyzer. 


In order to exclude the scatter of ion current pulses, the result was 
averaged over ten measurements (the ion energy was determined with the 
analyzer with an accuracy of 8%). 


Figure 2 presents the oscillograms of a typical signal from the photo- 
multiplier, whose magnitude was proportional to the current of ions having 
an energy of 1500 ev — a; b — represents the suppression signal of micro- 
waves having a wavelength of 8 m. 


Figure 3 shows the energy spectra of plasma protons (the coordinates: 


the distribution functions Ņ -- ion energy E;) when the plasma is heated by 
means of ion cyclotron waves for two voltages on the hf circuit (curve 1 -- 
Ue = 28 kv; curve 2 -- U, = 32 kv). The plasma density to be measured during 


heating was no less than 2-1013 cm3’. As may be seen from the figure, the 
energy at the spectrum maximum amounts to 2 kev. With an increase in the 
voltage U, on the exciting coil of the hf circuit, the mean ion energy in 

the spectrum increases proportionally to the square of the voltage U2. How- 
ever, during the heating pulse it remains almost constant, which points to /8_ 
large losses which are apparently caused by overcharging. 


For purposes of comparison, the proton energy spectrum is presented which 
was recorded when the plasma was heated by a rapid magnetosound wave with a 
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Figure 2 


voltage of 32 kv on the hf circuit (Figure 4). The ion energy at the spec- 
trum maximum is 150 ev in all in this case. 


As may be seen from the ion mass-spectrograms (Figure 5) when the plasma 
is heated by an ion cyclotron wave, there are only hydrogen ions HH, Ht, 3 


in the plasma. The presence of these ions is apparently due to the fact that 
the plasma electron temperature is low (20-25 ev). It is particularly inter- 
esting to note that all three types of hydrogen ions have approximately the /9_ 
same energy, although the resonance acceleration occurs only for 1° 


For purposes of control, the ion energy was measured by another method 
-- multigrid probe with retarding potential (Ref. 4). This made it possible 
to record the energy spectrum of electrons. The multigrid probe was placed 
at a distance of 10 cm from the magnetic mirror. The plasma density was 
greatly reduced by means of a diaphragm having several 0.1 mm openings. It 
was possible to separate an electron bundle or an ion bundle, depending on 
the sign of the pulling voltage on the first grid; the bundles were analyzed 
by the retarding potential which was supplied to the second grid. In order 
to decrease the possibility of ionization, a differential pumping was employed 
to maintain a vacuum of 6.7:10-" n/m? within the probe walls. 
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Figure 5 


The results derived from measuring the ion energy spectrum by means of 
the multi-grid probe confirmed the results obtained with analyzer measure- 
ments. The proton energy spectrum maximum was in the region of 2 kev when 
the plasma was heated with ion cyclotron waves. The electron energy remained 
low, and amounted to 30 ev in all. This result also coincides closely with 
the results derived from measuring the electron energy by the spectral 
methods. 


As would be expected, the measurements of the energy spectrum of ions 
and electrons by the multigrid probe, when the plasma was heated with a 
rapid magnetosound wave, showed that the electron and ion energies were 
approximately the same (150 ev). The ion energies obtained by measurements 
with the probe and the analyzer also coincided fairly well. 
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MEASUREMENT OF THE PERPENDICULAR ENERGY COMPONENT AND THE /10 


PLASMA DECAY TIME DURING HIGH FREQUENCY HEATING 
N. I. Nazarov, A. I. Yermakov, V. T. Tolok 


The correct determination of the charged particle temperature is of 
paramount importance when studying the effectiveness of plasma heating. 
One convenient method of determining the plasma temperature consists of 
measuring its thermal diamagnetism by an external diamagnetic probe which 
includes a column of the plasma to be heated. The method is based on 
measuring the difference between the strengths of the magnetic field out- 
side and within the plasma AH, which is a function of the gasokinetic 


pressure. For a plasma with small 8 ¢ = | » this difference is deter- 
0 


mined by the expression 


z = 4mp 
AH = Ho ~ Byithin = Hy (1) 


where p is the gasokinetic pressure. With a plasma having quasineutrality, 
we have 


p= nk (Tat + T se), (2) 


where n is the plasma density; k is the Boltzmann constant. Tii and Tje 


are the perpendicular ion and electron temperatures, respectively. Thus, 
by measuring AH and knowing the plasma density, we may compute its tempera- 
ture. 


The temperature was measured by this method on a "Sneg" apparatus 


Figure 1 


(Ref. 1). Similar measurements were performed in (Ref. 2). The plasma 
was either produced by a powerful ion cyclotron wave (Ref. 1), or by a 
rapid magnetosound wave (Ref. 3). Hydrogen in the 0.13-0.4 n/m? pressure 
range was used as the process gas. The waves were excited at a frequency 
of fọ = 107 cps, and the pulse high frequency power, transmitted to the 
plasma, was 100 kw. [In order to avoid transitional processes related to 
the sharp change in plasma density during the initial period, a coupled, 
high frequency pulse, whose envelope is shown in Figure 1, was supplied to 
the exciting coil. The first pulse produced a plasma, and the second pulse 
was employed to heat it. The duration of the pulses, the amplitude, and 
the interval between them could be changed independently over very wide 
intervals. 


The period during which the strength of the pulse magnetic field Hg /11 
changed was 24 microseconds. The strength of the field was selected so 
that either the ion cyclotron wave, or the rapid magnetosound wave, was 
excited resonantly at the time the second high frequency pulse came into 
operation. The plasma density was measured by a microwave interferometer 
at the wavelengths 8.2 and 4 mm. 


The quantity AH, which was caused by the plasma diamagnetism, was de- 
termined by measurements of the electromotive force (emf) by a diamagnetic 
probe. The probe consisted of two coils, one of which included the plasma 
column. The other coil was employed to compensate for the emf caused by a 
change in the strength of the confining magnetic field Hg. In order to 
eliminate the emf produced by the propagated ion cyclotron wave, a five-unit 
low-frequency filter was employed, which intersected all frequencies above 
3 Mc. The probe was located in the region of the "magnetic beach", at a 
distance of 30 cm from the edge of the exciting coil. In order to decrease 
the effect of attenuation of the diamagnetic signal, due to reflection of 
the charged particles from the walls, the diameter of the discharge chamber 
was increased up to 8 cm in the region of the diamagnetic probe, while the 


Figure 2 


diameter of the plasma with n >10!2 cm-3 equalled 3.5 cm. 


Figure 2 shows an oscillogram of the diamagnetic signal -- a, obtained 
when the plasma was heated with an ion cyclotron wave; b -- represents the 
interferogram of an 8-mm signal. The gasokinetic plasma pressure increased 
very rapidly (in~10 microseconds), and then barely changed until the end 
of the high frequency pulse. Since the plasma density changed very little /12 
-- (1.2 - 1.5) x 10!3 cm? -- during the high frequency pulse, and since 
it decreased slowly after it had ended (t~ 270 microseconds), it can be 
assumed that the diamagnetic signal was proportional to the rate at which 


the perpendicular energy component of the plasma L (Tii + Tle? changed. 


Figure 3, a shows an oscillogram of a diamagnetic signal. In order to deter- 
mine the plasma temperature, AH was calculated after integration of the dia- 
magnetic signal. The integrated diamagnetic signal is shown in Figure 3, b. 
The total value of T, thus obtained for resonance excitation of an ion cyclo- 
tron wave, amounted to 1 kev. The electron temperature, determined by the 
spectral method, was in this case 20-30 ev. Thus, the ion temperature was 
measured indirectly. The plasma temperature was determined by this same 
method when it was heated by a rapid magnetosound wave. In this case 

T = 200-300 ev. The values of T, obtained according to the diamagnetic 
signal, closely coincide with the measured energy of charged particles em- 
inating from the system along the magnetic field (Ref. 4). The small di- /13 
vergence between Ti and T), is caused by the fact that measurements of the 


plasma temperature according to the diamagnetic signal give an average (over 
the plasma column cross section) temperature. 


Ho 
Figure 4 shows the dependence of the ion temperature T; on me (Hp -- 
ci 
the strength of the outer magnetic field, Hei -- the strength of the magnetic 


field at which the gyrofrequency of a proton equals 10 Mc). It can be seen 
that maximum heating occurs with resonance excitation of an ion cyclotron 
wave. 


Figure 4 


The ion temperature barely changes during heating (see Figure 3). 
There is a very rapid temperature decrease after the high frequency pulse 
is recorded. This change in the ion temperature points to the presence 
of great losses. Since the ion temperature is considerably greater than 
the electron temperature when this method of plasma heating is employed, 
it may be assumed that one of the mechanisms for rapid ion cooling is their 
energy loss when colliding with electrons. However, for a plasma with 
De~ 2-10!3 and Te~ 25 ev, the time required to cool hot ions must be 100 
microseconds. In actuality, the cooling time equalled 10 microseconds. 
Therefore, there is no basis for assuming that this is the main loss 
mechanism. 


Another mechanism for rapid ion cooling may be their overcharging, 
since under the experimental conditions the highly-ionized plasma column 
was surrounded with a weakly-ionized cold plasma which was in contact with 
the discharge chamber walls. For a plasma with an electron density on the 
order of 10/3 cm-3 and a neutral gas pressure at the chamber walls of 
approximately 1.33+10-2 n/m2, the probability of overcharging exceeds the 
ionization probability. Consequently, the losses to overcharging may be 
considerable, and they continue to increase with an increase in the ion 


Figure 5 


energy. Figure 5 shows the dependence of the ion cooling time on their 
energy. The solid line corresponds to the time required for overcharging 
the ions for a density of the surrounding gas of ng = 6:1012 em-3. This 
density was selected so that a comparison could be made between the value 
obtained theoretically and the experimentally measured value TTi for 


T; = 200 ev. The nature of the dependence of tT, on Tį shows that under 
experimental conditions the energy losses are primarily caused by ion 
overcharging. 


Thus, these experiments enable us to draw the conclusion that a /14 


plasma may be heated to a temperature exceeding 1 kev by means of resonance 
excitation of an ion cyclotron wave. The limiting value T; is determined 


by the apparatus parameters. In addition, the results obtained provide a 
basis for assuming that -- when a hot plasma is insulated from the chamber 
walls by "vacuum interstratification" -- the time the plasma may be contained 
may be increased considerably, under the condition that there are no other 


loss mechanisms. 
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HIGH FREQUENCY ENERGY ABSORPTION BY A PLASMA IN ION CYCLOTRON 15 
RESONANCE IN STRONG, HIGH FREQUENCY FIELDS 


V. V. Chechkin, M. P. Vasil'yev, L. I. Grigor'yeva, 
B. I. Smerdov 


This article represents a continuation of studies we performed pre- 
viously (Ref. 1, 2) on high frequency energy absorption by a plasma in 
ion cyclotron resonance. 


As is well known (Ref. 3-5), the heating of a plasma by a variable 
field at a frequency which is close to the ion cyclotron frequency is very 
effective, if there are mechanisms leading to energy thermalization of the 
orderly motion of plasma particles in the field of an ion cyclotron wave. 
High frequency energy absorption by the plasma may occur, in particular, 
due to "close" collisions of ions which are in resonance with other types 
of ions, electrons, and neutral atoms in a cold plasma. This absorption 
may also be due to "collisionless" cyclotron damping which is caused by the 
thermal motion of ions in a high temperature plasma. In both cases, it is 
assumed that the amplitude of the high frequency field is fairly small, so 
that ions receiving energy in the wave field can transmit it to other parti- 
cles. If this condition is not fulfilled, nonlinear processes must arise 
in the plasma, due to which the ion distribution function is essentially 
distorted. Distortion of the ion distribution function by the ion cyclotron 
wave with a finite amplitude leads to a decrease in the wave absorption co- 
efficient down to a small value which equals, in order of magnitude, the 
absorption coefficient for pair collisions at a given temperature. 


The expression obtained in the case of a high-temperature plasma 


(Ref. 5) for a critical field strength of an ion cyclotron wave -- which 
leads to a significant distortion of the ion distribution function when 
it is exceeded -- can be written as follows 
Ric H 8nN,T; 1/3 

Eor~ eae Se (wr) (I) (1) 
where w is the wave frequency; ky = on, A -~- axial wave length; H -- con- 
stant magnetic field; t -- relaxation time of ions due to ion-ion colli- /16 
sions; Ni -- ion density; T; -- ion temperature. 


11 


In the case of a cold plasma, the critical strength of the wave 
electric field (at the absorption maximum) is as follows 


= Ti, (2) 
eYeff M 


When this value is reached, the velocity of the ordered motion of an ion 
liquid, with respect to electrons, equals the thermal ion velocity. M -- 
ion mass; e -- ion charge; Yeff ~- effective frequency of ion collisions. 


If the velocity of the relative motion of ions and electrons in the 
field of an ion cyclotron wave is greater than the ion thermal velocity, 
"bunched" instability may arise in the plasma, which is related to the excita- 


tion of high frequency (as compared with w= $2 ) longitudinal fluctuations, 


whose increasing increment is considerably greater than the cyclotron ion frequency, 
and the wavelength is considerably less than the ion cyclotron wavelength. 
The excitation of these fluctuations by an ion bundle moving in a direction 
which is perpendicular to a constant magnetic field has been investigated in 
(Ref. 6-8). 


Similar small-scale electrostatic plasma fluctuations, which take place 
in a wide frequency and wave number range, must lead to an increased exchange 
of energy between plasma ions and electrons (as compared with the exchange 
caused only by Coulomb collisions), and must also lead to a significant in- 
crease in all the transfer coefficients across a constant magnetic field. 


Let us examine certain results of experimental studies in this light 
(Ref. 9, 10). In these studies, the field strength of an ion cyclotron 
wave exceeds byat least one order of magnitude the critical field strength 
(1), and the time of ion-ion relaxation is either a comparable with the 
plasma decay time (Ref. 9) or considerably exceeds it (Ref. 10). For this 
reason, the strong high frequency energy absorption by the plasma and the 
ion heating, observed in these studies, cannot -- in our opinion -- be 
caused by cyclotron damping. In addition, the study (Ref. 11), which was 
carried out on the same apparatus as was employed in (Ref. 10), observed a 
rapid plasma decay if the high frequency power introduced into the plasma 
exceeded a certain critical value. This rapid decay was apparently caused /17 
by "bunched" instability. 


The phenomenon of anomalous plasma diffusion across a magnetic field, 
produced for a critical value of a high frequency field strength with a 
frequency close to ion cyclotron frequency, was discovered and studied in 
detail in (Ref. 2). In particular, this study found that increased diffu- 
sion occurs if the velocity, acquired by ions in an azimuthal high frequency 
field during the period between collisions, is comparable to the thermal ion 
velocity or exceeds it -- i.e., if relationship (2) is fulfilled. 


The study (Ref. 12) performed an experimental determination of the 
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increase in the effective frequency of ion collisions in a low-density 
plasma under conditions of ion cyclotron resonance in strong, electric, 
variable fields. It was found that the dependence obtained experimentally 
for the frequency of ion collisions on the high frequency field strength 
cannot be explained on the basis of a theory postulating ion collisions 
with neutrals. The assumption was advanced that such a dependence is 
caused by"bunched" instability which arises in strong variable fields. 


We can clarify the nature of this absorption by comparing the data 
obtained experimentally, regarding high frequency energy absorption by 
a cold plasma close to ion cyclotron resonance, with the theory advanced 
in (Ref. 5). Let us establish a relationship between the effective fre- 
quency of ion collisions and anomalous plasma diffusion in ion cyclotron 
resonance, which was studied in (Ref. 2). Finally, by making certain 
numerical estimates, we can show that both high frequency energy absorp- 
tion by the plasma, and anomalous plasma diffusion in fields with super- 
eritical strength values, may be caused by "bunched' instability produced 
in the field of an ion cyclotron wave (this instability was investigated 
theoretically in [Ref. 8]). 


Description of the Apparatus. 


Measurement Methods 


The apparatus which was employed for the study was described in de- 
tail in (Ref.1). The plasma was produced by pulse discharge with oscil- 
lating electrons in hydrogen, in a 7-10-2— 7 n/m2 pressure range. The 
diameter of the glass discharge tube was 6 cm, and the distance between /18 
the cathodes was 80 cm. The strength of the longitudinal, quasi-constant 
magnetic field could be changed between 4-10'—6.4-10° a/m. 


High frequency energy was introduced into the plasma by means of 
an artificial LC-line which was connected to the self-excited oscillator with 
a frequency of 7.45-°10° cps. During excitation at this frequency, 2.5 wave- 
length oscillation was applied to the section of the line which was slipped on 
to the discharge tube. This corresponded to an axial high frequency field 
period of 23 cm. By changing the anode strength of the oscillator, and also 
the connection between the line and the oscillator, it was possible to change 
ae amplitude of the high frequency azimuthal line current within 0.5 - 35 
a/cm. 


The oscillator was switched on for approximately 100 microseconds 
after the pulse of the discharge current had terminated. As is shown in 
(Ref. 1), the maximum high frequency power which could be absorbed by a 
plasma in resonance was 18 kw for a plasma density of 1.71013 em=3, and 
an aZimuthal line current of 30 a/cm. 


This article presents the measurements of plasma density, electron 
temperature, and the high frequency power absorbed by the plasma. The 
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electron plasma density and its change with time in the (1.7 - 0.25) 10!3 
cm? range was measured by means of an interferometer at a wavelength of 
8.1 mm = 0.81 cm. The experimental point corresponding to a density of 
3-10!3 cm? was obtained by extrapolation of the plasma density dependence 
on time toward large densities (Figure 1). 


The plasma electron temperature was determined according to the de- 
pendence of the luminosity intensity of the line Hg on the high fre- 


quency LC-line current [see (Ref. 2)]. The luminosity intensity of the 
line He was measured by means of a UM-2 monochromator and a photoelectron 


multiplier. Due to anomalous diffusion, the plasma decay time, for line 
currents on the order of 5 a/cm and above, was comparable with the electron 
lifetime between two collisions leading to excitation of a neutral atom. /19 
For this reason, for line currents which were greater than, or approxi- 
mately equal to, 5 a/cm, the electron temperature which was computed 
according to the line intensity Hg could be too low, 


The high frequency power absorbed by the plasma was measured with an 
all purpose meter of transmitted power which was described in (Ref. 13). 
The recorders for the current and the strength were connected to the line 
at the point where it was attached to the solenoid. By means of two inter- 
changeable current recorders with different sensitivity, it was possible to 
measure the power absorbed by the plasma, corresponding to a line current of 
0.5 - 10 a/cm. 


All of the measurements described in this article were performed at 
a hydrogen pressure of 0.133 n/m~. 
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Nature of High Frequency Energy Absorption 


Under the conditions of our experiment, the high frequency power 
absorbed per unit of plasma cylinder length is (Ref. 5) 


=FR ki (ki R) (em jof (X) sin?kyz, (3) 


where jọ is the amplitude of the azimuthal high frequency current in the 
coil per unit of cylinder length; R -- coil radius; a -~ plasma radius 
(it is assumed to equal the inner radius of the discharge tube in all of 


the computations); Kı ~- McDonald function; kı = = -- axial wave number 

(it is determined as the axial period à of high frequency current in the 
k 

coil); w -- oscillator frequency; ny = ll £(X) -- the function whose 


specific form depends on the nature of the high frequency energy absorp- 
tion. In particular, for collision absorption we have 


Xy ` We 
F(X) = effi (4) 
] (Xot o— o) +Y eff 


where 


Strictly speaking, relationship (3) is only valid in the case of /20 
long-wave fluctuations (or a small plasma filament radius), when k]| a«1 


and ky) a<< 1 (k, -- radial wave number). Under our conditions, both of 


these quantities are on the order of unity. However, as the computation 
showed, the error produced when equation (3) is applied to our case is 
small. 


Under the conditions of the described experiment, high frequency 
energy absorption by a plasma close to w = wy is "collision absorption" 
in the sense that it is described by (3), where £(X) has the form of (4). 
Formulas (3) and (4) may be used to determine the fact that, in the case 
of a hydrogen plasma, the absorbed power is at a maximum in the case of 


LSE = 2,44- 10~12N a (5) 
H 


0 
where — = —; Hg is the magnetic field corresponding to cyclotron 
Hn 
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resonance of ions at the oscillator frequency; Hp -- the magnetic field 
forwhich the high frequency power absorbed by the plasma is at a maximun. 


The dependence described by (5) is shown in Figure 1 (straight line). 
This figure also presents the points computed on the basis of experimental 
data regarding the shift in the absorption maximum from Hg for different 


plasma densities measured by an interferometer for a line current of 

4.5 a/cm. Over the entire range of measured plasma densities, the devia- 
tion of the points obtained experimentally from the computed dependence 
does not exceed the limits of measurement errors. 


In order to determine what collisions cause the observed high fre- 
quency energy absorption, a study was made of the dependence of the 
effective frequency of ion collisions on the neutral gas pressure and on 
the plasma electron density. It can be seen from (3) and (4) that the 
effective frequency of ion collisions Yefg can be expressed either as 


Poe AH -- where AH is the halfwidth of the resonance absorption curve 


2 Me 
AXwy 2 2 

-- or as jo, where A is the factor in front of jgf(X) in (3); Sm -- 
m 


the power at the absorption maximum. 


It was shown in (Ref. 1) that the power at the absorption maximum 
and the width of the resonance absorption curve depend slightly on the /21 
neutral gas pressure in the 0.133 - 1.33 n/m“ range, due to which fact 
the observed power absorption cannot be caused by ion collisions with 
neutrals. 


Figure 2 presents the dependence obtained experimentally of the 
effective frequency of ion collisions y,¢¢ on the plasma electron density 


for a line current 4.5 a/em. The crosses designate the points computed 
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according to the halfwidth of the absorption curves; the dots designate 

the points computed according to the absolute power at the absorption maxi- 
mum for a given density. The observed high frequency energy absorption is 
caused by ion-electron collisions. However, under the conditions for 

which the dependence shown in Figure 2 was determined, the electron tem- 
perature comprised 1 ev in order of magnitude (see Figure 3). The frequency 
of ion Coulomb collisions with electrons must be one order of magnitude 

less than Yefrf over the entire range of measured plasma densities. There- 


fore, the effective frequency of ion collisions, which was measured in the 
experiment, cannot be caused by Coulomb scattering of ions by electrons 
[the electron temperature, which was determined in (Ref. 1) according to 
Yeff under the assumption of Coulomb interaction of ions with electrons, /22 


amounted to 0.15 ev). 


In order to clarify the nature of the interaction between ions and 
electrons, and consequently the nature of the observed high frequency 
energy absorption by a plasma in ion cyclotron resonance, the dependence 
of the effective frequency of ion collisions on the line current density 
jo ~~ i.e., the dependence on the strength of the high frequency field in 


the plasma -- was determined. The quantity Yegg Was computed according to 
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the power at the absorption maximum for a density of 5.1°10!2 cm3, 
Figure 3 shows the dependences Yoff (jo), Tojo) and t(jo) -- where T is 


the plasma decay time -- on density 7.6-10!2 up to 2.5°10!2 em-3 for a 
magnetic field strength of 4.2-10° a/m, where the absorption reaches a 
maximum for a density of 5.1-10!2 cm-3, It can be seen that Yefg rapidly 


decreases with an increase in jọ(jọ< 1 a/cm). The frequency of Coulomb 
collisions between ions and electrons must follow this pattern, since the 
high frequency power introduced into the plasma -- and consequently the 
electron temperature -- increases with an increase in jg. The electron 


temperature, determined according to Yeff in the case of j,< I a/cm under 


the assumption of Coulomb interaction between ions and electrons, coin- 
cidesin order of magnitude with the temperature computed according to 
the line intensity He (see Figure 3). 


There is a sharp bend in the curve Yoff (jo) at the point jg < 1 a/cm, 
and with a further increase in jg, Yeff increases slowly, remaining at a 


level on the order of 2'10 secl and considerably exceeding the frequency 
of Coulomb collisions between ions and electrons for given electron tempera- 
tures. As can be seen from Figure 3, at this point a sharp decrease in the 
plasma decay time occurs. It was shown in (Ref. 2) that a decrease in the 
decay time is caused by increased plasma diffusion across the force lines 

of the magnetic field. In its turn, the diffusion is caused by an insta- 
bility produced in the field of an ion cyclotron wave having a large ampli- 
tude. 


It may thus be assumed that under the conditions of our experiment 
the anomalously large frequency of ion collisions in high frequency fields 
with supercritical strengths is caused by the more intense (as compared 
with Coulomb interaction) interaction of ions with electrons. The reason 
for this (and for anomalous diffusion) is "bunched" instability. The fact 
that absorption is "collision absorption" -- i.e., it formally satisfies 
equation (3) -- where f(x) in the form of (4) means in physical terms that 
the damping force (which is caused by an instability) of the ion motion /23 
directed toward the electrons is proportional to the relative velocity of 
ion and electron liquids in the field of an ion cyclotron wave. 


Comparison With the Theory of Ion 
Cyclotron Wave Stability 


We shall show that the values of the effective frequency of ion colli- 
sions in fields with supercritical strengths, computed either according to 
the halfwidth of the resonance absorption curve or according to the absolute 
power at the absorption maximum, as well as the plasma diffusion coefficient 
determined by the decay time, coincide in order of magnitude with the corres- 
ponding values computed for our case according to the theory of ion cyclotron 
wave stability (Ref. 8). 
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As follows from (Ref. 8), under the experimental conditions (Ne~ 
~5°10!2 cm3; H~ 4-105 a/m; Te ~ 1 ev) for longitudinal (electrostatic) 


high frequency fluctuations excited by an ion bundle and propagated almost 
across a constant magnetic field (cos* O< MŒ), we obtain 


Reo ~ Imo ~ (wei), (6) 


where We is the cyclotron electron frequency. 


The effective frequency of ion collisions may be regarded as the in- 
verse of the time during which an ion bundle with an initial energy per 


unit of volume N,; Muc (u -- the velocity of the relative motion of the ion 


and electron components in the field of an ion cyclotron wave almost equals 
the velocity of an ion liquid) excites the fluctuations (6) and transmits 
all of the energy of the ordered motion to the electron gas. On the other 
2 
mu 


hand, the electrons obtain the energy Ni 7 


from the ions during the ex- 


l 
citation time of the fluctuations (wti) /2, The braking time of an ion 
bundle is thus 


1 
iee 0) 


Substituting numerical values in (7), we obtain 107 sec for t, which coin- 
cides in order of magnitude with Yate ~ 0.51076 sec which was determined 


experimentally. 


Let us determine the diffusion caused by instability, employing the /24 
theory of nonuniform plasma stability (Ref. 14). The diffusion coefficient 
may be written as follows 


D~V?t, (8) 


where V is the plasma pulsation velocity; t = + _. the characteristic time 


of correlation disappearance. The increasing increment œŒ fluctuations (6) 
must be used as v. The pulsation amplitude may be determined by the condi- 
tion of balance between two processes -- one of which leads to an increase 
in the pulsation amplitude due to the development of instability, and the 
other leads to contraction of this amplitude due to nonlinear processes 
leading to oscillation damping. As a result, we obtain 


V=va,, (9) 


where Aj is the oscillation wavelength in the radial direction. Then the 
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diffusion coefficient is 
D ~x}, (10) 
or, since ku y v, where k is the radial wave number, we have 
D~ (w) Ë, (11) 


Under the experimental conditions, in the case of jo = 8 a/m u ~ 
~10ê cm/sec and D ~ 2.10" cm?/sec, which coincides in order of magnitude 
with D~ 5°10* cm*/sec, determined for the same conditions according to 
the plasma decay time in (Ref. 2). 


Conclusions 


l. The absorption of high frequency energy, which was observed in 
our experiments, close to the ion cyclotron resonance both for small and 
large amplitudes of the high frequency field can be formally described by 
relationships (3) and (4), which were derived for the case of "collision" 
absorption. 


2. Absorption is caused by the interaction between plasma ions and 
electrons. In fields whose strength is less than the critical strength, 
it is caused by pair Coulomb collisions. In fields whose strength is 
greater than the critical strength, the effective frequency measured 
experimentally of ion collisions is considerably greater than the fre- 
quency of pair Coulomb collisions. In this case, the anomalously large /25 
absorption is apparently caused by braking of the ordered motion of ions 
in the field of the ion cyclotron wave by high frequency longitudinal 
oscillations excited by an ion bundle. 


3. The value obtained experimentally for the effective frequency 
of ion collisions, and also the diffusion coefficient determined according 
to the plasma decay time, coincide in order of magnitude with the corres- 
ponding values calculated according to the theory of ion cyclotron wave 
stability. 


4. On the basis of the results obtained, it is natural to pose the 
question of "turbulent" plasma heating (Ref. 15) in ion cyclotron resonance 
-- i.e., brief excitation in the plasma of an ion cyclotron wave having a 
large amplitude, with subsequent thermalization of the ordered motion of 
plasma particles in the field of this wave by the high frequency longi- 
tudinal oscillations excited by an ion bundle passing through the electron 
gas. The heating period must be quite small, so that anomalous diffusion 
caused by plasma instability does not produce significant losses in plasma 
particles. 
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INVESTIGATION OF CONDITIONS PRODUCING A DENSE PLASMA IN A 
METALLIC CHAMBER AND ITS HIGH FREQUENCY HEATING 


O. M. Shvets, S. S. Ovchinnikov, V. F. Tarasenko, 
L. V. Brzhechko, 0. S. Pavlichenko, V. T. Tolok 


It was shown in (Ref. 1) that powerful, high frequency oscillators may 
be employed to produce a dense plasma in a metallic chamber. The present 
experiments represent a continuation of this study by investigating the 
conditions producing a dense plasma in a metallic chamber and by determining 
the main plasma parameters. 


The utilization of a metallic chamber as the operational vacuum body 
and the method of feeding the system from the hf oscillator have definite 
advantages and unusual features. 


1. The system has a low input impedance, which does not require high 
voltages when large, high frequency powers are introduced, and avoids 
several technical difficulties which accompany, as an example, the spatially 
periodic circuit proposed by Stix (Ref. 2). 


2. A good connection between the feed electrodes and the plasma is 
provided. A dense plasma with a cylindrical form is produced between the 
central electrodes, and has no direct contact with the wall of the discharge 
chamber. 


Investigation of the Conditions Producing a 


Dense Plasma in a Metallic Chamber 


The "Vikhr'" device (Figure 1), on which the study was performed, con- 
sists of a copper chamber 1 with a wall thickness of 2.5 mm, an inner diameter 
of 125 m, and a length of 2000 mm. It is placed in a magnetic field which /27 
can be regulated continuously between 0 -2°105 a/m. The configuration of the 
magnetic field may be varied, depending upon the method chosen to handle 
the plasma produced. Since we were interested in plasma heating by high 
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frequency fields, the magnetic field configuration was chosen in accordance 
with the conditions which were requisite for generation and absorption of 
ion cyclotron waves. The region of the "magnetic beach" was located in 

the center of the selenoid 11 producing the magnetic field. The residual 
gas pressure in the system did not exceed 1.33-10-" n/m?. The rods 3 are 
introduced axially into the vacuum chamber through the porcelain insula- 
tors 2. Aluminum electrodes 4 with a diameter of 50 mm and a iength of 

70 mm are located at the ends of these rods. The distance between the 
electrodes is 1000 mm. The feed rods 3 through the coaxial cables 5 and 
the capacitance C are connected with the coupling coil 6 of the hf oscil- 
lator 7. The oscillator power is on the order of 100 kw; the operating 
frequency is 1.82-10® cps. The ends of the chamber are closed with glass 
discs 8. The operational gas from the flask 9 is admitted into the chamber 
by the valve 10 through an opening in the glass disc 8. 


With the system employed for switching on the hf oscillator, the feed 
electrodes acquire a negative potential, due to the rectifying properties 
of the plasma; this negative potential can reach several kilovolts with 
respect to the chamber. This potential creates the condition for oscilla- 
tion of electrons along the magnetic field force lines between electrodes, 
similarly to Penning discharge. The oscillating electrons effectively /28 
ionize the operational gas and provide for a high degree of ionization up 
to a chamber pressure of 1.33+10-? n/m*. A core of dense plasma is formed 
between the electrodes along the chamber axis; the diameter of this plasma 
is determined by the electrode diameter. Due to the presence of a constant 
potential, the peripheral plasma begins to rotate according to the law of 
plasma behavior in crossed electric fields and magnetic fields. Ina few 
microseconds the applied high frequency voltage produces a discharge in /30 
the operational area, which is accompanied by a voltage decrease at the 
electrodes due to an increase in the oscillator loading (Figure 2, a). 
The central electrode acquires a negative potential, which is retained for 
the entire period of time that the dense plasma exists (Figure 2, b). 
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Figure 2 


Lines of admixtures from the chamber wall materials and from the 
electrodes -- CuII, AIII -- and a weak line CII (Figure 3) were observed 
in the discharge spectrum obtained by means of a ISP-51 spectrograph with 
a UF-85 camera (focal distance 1300 mm). 


When the outline of the line Hg was measured, it was found that under 


our experimental conditions there was an apparant Stark widening by the 
micropoles of the plasma. This made it possible to determine the charged 
particle density, by comparing the contour observed experimentally with 
the theoretical contour computed on the basis of the theory advanced by 
Kolb, Grim, and Shen (Ref. 3). For purposes of comparison, Figure 4 shows 
the observed contour I and the theoretical Stark contour II, computed for 
n = 2.0101} cm3. The figure also plots the Gaussian contour III with 
the halfwidth equalling the experimental value of the halfwidth 0.7 Be 
Close to the maximum, the line broadening was caused by the Doppler 
mechanism, and the slopes of the line are due to broadening by the Stark 
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mechanism. 


Measurements of the plasma density made it possible to determine its 
comparatively weak dependence in discharge on the magnetic field strength, 
which confirmed the assumption regarding the nonresonance mechanism by 
which a plasma is produced in discharge. At a pressure of 0.4 n/n, the max- 


imum plasma density was 2-101" =: For an approximate determination of 
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Figure 4 


the plasma density distribution over the radius of the system, different 
sections of the plasma were focused on the spectrograph slit. The most 


dense plasma (density on the order of 10!* ss ) produces a filament with a 
diameter of 20 mm; at a distance of about 30 mm from the discharge axis, 
the plasma density decreased by more than one order of magnitude. 


The electron temperature was measured with respect to the intensities /32 
of singlet and triplet lines of helium; it was (4 - 5) 10° °K. 


The dependence of the plasma density change on time was determined 
by probing the plasma with an ultrahigh frequency signal at a wavelength 
of 3 cm (Figure 5, c) and 0.8 cm (Figure 5, d) through the glass openings in 
the chamber. The measurements showed that in an optimum regime there is 


a plasma with a density of 2 1013 4 in the apparatus for 3.6 milliseconds, 


and with a density of > 1012 =} for 17 milliseconds. The form of the hf 


pulse is shown in Figure 5, a. The maximum plasma density was determined 
according to measurements of the relative change in the intensity of the 
spectral line Hg with time. The intensity of the spectral line He (the 


photomultiplier signal, Figure 5, b), if the main excitation mechanism is 
electron collision, is 


Av apAgattgn, < ¥54(0) ) 
ee Ago t Ag eek: 


where v,, is the frequency corresponding to the line Ha 3 Aik -- probability 


of spontaneous transition from the i level to the k level; 


( voa (0) ) =f vos (0) fe (0) do; 
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5, (v) -- excitation cross section (by electron collision) of the n = 4 


level. f£f,(v) is the electron distribution function with respect to 


velocity; [v%,(v)] in the case of Maxwell distribution of electron veloci- 
ties is the function of electron temperature which -- as measurements of 
the dependence T(t) have shown -- changes very little during the dura- 
tion of the hf pulse. Thus, the intensity of the line Hg must change 


proportionally to the product ne’ In our case, the neutral gas from 


the cold section constantly enters the discharge column, since 
Rp1 


x % 1075 sec (Rp1 -~ radius of the plasma column; Y F thermal velocity 


n 
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of the neutrals). Therefore, the neutral density na % const and, conse- 


quently, the intensity of the line Hg is proportional to the electron 


density ne’ At the moment when n = 10!2 5, the intensity of the line 


H, decreases by a factor of 200 as compared with its maximum intensity -- 


B 


i.e., the maximum plasma density in our case is 2°10!" 4, which coin- 
cides with the contour of the line Hg determined according to Stark 


broadening. The plasma density determination based on the intensity of 


the line Hg» at the moment when a signal with a wavelength of 0.8 cm 


begins to pass through, coincides with the result of microwave measure- 
ments. 


N 
w 
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Thus, we have studied the conditions producing a dense plasma in a 
metallic chamber. A plasma with a density on the order of 10!" 4 
and an electron temperature (4 - 5) 10° °k was obtained in the experiments. 
The weak dependence of plasma density on the magnetic field strength points 
to a nonresonance mechanism by which the plasma is produced. The aim of 
our subsequent experiments was to investigate the heating of the plasma 
obtained by the generation of ion cyclotron waves. 


High Frequency Heating of a Dense Plasma 


in a Metallic Chamber 


The next stage of our investigation was to study the possible heating 
of ions by generating ion cyclotron waves in a dense plasma produced by a 
hf oscillator in a metallic discharge chamber. The "Vikhr'" device was 
modernized so that, when the ion cyclotron wave was generated at the ends 
of the coaxial close to the magnetic mirror, it was propagated into the 
center of the discharge chamber, where the magnetic field strength was de- 
creased to a value equalling the cyclotron value for protons, forming the 
region of the magnetic beach. When the wave was propagated along the system 
axis and approached the region of the magnetic beach, its velocity decreased, 
and it was damped, transmitting its energy to the plasma ions. In order 
to propagate the wave in the region of the magnetic beach, it is necessary 
that the magnetic field strength per wavelength change by several percents 
-~ i.e., a smooth change in the magnetic field strength over the length of 
the system is requisite. If the opposite is true, the wave will be reflected. 


The residual gas pressure in the system did not exceed 1.3°107* n/m?. 
The experiment was performed with a mixture of two gases -- hydrogen and 
helium or hydrogen and argon. The operational pressure was established 
by a stationary regime of the valve operation in the 0.1 - 0.8 n/m? 
range. At a pressure of 0.4 n/m?, the hf oscillator, operating at a 
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frequency of 1.82+10© cps with a power of 150 kw, produces a plasma with 
1 


a density on the order of 10!" n? With a magnetic field strength of 
H = 1.25 Hy A (H, a53 the cyclotron value of the magnetic field strength 
for protons) waves are generated in the plasma with A = vac = 6.6 cm [34 
ec 
(Agac 77 Wave length in a vacuum; n= ~y Fa i~  —- refractive index 
YA (ay 
of the medium; v, = a: eee Alfven velocity; H -- magnetic field; w -- 
Y4Tp 
operational frequency of the oscillator; Weg m7 cyclotron frequency of 
proton rotation; p -- mass density. 


This system has several advantages. The pulse which is initially 
transmitted to the ions is perpendicular to the vector of the outer mag- 
netic field strength, which facilitates the retention of ions in the cork- 
screw configuration of the magnetic field. Since the system has a low inped- 
ance, energy is readily introduced into the discharge chamber and is trans- 
mitted directly to the plasma ions. An increase in the density and diameter 
of the plasma does not make the conditions worse for wave generation. At 
the same time, the spatially periodic circuit for introducing hf power 
into the plasma, which was advanced by Stix (Ref. 2), loses any physical 
meaning with an increase in the plasma density and diameter. When a mix- 
ture of two gases is heated, or when it is necessary to heat the plasma 
electron component simultaneously with the ion component, this system makes 
it possible to introduce the power of two oscillators operating at different 
frequencies. 


As has been pointed out, an increase in the plasma density does not 
impede wave generation, since spatial periodicity is not given externally, 
but is established as a function of the plasma refractive index, and may 
be small (several centimeters). However, the possibility of periodicity 
is not excluded, if rings are placed endwise on the coaxial at different 
distances from each other; these rings will introduce a perturbation, 
creating a specific periodicity along the system axis. One unusual feature 
of the discharge is the fact that its nucleus, consisting of a plasma which 
is almost entirely ionized, is surrounded by a plasma having a low density 
and a neutral "housing". During generation and absorption of ion cyclotron 
waves, the dense plasma nucleus is heated, and the cold plasma surrounding 
it with a low density contributes, in all probability, to the suppression 
of channel instability. 


In the experiments described, the ion temperature of the plasma and 
the charged particle density, which were averaged over time, were deter- 
mined by means of optical methods. The change with time in the intensity 
of the spectral lines for hydrogen and admixtures was also studied. A 
ISP-51 spectrograph was placed in such a way that the central portion of /35 
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Figure 6 


the cylindrical plasma column, lying in the region of the magnetic beach, 
was focused on its input slit. Oscillograms showing the intensity of 

the spectral lines for hydrogen and helium or argon, added to the chamber 
in small proportions to the operational gas, showed that regular oscilla- 
tions are produced in the intensity of the spectral lines with a frequency 
on the order of 20 kc, close to the magnetic field strength corresponding 
to the cyclotron value for protons (Figure 6). These oscillations appear 
during wave generation in the plasma, and are caused by the eccentric ro- 
tation of the dense plasma filament as a whole with respect to the chamber 
axis, in accordance with the drift law in crossed, radial electric fields 
and axial magnetic fields. The direction and frequency of the filament rota- 
tion was determined by two photoelectron multipliers oriented towards the 
end of the chamber and placed at a radius of 3 cm from the system axis. 
One of them was shifted along the azimuth. The oscillation phase of the 
light intensity was thus changed. 


One interesting feature was discovered when the ion temperature was 
measured according to the Doppler broadening of the hydrogen and additional 
gas lines. The width of the hydrogen line depended comparatively little 
on the magnetic field strength, and the additional gas lines were broadened 
considerably when the magnetic field strength was close to the cyclotron 
value for protons. The ion temperature of this gas, determined for the /36 
optimum operational regime of the apparatus, amounted to 250 ev (2.5+10© °K). 
Measurements of the halfwidth of the line H, showed that the hydrogen atom 
temperature was below the temperature of thé additional gas, while there 
was considerable Stark broadening of the line contour, corresponding to a 


4, 


plasma density of 2°10! There was also Doppler broadening of the 


3° 
cm 

admixture lines (copper, aluminum, oxygen, carbon, nitrogen). When the 
temperaturesof different admixtures and additional gases were measured, 
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Figure 7 


no ion temperature dependence on their mass was detected. The heating 
of the additional gas ions was terminated when only one additional gas 
was admitted into the discharge chamber up to the previous operational 
pressure. One of the characteristic dependences of helium temperature 
(when it was added to the chamber in a small amount) on magnetic field 
strength, measured over the halfwidth of a helium line with a wavelength 
of 4921.93 A, is shown in Figure 7. The helium atoms acquire a maximum 
temperature at a magnetic field strength which is close to the cyclotron 
value for protons. A certain temperature increase is observed as a mag- 
netic field strength is approached which equals the double cyclotron value 
for hydrogen ions. 


If it is assumed that the halfwidth of the H, line is determined by 


B 
Doppler broadening, the temperature of neutral hydrogen is considerably 
lower than the plasma ion temperature. However, under our conditions 

(in the case of Te = [4 - 5] 10° °K), the lifetime of neutral hydrogen in 
the plasma, with respect to the ionization process, was small as compared 
with the time of Coulomb collisions. Therefore, neutral hydrogen cannot 
acquire energy equalling the ion energy. Since the frequency of hydrogen 
atom collisions with electrons is greater under these conditions than the 
frequency of collisions with ions, the Stark mechanism is the predominant 
mechanism leading to the broadening of the line of the residual neutral 
hydrogen. 


In an optimum operational regime in a mixture of two gases -- hydrogen 
and argon, the dependence of the argon ion temperature on the magnetic 
field strength is resonant in nature with a maximum close to the magnetic 
field strength corresponding to the cyclotron value for protons. The 
maximum temperature of argon ions in the experiments was 2.5*10© °K, and 
for electrons -- 5-105 °K. In order to determine the distribution of the 
plasma ion temperature over the radius of the system, different sections 
of the plasma were focused on the spectrograph slit. The hottest plasma was 
located in the center of the filament. The temperature rapidly decreased /37 
over the radius (approximately 5 times greater along the axis than on a 
radius of about 4 cm). Figure 8 illustrates the temperature dependence of 
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Figure 8 


plasma on the applied high frequency voltage. The increase in the mag- 
netic field strength and the high frequency voltage -- i.e., the hf power 
introduced -- made it possible to obtain a higher plasma temperature. 


The mechanism by which the proton energy is transmitted to the addi- 
tional gas has still not been definitely clarified. Since the gasokinetic 
pressure in these experiments may exceed the magnetic pressure, the possi- 
bility is not excluded that centrifugal instability may be produced, which 
can lead to transmission of proton energy to the additional gas and can 
lead to its heating. 


In certain operational regimes of the apparatus, we observed genera- 
tion of a rapid magnetosound wave at a magnetic field strength which was 
less than the cyclotron value for protons. However, conditions were not 
favorable for studying it at the existing oscillator frequency (1.8210 


cps). 


Thus, in all probability, these experiments illustrate the feasibility 
of high frequency heating of a dense plasma consisting of two types of 
ions by resonance generation of ion cyclotron waves for one type of ions. 
The mechanism by which energy is transmitted from one type of ions (pro- 
tons) to other ions (helium, argon, admixture) requires further study. 
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HIGH FREQUENCY PLASMA HEATING 
K. N. Stepanov 


As is well known, a rapid decrease in the frequency of Coulomb colli- 
sions of plasma particles leads to a decrease in the rate of ohmic plasma 
heating both by constant fields and by high frequency electric fields when 
there is a temperature increase. Therefore, great hopes have been expressed 
(Ref. 1, 2) for achieving plasma thermonuclear temperatures by utilizing 
different effects of collisionless energy absorption of high frequency 
fields by a plasma. In principle, this amounts to Cherenkov, or to cyclo- 
tron, absorption and radiation of waves by electrons and ions. 


However, collisionless energy absorption by a plasma leads to great 
distortion of the velocity distribution function of electrons and ions and 
to attenuation, or even complete discontinuance, of absorption. It also 
leads to every type of plasma instability. If collisionless plasma heating 
is effectuated by weak electric fields, so that Coulomb collisions provide 
a Maxwell distribution function, the heating rate is the same as for ohmic 
heating. The time required to heat the plasma up to thermonuclear tempera- 
tures with weak fields is less than the time required for containing the 
plasma in a thermonuclear reactor, utilizing the reaction D + D, witha 
positive energy output. Therefore, under the condition of producing a 
stable, plasma configuration, it is primarily possible to achieve thermo- 
nuclear plasma temperatures with slow heating by weak, constant fields or 
by high frequency fields, which in themselves do not lead to strong plasma /39 
instability. 


On the other hand, for rapid plasma heating it is very tempting to 
employ strong electric fields, under whose influence electrons or ions 
in the plasma acquire a large directed velocity (Ref. 3 - 6), as well as 
strong electron bundles (Ref. 7) or the collision of plasma clusters 
(Ref. 3) ("turbulent" methods of plasma heating). Due to the development 
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of bunched instability in the systems, there is a rapid increase in the 
energy of high frequency oscillations leading to braking of the acceler- 
ated particles, to an energy increase with respect to motion ("tempera- 
ture"), and to an energy exchange between different plasma components. 


This article presents a brief survey of the collisionless, high fre- 
quency methods of plasma heating, and it also compares them with the 
ohmic heating method. 


Ohmic Plasma Heating 


When investigating the processes of plasma heating, we shall disre- 
gard energy losses due to cyclotron radiation of electrons, which is ab- 
sorbed by a plasma if its dimensions are fairly large, as well as losses 
due to plasma thermal conductivity on the chamber walls -~ which is also 
a surface effect. In addition, we shall assume that the residual gas 
pressure is small, so that we can disregard energy losses due to over- 
loading and radiation of admixed atoms excited by electron collision. 
Under these conditions, the energy losses by the plasma are determined 
only by braking radiation of the electrons. 

When the ohmic method is employed to heat a plasma by a "constant" 
electric field with the strength E, the electron energy increase is deter- 
mined by the following equation 


T= R, (1) 


where w = EA is the mean electron energy; T, ~- electron temperatures; 


e 
Ng Q4. = sE? (2) 
-- Joule heat liberated per unit volume per unit of time; ng -- electron 
2 — nan 
density; o = maa Ce plasma conductivity; eo Van -- time of 
De 4Y QangetA 
electron mean free path; A -- Coulomb logarithm. The intensity of elec- 
tron braking radiation (Ref. 8) is 
327/27 „e° 
Q_ = 32V Tee (3) 


3 V am, fic'm, K 


If we take into account the transmission of energy to ions, then we 
must replace Q, by Ż Q. in the case of T, ,% T;. This effect, as well as 


other effects which do not change in order of magnitude, is not taken 
into account from this point on. 
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Under the influence of the electric field E, the electrons acquire 
a directed velocity 


sE 
u = —. 
eng 
engVe 
If the field strength E exceeds the critical value (E> » where 
0 


7 
u= E is the electron thermal velocity), then all of the electrons 
e 


are carried along by the electric field in a continuous acceleration regime. 
The motion of the electron gas with respect to the ions leads to the phe- 
nomenon of bunched instability which is related to the buildup of longi- 
tundinal high frequency oscillations in the plasma (Ref. 9). The inverse 
influence of plasma oscillations on the electron motion leads to electron 
braking -- i.e., to anomalous plasma resistance, and also to increased 
radiation of radio waves, which considerably increases the plasma thermal 
radiation (Ref. 10, 11). 


Plasma bunched instability can arise in the case of E << Eorp» i.e., 
in the case of u << vg. For example, in a very non-isothermic plasma, 


in the case of Tẹ >> T4 (T; -- ion temperature) the electrons build up 


sound oscillations (Ref. 12), if the velocity u exceeds the speed of sound 
T, 
v= z - In the case of T, S Ty, the electrons build up longitudinal 


ion cyclotron oscillations, if u > 10,5 (vy = nee thermal ion velo- 
i 
city) (Ref. 13). However, it may be expected that -- since only a small 
group of resonance electrons participates in the oscillation buildup -- 
the formation of a "plateau" in the electron distribution function will /41 


lead to a decrease in the increasing increment ( Ref. 14 - 18), and nonlinear 
effects will lead to stabilization of these oscillations (whose amplitude 
will be small) (Ref. 16 - 18). 


Expression (2) may be employed only in the case of E << Ecore Assuming 
that E = Ecer, where a << 1, we obtain 
aT, 


Q~ (4) 


e 


Let us substitute expressions (3) and (4) in (1). We then have 


—9 aoas 
dw = 5 1o-"V 7) n kevfsec. (5) 
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where Te -- in kev, ng -- in cm3. The heating process terminates, 


iu. O in the case of | 


dt 
T ~ 5. 105a? kev 
max . (6) | 


It thus follows that even in the case of heating with weak fields (a = 
= 107! - 107%) it is possible to achieve thermonuclear temperatures 


(T ~a 50 kev). 
If Q. << Q,, then -- disregarding Q_ as compared with e€ Q, in (1) 


and taking into account (4) -- we obtain 


eee ese (7) 


where Tg and tg -- are the initial values of T, and Tẹ» The influence of 


braking radiation on the heating process is significant in the case of 
Te ~ Tax? i.e., in the case of 


5- 105 \3/ 
t~ a(r) 2. (8) 


For example, in the case of Tp ~ 100 ev, no a 1015 cm-3, to% 2:10 8 
sec and q % (1/30), the temperature T, = 50 kev is achieved during the 


time t ~ 0.2 sec 


and the electric field strength changes between 0.3 - 3:107" v/cm. The 
energy exchange between ions and electrons takes place during the time 


te 0.5 sec 


Mi 
Tie ~ Te m, , (10) 


so that the separation between the electron temperature and the ion tem- /42 


perature is small. 


The heating time (9) is not large from the point of view of producing 
thermonuclear reactors with a positive energy balance. The time the plasma 
is contained in such a reactor, when employing the reaction D + D, must be 
greater than (Ref. 19) 


i*m sec. (11) 
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For the example under consideration (ng % 1015 cm~ 3) t* ~ 10 sec, 
which exceeds the heating time of the ion component by a factor of 20. 


Cherenkov Ion Heating 


In the case of high frequency plasma heating employing the method 
advanced in (Ref. 20), oscillations of an axial magnetic field H, which 
is produced by azimuthal electric currents, lead to a variable azimuthal 
electric field E, in a plasma cylinder located in a strong longitudinal 
field Ho. Radial, drift oscillations of the plasma arise, due to particle 
drift in crossed fields Eo and Hg, which leads to density oscillations -- 


i.e., to the appearance of sound (more precisely, magnetosound) waves. 
(It is assumed that the wave frequency w is considerably less than the ion 


eH 
gyrofrequency wj = — and tne wave length is considerably less than the 


i 
v 

Larmor radius of ions pį = i with thermal velocity). 
Wi 


Ions having the velocity v, along Hj, which is close to the phase 

velocity of a wave Vo Ea vigorously interact with the field Ey: If 
It 

the wave phase velocity is on the order of the ion thermal velocity, the 
number of resonance ions is large, and there is strong wave absorption. 
(As is known, in the absence of a magnetic field sound oscillations in a 
plasma, in the case of Te < Ty, cannot be propagated in general, due to 
strong Cherenkov absorption by ions (Ref. 21); in the case of Hg #0 and 
Te < Ty, magnetosound oscillations are also damped during one period, if 
Vo a Vi [Ref. 22]). 


Let us investigate Cherenkov absorption by plasma ions of the energy 
of an electromagnetic field produced by azimuthal electric currents, /43 
which take the form of a moving wave and flow into the coil placed on a 
plasma cylinder having the radius a: 


je = ja cos (k ız — wt yè (r — a). (12) 


If ky a< 1, then the current (12) produces a variable longitudinal 
magnetic field with the strength 


H; =H cos (kız — vt), (12") 
where H = mg The strength of the aximuthal electric field is 
E EE rsin (kız — ot 
P= — aa? 1z — of), 
kyc 


where ny = uF is the longitudinal refractive index. 
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In the absence of damping, the value of the energy flux, averaged 
over time, in the plasma per unit length 


equals zero. It is apparent that, in the presence of damping, the energy 
flux in the plasma equals, in order of magnitude, 


aC fjg % 
San A 


[aa l (13) 


where x is the correction to the "transverse" refractive index n,, caused 
by the Cherenkov oscillation absorption by plasma ions. 


The dispersion equation for a magnetosound wave has the form (Ref. 22) 


2 2 Ë 1 + in? o 
nn = v2, + a7} (14) 
where the coefficient n ~ 1, if Vo ~ vy. The component ~ i takes into 


account Cherenkov wave damping in an ion gas. 


Let us investigate a plasma with a small gasokinetic pressure 
8nnoT ie H 
at <1} In this case, the Alfvén velocity Va=——%— is consider- 
0 V innom; 
ably greater than the ion thermal velocity v; and the speed of sound V,. 


Since ni~ => in the zero approximation it follows from equation (14) that /44 


2 
v 
n| = inj. In the following approximation we obtain: "1 = in (14 bint #1), 
i.e. 
x 8anoT; 
[ni | H? (15) 


Taking into account expressions (13) and (15), we obtain the following 


formula for the mean increase in plasma ion energy oe = aay (Ref. 23) 


dw H 
a~ (a) oTi (16) 


In actuality, the energy is absorbed by resonance ions with v, % 


x Vres = a These ions may be regarded as magnetic dipoles with the mag- 
1l 


netic moment KE mo _ Ti which moves in the wave magnetic field with a 
2Hy Ho 
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strength which slowly changes in time and space. Within the frame of 
reference in which the wave is at rest, the equation of dipole motion 


aH, 


Oz 
leads to the law of energy conservation 
m;z? 
aT + pz = const. (18) 
Ions with velocities vy in the Ki -Av<vy a + Av range, where 


Av = = (19) 
i 


are trapped in the potential well and effectively interact with the wave. 
It is apparent that the number of resonance particles per unit volume is 


Av H 20 
E (20) 


Due to the shift of ions under the influence of the force -ų aa 


> 
the distribution function in the resonance region is distorted during the 


time Thonlin.? during which a trapped particle covers a distance on the 


order of the potential well width, i.e., 


1 
TrTel~ kh, Av’ (21) 


Relaxation of the distribution function due to collisions in a narrow re- 


gion of the width’ Av close to vj = YE occurs during the time period 
Av\2 (22) 
*rer~ (5) ’ 
where _ 
3m; T 
u= =. 
4 V 2k nget A 


It is apparent that the ion distribution function will be distorted by an 


insignificant amount if Thonlin. >? Trey: The critical value of the varia- 


ble magnetic field strength H = H,,, at which the distribution function dis- 


tortion becomes significant, is determined from the condition Thonlin. ~% 
œ% Trel» Taking into account (19), (21), and (22), we obtain (Ref. 23) 
H 
ke~ (ony (23) 
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Expression (23) for H,, can be obtained by another method. An in- 


crease in the perpendicular velocity component (or a change in the mag- 
netic moment) in the presence of a variable magnetic field is determined 
in the drift approximation from the following equation 


- Yo oH, 
Ua = Hok, z” (24) 


The time of nonlinear distortion of the distribution function with 
respect to the velocity v} under the influence of the field Ey -- i.e., 


under the influence of the force V nie in (24) -- is determined according 
z 
to the relationship 


1 


T in ~z 
nonlin k, Ao, (25) 
where 
E ee U~ Ui 
o ki (26) 
If the ions obtain a velocity increase. Av,, then the relaxation [46 
time of the velocity distribution function is determined by 
Ao, \2 
Trel veil | é (27) 


The critical value of the magnetic field strength H = Hoer» at which 
great distortion of the ion distribution function f(v, vj) occurs close to 
w : : ‘ gee r 
Vy = KP which is determined from the condition Thonlin. ~Trel’ coincides 


with the value of (23) in the case of k Mi It is also apparent that the 


distortion of the velocity distribution function of the particles close to 


WO, be ees ; ; : 
vi = ky is not significant, if the energy acquired by resonance particles 
Ww 


having the velocity vj (in the a Av); < vy < R + Av, interval), during 


the time between two collisions 1; 


~ 3, 
dw dw i H 
oe ~~ |>] oTr 
di x dit ` Bo, @ (2) ifi, (28) 


is small as compared with their thermal energy T;. The critical field 


no VTi,» coincides with (23). 


strength H, determined from condition t T 


res 


i 
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During heating with fields H = oH, under the condition om Wi (it is 
apparent that the fulfillment of this condition requires that the fre- 


quency be increased proportionally to Y T; as the temperature increases) 


ac SIL ll CT (29) 
dt ~ a, (we) % (o) t NT: 


where Tg, Tg and wo are the initial values of Te» Te and w. 


Assuming that w~ T; and integrating equation (29), we obtain the 
following expression for the ion temperature (Ref. 23) 


a? */ss a?t wh 
T= Tofi t =| =T,[1+ a] : (30) 


To (oto) G To (oto) / 


Plasma heating by means of ion Cherenkov resonance may be intensified, /47 
if -- instead of one wave -- several waves are employed with phase veloci- 
ties differing by 2 - 3 Av, so that all particles in the -vi < VI < Vi 


range may be resonance particles. In order to do this, it is eee to 


produce a wide wave packet with a phase velocity scattering of iG a Vee 


It is apparent that the total number of such waves is iy a It is also 
evident that in this case the critical strength of the variable magnetic 


field responsible for great distortion of the ion distribution function over 


> 
the entire |v] < vy range is 


U; aan oF H 
Hor, tot.~ to Her~ V Heo ~ (ony i (31) 


In the case of heating by the field H, ~ aH (A ~ a?Hor) in the 


cr. tot. 
case of a wide wave packet, we have 


dw v; 7\ oT. ~ alla 
dt Bo i. ers (32) 
Integrating expression (22), we obtain 
T= T(1+ 2)" (33) 


Formulas (32) and (33), which determine the ion temperature during the 
heating of plasma ions under optimum conditions, when practically all the 
plasma ions participate in absorption of the high frequency field energy, 
are similar to formulas (1), (4), and (7), which determine Joule electron 
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heating by subcritical "constant" fields. 


In the case of a = 0.01, nọ% 10!5 1/cm3, To = 100 ev, Hg 10°C, 
and wg% 107 sec}, we obtain: To% 10-6 sec, T; = 50 kev during the time 
ta% 1 sec, while in the case of t = 0 H, ~500G, Ey © 50 v/cm at the end 


of heating w ~ 20 wo ~ 2*108 sec”! ~ 0.2uy and H, ~ 20G, Ey ~ 40 v/em 
l 
T 
Ti 


is the fact that energy is transmitted directly to the plasma ion com- 
ponent. 


t=0 barely changes). One advantage of this method 


In the case of a very non-isothermal plasma (Te >> Ty), as was shown 


in (Ref. 24), a resonance relationship between the outer circuit and the /48 
plasma can exist in the case of Ve% Vs. In this case for Te < 10T; ion 


absorption is still significant (on the order of electron absorption), and 
resonance Vo = Vs can exist. The energy absorbed at the maximum increases 


by 10 - 100 times as compared with (16) (Ref. 23). 
An expression for the high frequency power absorbed by a nonuniform 


plasma cylinder in the case of ah v v4 was obtained in (Ref. 25). 


As of the present, the effect of Cherenkov absorption of a magneto- 
sound wave by plasma ions has not been studied experimentally. 


Cherenkov Electron Heating By the Field of a 


Magnetosound Wave (Helicons) 


A rapid magnetosound wave can be propagated in a plasma with a large 


density (22>, , where o= fem -- Langmuir frequency, We = eng 
Me Pac 
e 


electron gyrofrequency); the refractive index of this wave in the wj<<w<<w, 
frequency region is determined by the expression 


Oo e a 
"Te, Va a (34) 


(in this frequency region, rapid magnetosound waves are called "whistling 
atmospherics", or simply "atmospherics", "whistles", and "spiral waves"). 


Cherenkov whistle absorption by plasma electrons is weak not only 
Ww 


kil >> ve, when the damping coefficient is exponentially small (Ref. 26, 


for 
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27), but also for — v In the last case (Ref. 28, 29), we have 


ki < e’ 
a Fie 
ao w, (35) 


e 


Since the whistle field readily penetrates a dense plasma, these 
whistles may be employed to heat the plasma electron component, while 
electrons located within the plasma cylinder will be heated up due to /49 
the comparatively small absorption coefficient (35). 


If ky a v1, then -- substituting (35) in equation (14) -- we obtain 


~\2 Ve 
dw H T o, 8noT, 
n C lo a) (36) 
n Oe 8xnoT, 1/2 ` . S 
Since the factor oe differs very little from unity in order 
0 


of magnitude, this formula coincides with formula (16) (substituting Te 
by T; i). However, since w >> wy in the case under consideration, and since 
formula (16) was obtained for the frequencies w << wy, for the same values 


of H electron heating by the whistle field takes place much more rapidly 
than Cherenkov heating of ions by a low frequency field. 


The change in the electron distribution function with respect to vy 
by the quantity 


it 
W 
Av, ~ Ue Hy Ry, (37) 


is insignificant, if the time of nonlinear distortion of the distribution 


function Tyonlin, ~% is a little greater than the relaxation time 


E&L.. 
k Av 


2 
v 
Trel “re(—4) «e The critical strength of a variable magnetic field 
v 
e 
H = Hor» at which the collisions cannot equalize the distribution function, 


is determined from the condition T 


nonlin. % trel? 
Ho [EtA Ho ve AnnoT .\*/ 
cs wae.) (eye A] oe 


This expression for Hey is obtained under the condition that resonance 


‘o 
electrons ee AO < A + Ava), which acquire the following energy 
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per unit of time 


dw 
dt 


~ a, Ye 
~ dw . Pe H T ®e 8xioT, 
res @ My Ay Of ele” a (39) 


If H = Hoy during the time ~% Te, acquire energy which is on the order of /50 


the thermal energy. 


During electron heating by whistle fields with H = aH.y, the energy 

acquired by one electron on the average per unit of time is 
dw aT, W HE a 

dt ~~ ; E : aT, (40) 


: He He \"* 
Assuming, for purposes of simplicity, that = mr} ~1, we obtain 
Oe Bano, 


a?t Ys 
to (wt) '/* 


(41) 


Under optimum conditions, when a wide wave packet is employed and when 
a large portion of electrons with a velocity of vyl < Ve are resonance 


electrons, electron heating by a field with Hp %aHer, tot, vo crHo is 
determined by the following expression 


dw aT, 
d n’ 


gnm, a 


T.=T(1+% 


These formulas describe Cherenkov electron heating by a whistle field under 
optimum conditions, if all the plasma electrons participate in energy absorp- 
tion. They are similar to formulas (4) and (7), which may be used to deter- 
mine the electron temperature increase during ohmic heating by fields with 
subcritical strength. 


The study (Ref. 30) was devoted to a theoretical investigation of 
Cherenkov electron heating by the electric field of a rapid magnetosound 
wave propagated in a nonuniform plasma cylinder. Cherenkov whistle ab- 
sorption was determined experimentally in (Ref. 31). 


Ion Cyclotron Resonance. 
Absorption of Alfvén Wave 


Ions are heated most effectively by the field of a high frequency wave 
under conditions of ion cyclotron resonance. If the frequency of a wave 
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propagated in a low-pressure plasma is close to the ion cyclotron fre- 
quency, then there will be a large number of particles having the velocity 
Vil, which is close to the particle resonance velocity /51 


om 


Pres ge (43) 


and effectively interacting with the wave field. Therefore, the wave 
energy absorption by ions will be large. Cyclotron damping of waves is 

the inverse of cyclotron radiation of waves by charged particles in a 

` magnetic field. This effect was first studied in (Ref. 32), where they 
investigated the damping of magnetohydrodynamic waves propagated along the 
Magnetic field. A study of the cyclotron absorption of waves having a 
frequency of w % wj was pursued in (Ref. 33) [see also the studies (Ref. 20, 
34, 35)]. 


As is well known, in the frequency region on the order of wy ina 
cold plasma (Ty = 0) there are two branches of oscillations corresponding 
to an Alfvén wave (which can only be propagated in the case of w < a;) 
and corresponding to a rapid magnetosound wave. As the Alfvén wave fre- 
quency approaches wi, its refractive index and cyclotron damping coefficient 
increase. In the case of wy - w< kjj vy, the propagation of this wave by 
strong damping is impossible: Re k~ Im k~ 1/6,;, where 65; is the depth to 


which the field penetrates the plasma (Ref. 28, 36, 37) 
ana L, (44) 


A magnetosound wave is absorbed eee in the case of |w -—al< ki Vi- 


Let us first examine cyclotron heating of a plasma cylinder by an 
Alfvén wave which is strongly damped. The azimuthal currents (12) excite 
this wave if kı 6; ~ 1. Also assuming that ky a~ 1 and wy — w< ky Vis 


we obtained the following expression from formula (13) for the mean energy 
acquired by one ion per unit of time 


~ ~ \2 H? 
dw He H oT 
dt ~~ Brg  \Ha/ ?* 1 Bun gT, * (45) 


It thus follows that the electromagnetic field energy which is accumulated 


2 
in the plasma i is absorbed during a period of time on the order of 


vi 


In addition, a comparison of (45) and (16) shows that, for the same /52 
variable magnetic field strength, the energy absorbed by the plasma in 
the case of w uw; is considerably greater than in the case of Cherenkov 


45 


resonance. In the first place, this is due to the difference in the fre- 


quencies oe) w and, in the second place, it is due to a large factor 


dt 
H2 
——— in expression (32). On the other hand, for the same values of T 
8mnoTi 


and k, the electric field strength in the case of cyclotron resonance is 
Ws 
times greater than in the case of Cherenkov resonance (cher -- 


“Cher. 
the wave frequency under Cherenkov resonance conditions). 


In order to determine the number of resonance ions responsible for 
energy absorption during cyclotron resonance, let us investigate the parti- 
cle motion in the field of a flat cyclotron wave: 


H, = H sin (kız — vt), H, = —H cos (k yz — vt); 


H H. (46) 
E, = ——cos(k,z— ot), Ey = — — sin (k ız — wt), 
ny ny 
k ye E x . ° 
where fy ar peq7 aed is the refractive index. 
AV 9/9 — 


In the absence of a wave, the ions move along a spiral: 
ne Vo = (oy cos (wit + Po), — 0% sin (wit + Po)» vi) 
When there is a weak field, the velocity perturbation is determined 


according to the equations of motion 


~ 


. eH 
uU, = — — 7° 

z mn, cos (kız wf) -F Wy} 

. eH . 

U ga — — WU" 

y nn sin (kız wt) 0x5 

ox 
-e e ev, H 
om, 0 0 g = 
z= ne (02H, vg Hs) = — —— cos ®, 


‘ 


where Ọ = kyz + (v: — o) t + Yo 


Assuming that u = Va + ivy, from the first two equations we obtain 5 
F eat _ eH i (k yz — wf) 
u + ivu = mai i e è 
Th = -iwjt a 
us, assuming that u = uge » we obtain 
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ao 


~ 


eH 


(> — P.) 
mi” y = g 


U = — 


The third equation can be written as follows 


iS 
5 ULAR (47) 
P + agp COSP = 0. 


Let us first investigate nonlinear distortion of the distribution func- 


tion close to v|j = Ves? caused by a variable magnetic field. The law of 


energy conservation follows from (47) 
ev? Hk 


l 
me sin ® = const. 


1 s 
zP + 


The time of nonlinear distortion of the distribution function close to 
Vil = Vreg equals the potential well flight time of trapped particles 


on 
vH, 
mk yc 


having the velocity ~ Avj y 
Tnonlin aoe 
Ry Avy Š 


We obtain the following expression from the condition T li T 
Avj \2 nonlin. ~ rel~ 
~ aa for the critical value of the magnetic field strength 


Vi 
Ho [BTN 
Ber~ cate H j’ 


(48) 


This expression may also be obtained from the formula for the nonlinear 
decrement of cyclotron wave damping, which is determined on the basis 


of the quasilinear theory (Ref. 38) in the case of V es >> vy, if we set 


s~ Vy and kı sl- in formula (25) of the study (Ref. 38). 


vV. 
re . 
õi 


However, the influence of the accelerating field È on the distribu- /54 
tion function change in a plane perpendicular to Hp is considerably 


stronger than the influence of a variable magnetic field on the distribu- 
tion function change with respect to velocity along Hg. 


Let us determine the time of nonlinear distortion of the distribution 
> 
function by the field E: 
T ra 
nonlin ~ yao,’ 


where 
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= = ‘fs 
2 
mE Vi fo 49 
T Bene ace a 


Av, \2 
P L : ~ 
Equating T onlin: and Trel ~“ (= in the case of Ha Hoy» we obtain 
i 
Hw Hy [Eno TN 
cr (#t,)*/a H? i (50) 


We obtain the same expression by assuming that only a group of resonance 


particles with Vee > AVL < Vil < Vreg + AVL participates in energy absorp- 


tion. Resonance particles acquire the following energy per unit of time 


a 3 2 
dw H /s a H? le 
d es \Ho AE ° 


In the case of H ~ Hery» during the time a Tį these particles collect 
the energy ~ T;. The critical field strength (50) is considerably less 


than (48). This means that the influence of nonlinear effects caused by 
the electric field is manifested for smaller field strengths than is the 
case for the influence of a variable magnetic field. 


In the case of H = OH.» we obtain the following expression for ae 


dw aw, H 3 "e a°To Hg i Te 5 
a o jn T) Sosa T) (T) ” ee 
Thus, we have 
H? 1/7/15 aĉt T 
Ty = Tfi Eeen To = 5 1/3 lest] l me rfi + To (oto) f (52) 


Let us now present a numerical example. Let us set ng v 1013 cm3, /55 
Ho ~ 5°10°G, w; % 5:10’ sec™l! and a% 1. Then in the case of T; = 10 ev, 
dw 


we obtain: 1; % 31076 sec, H ~ 2G and at ~ 300 kev/sec; if Tá 100 ev, 


then tį; 107" sec, H~ 0.46 and oe ~% 30 kev/sec; if T; ~ 1 kev, then 


Ti % 3-10-3 sec, Fa 0.1G and @ 4 3 kev/sec. 


re 


In the case of Tp % 100 ev and Ha Hoy» the temperature T; v 10 kev 


is achieved during the time t ~% 3 sec. 
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Cyclotron heating by fields with subcritical strengths may be in- 
tensified, if several waves (wave packet) are employed having phase 
velocities which differ with respect to 2 ~ 3 Av]. The number of such 


waves is ae and the critical strength of the total magnetic field is 
Av 


H ~ Y Hory» where Hep is determined according to (50). In the 


cr. tot 
case of heating by the field H, ~ oH... tot’ 


tain (32), and for Tį we obtain (33). 


in this case for © ve ob- 


Cyclotron wave absorption in a plasma cylinder with a constant 
(over the cross section) plane and temperature was analyzed in (Ref. 33, 
39). The absorption of long wave oscillations in a nonuniform plasma 
cylinder was studied in (Ref. 40). 


In the case of a plasma with a large density and magnetic fields 
with a high strength, the skin depth 5; (44) is small. When ô; is less 


than the plasma radius, plasma heating by an Alfvén wave is ineffective, 
since the wave energy is absorbed only by ions located on the periphery 
of the plasma cylinder. 


In order to avoid this difficulty, Stix (Ref. 41) employed the in- 
genious idea of "magnetic beaches". If the plasma cylinder is placed in 
the field Hg with a slowly decreasing strength, then cyclotron damping is 
exponentially small in the region wy - w>>kjj vy. An Alfvén wave readily 


penetrates the plasma, and they may be excited resonantly. When the wave 


is propagated along Hj, the difference wj - w decreases, while the refrac- 
2 


ae ae — m? 
tive index n|| ~ 1/ w; —- w and the damping coefficient e~ exp [aa | /56 
wre 
increase. If the magnetic field decreases smoothly, the reflection co- 
efficient will be small, and the wave will be absorbed at the approach 


to the region of strong cyclotron damping, where wy - w < kil Vi (region 


of the "magnetic beach"). The field behavior close to the magnetic beach 
was determined in (Ref. 42, 43). Absorption of Alfvén waves at the magnetic 
beaches was determined experimentally in (Ref. 44). Experiments performed 
at Princeton closely coincide with the theoretical computations of cyclotron 
damping in a linear approximation [see the summary in (Ref. 45)]}. On the 
other hand, in several devices for ion cyclotron heating (see, for example, 
[Ref. 46]) the variable field strengths are on the order of, or even con- 
siderably larger than, the critical strengths. Under these conditions, 

we must expect a decrease in the field energy absorption, as compared with 
the case of weak fields. 
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The above-mentioned difficulty entailed in high frequency heating 
of a plasma having a large density by a strongly damped Alfven wave is 
unimportant for a rapid magnetosound wave. A magnetosound wave is ab- 
sorbed in the case of w % wj to a considerably lesser extent, and can 


penetrate the plasma readily. The damping coefficinet of this wave in 
the case of |w - œ| < k] vy is on the order of (Ref. 28) 


L 
Snn T, 
x oti 
žy e (53) 


where nj % = 


If kya ~ 1, we obtain the following expression (Ref. 39) for the 


energy absorbed on the average by one ion per unit of time from the formula 
(13), taking into account (53) 


xs \2 1/3 
dw |H Ho 
dt (i) orat] : OY 


Plasma heating by a magnetosound wave is greatly intensified, if the 
frequency w uw; coincides with the eigen oscillation frequency of a plasma 


cylinder w In this case, the field strength in the plasma increases /57 


res’ 
3 —< 8x97; . 
in the case of lores - w | sw —s as compared with the non- 
H 
0 
resonance case, by a factor of v H .: 
8nngT; 
/ B (55) 
H: ~ H 
BnngT, 


It is apparent that we then have 


~ \2 2 3/2 
dw [H Ho 
We ~ (#,) oT: e 2 | (56) 


Multiple resonance: w = 2w;, may be employed to heat a plasma with a 


large density by a magnetosound wave. The damping coefficient of the mag- 
netosound wave in the case of |w - 2u,; |< ki; vz is (Ref. 28) 


tn (et ): (57) 
Rye 


where asy Te is the ion Langmuir frequency. It follows 
i 
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from (57) that in a plasma having a large density Q; > ky ¢ resonance 
at the multiple frequency w % 2w; is more advantageous than at the main 


frequency. 


The heating of a plasma cylinder which is uniform across its cross 
section was examined for wx nu; (n = 2,3, <...) in (Ref. 47), and the 
case of wavelengths (kj a << 1)in a nonuniform plasma was investigated in 


(Ref. 48). Cyclotron absorption of magnetosound waves in the case of 

w = wy and w = 2w; has not been determined experimentally as yet, although 
the inverse effect -- cyclotron radiation of ions in a dense plasma -- 

was studied recently (Ref. 49). 


The statements presented above illustrate the following: 


(1) Cherenkov and cyclotron plasma heating with weak fields, when 
distortion of the ion distribution function is compensated by collisions, 
occur at the same rate as ohmic heating by a "constant" electric field, /58 
whose strength is less than the critical strength; 


(2) Plasma heating by weak fields up to thermonuclear temperatures 
(Tj ~ 50 kev) takes place over a long period of time. However, this time 


is less than the time for containing a plasma in a thermonuclear reactor 
with a positive balance; 


(3) When heating is performed with fields having subcritical strengths, 
it is more advantageous to employ a plasma with a great density, since col- 
lisions occur more frequently in it, the critical field strengths are larger, 


and the heating time is less than Cnn. 
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DIELECTRIC CONSTANT OF A PLASMA IN A DIRECT PINCH MAGNETIC 
FIELD AND IN A DIRECT HELICAL MAGNETIC FIELD 


V. F. Aleksin, V. I. Yashin 


~~ 
a 
© 


Dielectric Constant of a Plasma in a 
Direct Pinch Magnetic Field 


The study of electromagnetic oscillations in a nonuniform plasma 
has aroused a great deal of interest recently. This is primarily due to 
the discovery of several instabilities. A great number of these studies* 
has been devoted to studying the electromagnetic properties of a slightly- 
nonuniform plasma located either in an almost uniform magnetic field with 
parallel force lines, or in an axially symmetrical magnetic field with 
helical magnetic force lines. 


In conjunction with these cases, when investigating the problem of 
controlled thermonuclear synthesis it is very important to study the plasma 
electromagnetic oscillations in pinch magnetic fields and in helical mag- 
netic fields and in a stellarator with helical current winding (Ref. 3), 
which has a more complex structure than magnetic force lines. 


As is well known, the electromagnetic properties of media are de- 
scribed by the dielectric constant tensor or the electroconductivity tensor 
related to it. In slightly nonuniform media, one can introduce quantities 
which are similar to the electroconductivity tensor and the dielectric 
constant tensor. In contrast to a uniform plasma, these tensors depend on 
spatial variables in wave vector space k and the frequency w. 


If E(k,w) is the electric field strength, then the density of the 
current induced by this field can be written as follows 


ja(r, t) = f dkdwoug (k, w, r) E(k, w) ef =op, (L) 


where Tyg lk, w, r) is the electroconductivity tensor of a nonuniform plasma, 


which is related to the dielectric constant tensor by the well-known rela- /61 
tionship 


Eap (k, w, r) = Sap + H op (k, o, r). (2) 


The explicit form of the tensor oag tk, w, r) may be found by different 


* A detailed list of the literature is presented in the review articles 
of A. B. Mikhaylovskiy (Ref. 1), A. A. Rukhadze, and V. P. Silin (Ref. 2). 
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methods. In our case, it is advantageous to employ the method advanced by 
VD. Shafranov (Ref. 4). In order to do this, it is necessary to determine 
the trajectory of the unperturbed particle motion and to solve the kinetic 
equation by the characteristics method. 


In the absence of a balanced electric field, the unperturbed charged 
particle motion may be described by the equations 


d d 
zaa; G=v, (3) 


where B is the magnetic field strength in the plasma; e and m -- particle 
charge and mass, respectively. In the cylindrical coordinate system, the 
strength components of a nonvortical, direct pinch magnetic field have the 
following form 


B, = Bo D% gn sin azn; 
n=l 


Bz = Bo ( + y fn cos anz); (4) 


nal 


where Bọ is a uniform magnetic field; a = 2; £,() and 8, (T) -- functions 


of the coordinate r which are connected by the relationships 


df d(r 
pe a nagn; ae) 


dr 


= rnafy. (5) 


In the presence of a plasma, the field components, which are related to 
the pressure gradient and the longitudinal current and which may be found 
from the equations of plasma equilibrium, must be added to the nonvortical 
field (4). 


A solution of equations (3) 'in the general case of arbitrary fields /62 
entails considerable difficulties. Let us examine the case (which is 


of practical interest) of a low-pressure plasma [86 = ae << 1) without a 


longitudinal current, in a pinch magnetic field (4) having a large uniform 
component Bg(f,~ g,~ 6 << 1). 


In solving nonlinear equations (3), the presence of a small parameter 
ô makes it possible to employ -- along with the drift approximation -- the 
method of averaging (Ref. 5) when solving drift equations for flying parti- 
cles. The number of particles which are blocked is small, and their contri- 
bution to the tensor og may be disregarded. Avoiding cumbersome computa- 


tions, let us derive the final result of solving the equations of motion (3) 


2 
within an accuracy of terms on the order of = (5 = go2) inclusively 
B 
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r(t) Spay a (cos na (v ié + z) — cos naz] — 


a=1 
vi n . ” 
S [sin (y — 9st) — sin y); 
20, + v? 


2v irog 
n=l 


=p HE+ ga [sin na (v é + z) — 


— sin naz] ta — zz [cos (7 — ost) — cos 7); 


(6) 
AA : x 
2(f) =z +v — z [sin na (v jé + z) — sin naz]; 
2nav i 
nal 
v =v de vfs [cos na (0, t + z) — cos naz]; 
1 
v, (t) =v, + 5 z Vafa [cos na (vf -+ z) — cos naz]. 
n=l 
where v| and vj are the velocity components which are longitudinal and 163 
transverse to the field, respectively; y -- the initial phase of particle 
rotation around the center of a Larmor circle; Vo -- the averaged velocity 
component related to the drift of the Larmor circle center: 
2 2 4 
I, = ihe PL us! 
Upma Top u, + rog u+ ro poi far 
rp, (7) 
u, =— $ D gaga -4 D hho T r 


n=ì 
n=l 


For purposes of simplicity, we have omitted the index 0 for quantities 
taken at the initial moment of time t = 0O in formulas (6) and (7). We 
shall also do this from this point on. 


By solving the kinetic equation, with no allowance for close collisions, 
by the method of characteristics, we can find the expression for the electro- 
conductivity tensor in a cylindrical coordinate system from the equation 
for the induced current density 


Sap (k, œ, r) = =D f dvv. Jawol (s w mae ka + (8) 


ug(t) 
tre — ays 


[v (2) e2],8p, } e~it + kL rdt, 
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where Fp is the equilibrium distribution function of the particles; k -- 


the wave vector with the components k}, k,, kg (ky = a m ~- whole numbers); 


z -- the sum with respect to ions and electrons. The electric field may 
e 
be determined by the expression 
E (r, o) = X f dk,dk:E (kr, ka, m) i*r tH z+, (9) 
m 


For a slightly nonuniform plasma in a strong magnetic field, we may 
write 


P= (1 + živi VY zg) Fw Ui, V), (10) /64 


for the equilibrium distribution function within an accuracy of terms on 


Vv 
the first order of smallness with respect to the parameter aa (vr -- 


thermal particle velocity, a -- characteristic dimension of the nonuni- 
formity for the main plasma state). Here, t -- unit vector in the field 
direction; f -- arbitrary function of the velocities v], vj and of the 


variable ¥ which is an integral of the drift equations. Within an accuracy 


of the terms =, the quantity Y coincides with the integral of the force 


line equations 


Y =r 4 2r 2o cos naz. (11) 


n=l 


Due to the smallness of the parameter ô, we may approximately compute 
the function f which depends only on the coordinate r. We should note that,in 
order to avoid this, in a plasma with the selected distribution function 


2 
div j #0 we may add the component vB x (Eig vfa to the function Fo, 


where integration is performed along the magnetic force line. However, as 
may be readily seen, in our approximation this term is small. 


Let us select the Maxwell distribution of particles with a nonuniform 
equilibrium density no and the temperature T, which depend on ¥, as f: 
_ m \" m tor (12) 
F= no (5%) exp(— g - 5 *4), : 
Retaining the important terms in (8), after simple transformations we ob- 
tain 
: o{ AA 1 2 
Sap (k, w, r) =in Ve fA AE vi, V) Qigdv,dv;, 
i e : (13) 
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where the following notation is employed: 


A 2 es 
M =i ÈJ: (n) Jela) Ba s (n zten =) > 


m = {Mela [iie ae (a to) PI 
an 2nav* ws à J? 
2 2 2 
na ky (2; +oi)ea _ Fifa 
20,2 2nav? 


$n = arctg 


k, En 


Ý kivY) T 3 əm ð 
ra (marta T 


(14) 


0 


Qn 
Qag = -áf dv. f Up (£) A (£) dt; 
0 


A (t) = exp {— i (w — nasu, —kyvg — hv.) t — È X 
x [sin (y — wst — ) — sin (y — ẹ)}l}; 


k k T ARTE 
} = arctgz; p= a, hi =V RFR. 


B 


The velocity components which are included in the expression for the ten- 


sor Qag in a cylindrical coordinate system have the following form: 
v, (t) = 1 cos (y — ost) + D ads sin na (v,f + 2); 
n= 


d (t , 1202 +02 | — 
r0 =v, sin (y — osf) + D, E fn cos na (oat + z) + O53 (15) 


n=l 


a 
v: (t) =v, — Dia fa cos na (vit +2). 


n=l 


By employing them, we can obtain the explicit form of the tensor Qag 
Qag = Dar (16) 
where the vectors q and q have the following form 
TEA [2u, (© cosp + iy (8) sin o] Ay Ep ©; 
g= 0 |E Jp È) sin p — iJ, (È) cos o] +T O; 


v? B,—B 
g= g — at e)n O; 
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g- = 01 [Ze © cosy — iy © sing) P —F 0, J, O Deeper, oe Lie 
Je = 01 [2 Jp Osing + pcos] ch + Dodot) M; 


a v? 
dz = Vado (ep? — o To © Shien’ e; 
c® = (o — pog — nasu, — lav, — ke0 — kav p) 
(Jp (8) is the Bessel function; the sum with respect to 1 is taken from 
-æ to +”, and it is thus assumed that Bg =T 8? fL = j and fo = 0). 
As is well known, the electroconductivity tensor may be employed 


to obtain the expression for the polarizability vector x which charac- 
terizes the density of the induced charge p, 


p = | dkdue! — YF, (k, w). (18) 
The components of the polarizability vector have the following form 
AA s ~ 
Xa siz Ne f dord ML dF (0s, v, V D Oge (19) 
i,e p=— o 


The expressions for the electroconductivity tensor and the polariza- 
bility vector may be simplified in the case of electromagnetic oscilla- 
tions, whose wavelength is much greater than the pinch modulation depth 
(na << 1). In this case, we must get rid of complex sums and products, 


since the terms with s = s' = 0 and the operator M = 1 will make the 
main contribution. 


Plasma Dielectric Constant in a Dielectric Helical 


Magnetic Field of a Stellarator 


The strength components of a nonvortical magnetic field from a helical 
current winding with the finite step L have the following form (Ref. 6) 


B, = B, D En sinnd; Be= Bo 2 fa cos nf, (20) /67 
n= 
Bz = B, — arBe, 


where Bp is the strength of a uniform, longitudinal magnetic field; 


0 = > - az; a - 2, Sy (x) and f,(r) -- functions of the coordinate r, 


which are related by the following relationships 
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(8a) = nfa(1 + a?r?); 
(rfa = ngn. 
Just as in the case of a pinch. field, when deriving the explicit 
form of the electroconductivity tensor, we investigated a plasma with a 


small gasokinetic pressure (B << 1), without a longitudinal current, in 
a helical magnetic field having a large axial component 


(21) 


(fa ~Ba~ 8K 1). 
i 2 
Within an accuracy of terms on the order of £, the particle tra- 


jectories and velocities may be described by the following expressions 


r(f)=r+ Yi%[cosn (0 + mt) — cos no] — giint — 
n=l 


— ost) — sin y]; 


sO=—e 24 F fa [imni —sinn( + ea 
n=l 


x [cos (y — gt) — cos y]; 


z(t) =z+ vt + z xt [sin n(o + a) — sin no; (22) 


u, (t) = v4 cos (y — vst) +o y gasinn (0 + w), 


a! 


7D Son sin (y — ost) + v; Yprcosa(ag 2 2A + os 


n=ì 


vz (t) = vz = ar Mia cos nfo Ta 
n=} 


where Vos Vz and Vg = Vo - arv, are the averaged velocity components re- 


lated to the motion of the center of a Larmor circle along the averaged 
force line and drift in a nonuniform field: 


ty = orto, + Sey + hug + 3s 
i 


rop 


o talg, (23) 


raog 


2 =v, —4 o1 DR + e) — 


n 
Here x is the torsion angle of the magnetic force lines: 


o 


= son xf + g — 2 nfn + a®rf2) (24) 
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œ 


i = g X (fè + 8 — 2nfagn — arnaga); 


4nrp’ 
uz = ty (f2 + gt — 2nfnga — a?r?nfaga) + Be? (25) 
a 
1 


The distribution function Fp in a slightly nonuniform plasma may also 


be determined by (10), in which ¥ describes the equation of magnetic sur- 
faces: 


g, 
Y = r? —2r Y= cos nð. (26) 


a 


We again select the function f as a Maxwell function, with the density and 
temperature dependent on the magnetic surfaces Y. 


As a result of the computations, we obtained the following from the 
general expression for Fug (8) 


5 AAI a z 
Cup aia Je (MEET. ow Y i (27) 
ie p=% 


where 


A 2 ss’) nd baa 
= TLD Ide tage etal, 


n=l s,s’ 


i 2 2) 1/3 
1 2 2 262 k, v 
„=i [Metter (IgE Sar) | ; 


2 
Pa = arctg cet — or ar); 
? 
` p 7’ e B, 
je Uy [2 Jp (È) cos + i; @sing| + g 21T E) 
à P i ay oa By - 
g= 01 [2 Jo @)sing— il, O cos] + Delo O + oO 


SR v B 
q= (Te — ar z+ pt) V0: 
R (28) /69 

Ca [2 J (8) cosy — iJ (8) sing | —49 44, (0 x 


1 
x y gic el”; 
mig 
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Ge = 01 (F Jo Osiny + if, E) cosp) + Gola © oS + 


+ 72st) Shovel 


ton Aden (0) aroi, J, (t) X fic et; 
qz = U2 p (E) cp gji 40, p È) 7 ficp , 


= = nsvg lo m 
The quantities Kgs kj» Y, & and the operator L are determined by the ex- 
pressions (14); it is assumed that S- = -8> -£5 f and fọ = 0 in the 


sums over the index 2. Correspondingly, the plasma polarizability vector 
in a helical magnetic field is 


pe AA ` ~ 
Xe ain Yet [MLZ flor, v Y) X J È) qudv,do’. (29) 


p=—% 


Just as in a pinch field, expressions (27) and (29) may be simpli- /70 


field, if the oscillation wave length is much greater than the difference 
between the maximum and minimum radii of the magnetic surface (na << 1). 


We may employ the expressions obtained for Sug and X, to study the electro- 


magnetic oscillations and plasma stability in pinch magnetic fields and 
in helical magnetic fields. This will be the subject of future research. 
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SECTION II 


LINEAR PLASMA OSCILLATIONS /71 


KINETIC THEORY OF ELECTROMAGNETIC WAVES IN A CONFINED PLASMA 
A. N. Kondratenko 


Problems of electromagnetic wave propagation in a confined plasma 
are of considerable interest to studies on methods of plasma heating, 
acceleration of charged particles, plasma diagnostics, and other possi- 
ble applications. The hydrodynamic theory of plasma wave guides for slow 
waves has been studied in great detail [see, for example, the articles 
(Ref. 1 - 4)]. However, the hydrodynamic theory does not encompass the 
important phenomena related to the particle thermal motion -- for example, 
wave damping which is particularly great at smali phase velocities. 

Since the phase velocity of a propagated wave Vg is less in the wave 


guides of slow waves than the speed of light, and since it may be compara- 
ble to the mean thermal velocity of electrons vp, or ions vp;, the neces- 


sity of a kinetic examination becomes readily apparent. 


On the other hand, the confinement of a plasma leads to a new type /72 
of wave -- surface waves -- whose damping, as was shown in (Ref. 5, 6), 
is proportional to the thermal velocity of plasma electrons for vag << V9. 
In a nonconfined plasma, where there is no surface wave, the damping of 
the longitudinal three-dimensional wave is exponentially small (Ref. 7). 


Formulation of the Problem 


Let us investigate the propagation of slow electromagnetic waves in 
a plasma layer which is 2a thick in one direction, and is not confined in 
the two other directions. As is known (Ref. 4), waves propagated under 
these conditions are surface waves when there is no magnetic field. 


A self-consistent system of equations describing these processes con- 
sists of the Maxwell equations 


1 OH, 
rotE=—=2- a; (1) 
_ 1 OE , 4x, 2 
ROU Se op rr] 4) 


and a linearized kinetic equation for the deviation of f, from the equili- 


brium distribution function fo, of the a type of particles (a = i -- ions; 
a = e ~- electrons), in which we shall disregard particle pair collisions 
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D nen f 2 


E+ Liv, HI) =0. (3) 


The standard notation is employed in equations (1) - (3). We shall select 
the coordinate axes so that the z axis coincides with the direction of 

the wave propagation, and the x axis is perpendicular to the layer. The 
YZ plane is the plane of symmetry. 


Equations (1) ~ (3) must be supplemented by the boundary conditions. 
Let us assume that fy = ft + fy, where Et is the distribution function 


for v, > 0, and fo -- for vy < 0. We shall select the conditions of the 
mirror image (Ref. 8) (the final result does not depend quantitatively 
on the reflection condition) as the boundary conditions for the function fò 


fi (44, ve >0, v2) =fa (4a, ox <0, v). (4) 


We obtain the boundary conditions for the fields from the Maxwell equa- 173 
tions (1) and (2), integrating them over an infinitely thin layer which 
encompasses the plasma-vacuum boundary: 


FP1( +a) = Eyac¢ta), HBH +a) = HY% t a). (5) 


Dispersion Equation 


We can write the dependence of the distribution function fq on 
time and the cordinate z for the fieldsin the form exp i(k3z - wt). If 


fo, is the function of energy, we obtain the following from equation (3) 
with allowance for the ee conditions 


f(x) = ae] Ba) get tE Oa Te 
à F iy (a-++x) a ðfoa (6) 
Mon o cos y (a —}) — 


ee 


— iE, oz * sin y (a — 9l. 
where 


= kv,— o 
Uy 


In order to compute the currents j, and j,, it is necessary to deter- 
mine the sum and the difference ft + f,- Employing the values (6), we 
obtain 
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2i 0 
RHR =z. Te Fae Kile t) dt — 


—a 


ape ged. z | BO Kal E) dt; 7) 
fi-pe—pe geet. 8 eke tee 
ret ( nex oe n 
where tmas a XLE, 
, ane cos y (x — a) cos y (E +a), x > & (9) 
HDE ames ae os (10) 


The kernels Kı and Ky, which are expanded in Fourier series, have 
the following form 


A ‘COS a ae x 
4 
Kejsi 42 (11) 
n=0 
sina af in tn x 
7 
K: (x, p-s va , (12) 
where a, = ai the prime over the sum indicates that the sum term corres- 


ponding to n = O must be multiplied by 


Employing the values (7), (8), (11) and (12) and using equation (1) 
to change from the fields Ex, E, to the fields E,, Hy» we may obtain the 


formulas for determining the currents: 


G sina,* at (= (p+ 
oe an $E: cos oat wane oe Oe 
P ydos Af, , 
+ iĝ fa sin ont Je E — až i J ? (13) 
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xN 
N 
E> 


, XY cosa, x 2ie? f v,dv, k dfoa 
ie= y [inom a, i fm (ey h- 
0 


a=0 


d 
Z a zt) — Ban $ Hsin nidi y Ji = (aca a A 


aai" 


where B = —. 


The solution of the Maxwell equations (1) and (2) for determining 
the fields E, and Hy leads to the following equations 


dE, ok Ank 

Ee iB P) Hy = Sis (15) 
dH, o ; 
Ge Hip Be = Sie are) 


We may write the solution of the integro-differential system of 
equations (13) - (16) in the following form 


E- (x) = SE zn COS Ox; H,(x) = alt yn SİN Arx. (17) 


We obtained the following value for the Fourier components E p? 


Beas (—1" 3 Ha) Ab f thot —F) p (18) 
where 
a= (i2 A) [A +80 p| + (Aa — an) (Aa + aa; (19) 
P 14%, foa. 
A, att i Ka EA E a? ðo, ’ 
—co 0 n 
2 May ; ( do, fofoa q foc). 
D e fa e we (E e): (20) 
2 h T do, foa 
iaaa “Sat Soto.) ty, T 
inol , C d fon a? oa 
A, =2- Do f odo | os hra = PE . a 
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4ne2p 
Here g2 = a0, po —-- equilibrium ion density which equals the 


m 
a 

equilibrium electron density; fia is the density normalized to unity 
OE a dv = 1). 


Selecting Maxwell distributions as the equilibrium distribution /76 
functions fog? we may write the equation for the plasma impedance 


(21) 


a=é 


E (a) 2 byw’ : Ve 
iH @= eA aTe 1— 25 Qia . 
We thus have 


e 


Q? a? ST 
A= — k 1 Ee (ent fon -e(1- #a.)]- 
3 


a=i a=i 


e 2 A 2 
-a(1=r 5e) 1— YB (Qie+ Qa) |; 
a=i amt . (22) 


Qia = 2 on f dz (oa — bz) e7? [loa — 62)? Ia (2) — Val; 


Q= Zo § dzz? (oa — bz) 7l (2); 


(23) 
Qaa = 2 aab § dzz (0a — bz)? e7 Ta (2) 
[s] a ; v = oF, 
Č apta? m, —- the mean thermal velocity of a type of particles; 
b="; h(g) = (ee 
a, (Sq — 52)? — x? J ° 
Equating the plasma impedance (21) to the vacuum impedance 
E, (0) c £ 
AOT Ekl — BP, (24) 


we obtain the dispersion equation 
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; © o2 a 
Zk x] i —6 i idi = (1—§?)?. (25) 


n=0 ast 


The sum (25) consists of the real and imaginary parts, since each Q con- /77 
tains the real part which equals the main value, and the imaginary part 
which is proportional to the residue of the integral I, (2). 


High Frequency Oscillations 


kare 
w 


Let us set 0n =Q, <1, m>. 


We shall employ the following notation 


bases g? 

i= ere X p[i-e(i že): (6) 
n=0 
a. T 

Pa = Im) Lji—e( - 3a) (27) 
n=0 


Since a, are evenly included in the sum (26), we may change from the 
sum to the integral (Ref. 8) 


i c a o? 
Pi =i ks | a hı —e(ı a Qı (@)| dq, (28) 


c 


in which integration over the contour C is performed from -~ to +œ, passing 
around the singular points q = a from above. Let us deform the integration 


contour, and let us enclose it in the upper half-plane. Let us set 
qo = ik, -- the value of q at which A(q) = 0. We then have 


1 — Be 
Py =k, . —; cth ka, (29) 


since the integrals Qj have the following principal values 0} (j = 1, 2, 3) 
when the electron pressure is disregarded: 


, 3 s g2 
Qi =Q:=1, Q3=0 (=1-5 =a J. 


When computing the sum P2, we should note that -- since da5 


>> 1, the integral residue 


=b >> b -- for small n, just as for o 


k3’Te 


I(z) is exponentially small. Therefore, components with n which is 
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larger than a certain ng (i.e., shortwave components of Fourier expan- L78 | 
sion), which may be determined from the condition o, ~% 1, make the main 
contribution to the sum P2. For such n, b2 << 1, and Qj have the 

following values (we shall omit the index e): 


Q: = 2| —1 + ozh (oferat — ive)|: 

Qz = oe” (25 ena — iVa); (30) 
“0 

Qs = 206 |— + (ot — je (2ferae a ivz)| 


Since Q} % 02, Q v0, 23 ~% o, in the case of b? << o2, b? << 1 for A, 


which is determined by equation (22), we obtain the following expression 


2 o? 
s= àli =% o). (31) 
Let us employ the following notation: CH = ReQi, Qi = Im Qı. We 
then have 
2 o? oF oQ; 
Pa = — 7 hs ae \ | o2 .\? /@ ,\2° 
(1-7) + (Ze) (32) 
n> Ny 


Since ng >> 1, we may change from the sum to the integral 


P: = ak Se, BP re ’ ooo he 
T Va w? w g2 ,\2 Qg? ,\2° 
7 bea <i) +(% ci) Ga 


We may disregard the quantity a in the denominator of the integral, 
and for Qi we may employ its value in the case of o V1. In addition, 
making an exponentially small error, we may multiply the lower integration 
limit by ©. We finally obtain 


2 plat Ée, zq SOY (34) 


UB cthea +i- EE. Saw. (35) 


Let us assume that k3 > k3 = k3 + iô, 6 << k3. Then the real part /79 
of equation (35) produces a relationship between the phase velocity ß and 
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the frequency of the propagated wave w: 


1 
— B2)? al 
CEY cth kia = (1— B}, (36) 
and the imaginary part produces damping: 
2 Bee QN Ë 
ye o pala, p ,thke a" (37) 
“Tae "19 aire 


Formula (37) can be considerably simplified in the case of large 
(ky a >> 1) and small (k| a << 1) layer thickness. We have excluded 


B2 from equation (36), and shall substitute it in (37). In the case of 
kj a >> 1, we obtain 


zr o% dT ep? è?’ (38) 


where k = = The damping determined by this formula coincides with the 


damping found in (Ref. 5, 6) for a reflection coefficient of P = 1. In 
the case of k] a << 1, we have 


1 
p2 a het =! at CL eae. (39) 


It may thus be seen that wave damping is increased when there is a 
decrease in the plasma layer thickness. 


Ion-Sound Waves (Vīe >> Vg >> Vri) 


The dispersion equation (25) may be considerably simplified, if we 
set the speed of light c = œ. In this case, we have 


a 
b= 28 sifi -DA to +20. (40) 
a=i 


A A 
We shall employ the following notation: A" = Re -k2 A' = Im k2" 
3 


Then the real P3 and the imaginary P, components in the sum (25) have 
the following form 


2 Boe SAP 2 
Pee 2s at pan’ (41) 
n=0 
2 A’ 
P, ane akz A” + A” (42) 
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We shall employ the same procedure in computing the sum P3 as for 
computing P], but we shall take the fact into account that in the case 


of oe << 1 we have 
Ref (4) + 3 Q: (9) + 203(0| = 2a (43) 
We obtain 
ene. (44) 
where Cab: be eae -i n=. 


Employing the integrals (23), we find 


2V x03 2 
PERS areas a ee ee (- rep) (45) 
(1 + 6)? 
o2 
e 
Since Ta Be << 1 for any n, we then have 
2 
,__ 2Vxb g? o Nee TO er NS 
am ae (ea) © FP Fi (nea) | Gi 
(14-67)? 


Consequently, both electrons and ions make a contribution to the wave 
damping. The ion component Pua of the sum P, has the same form as P, with 


a replacement of the indices e > i: 


| ne ae ea LD (47) 


In computing the electron component P,e of the sum Pis and in computing 


Ppi» we shall set 4A"? >> A'?,. We then have 


( ea)? 

* {itso} 2 

y l—e 2 

TE =i A Sal ie J ‘st 2° (48) 


Thus, the dispersion equation (25) has the following form for ion- 


sound waves 
k3 cthza 


e+ (Pu + Pa) = 1. (49) 


We may thus find the expression for the spatial damping decrement 


72 


/81 


—2,(P 
3 = k3(1 + &) ae L ute . (50) 


fe, +s 


sh?xa 


In a nonconfined plasma, where there is no surface wave, the damping of a 
three-dimensional ion-sound wave is determined by Cherenkov absorption of the 
wave energy by plasma electrons, and the ion contribution to the damping 
is exponentially small. As may be seen from formulas (47) - (50), plasma 
ions also make a great contribution to the damping for a surface ion-sound 
wave. The physics here is the same as for high frequency oscillations -- 
the Cherenkov absorption of the wave energy by plasma ions is particularly 
significant for the short wave components of Fourier expansion of the 
surface wave. 


Let us study equations (49) and (50) in special, different cases. If 
k3a(l + &) /2 >> 1, then we obtain from equation (49) 


2 
ical ie (51) 
Eine 
4 G 1— e [* ne 2\2 | (52) 
6, ee 1 — £7)“ | eif 
yee a e a a A 
k3a 


since for large k3a the sum included in Pre approximately ar a 
T 


It can be seen from formula (52) that the ions make a particularly {82 
large contribution to the damping, if e4 > -1. If k3a(1 + &) le << 1, 


then ane 
poeta lise ||, 
ST te} |’ (53) 


d= katt È $ (Par + Pae). (54) 


The following limiting cases are possible: (1) At a ee plasma 
= + <1 <1), 


electron temperature or a small layer thickness, when — Bet 


the damping decrement is 


(55) 


var + ate 
5 E; OTe 


1 
(in the case of k3a << 1 the sum included in P,e approximately equals D: 


(2) For a large plasma ee thickness or an electron temperature which is 


>l!) 


not too high when PE a 


Te 
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7 
y 1— 2 awole 8 
b= 1/ 2+. dem a| (56) 


R je Ë Ü | ao e? V 898 


We cannot assume that a > œ here, due to the condition of the initial ex- 
pansion. 

As may be seen from formulas (52), (55) and (56), the damping decre- 
ments of a surface ion-sound wave are large in different cases. Plasma 
ions make a significant contribution, and sometimes the main contribution, 
to the damping. The wave phase velocity is decreased with a decrease in 
the layer thickness; therefore, the ion contribution to the damping in- 
creases, while the electron contribution decreases. 
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KINETIC THEORY OF A SURFACE WAVE IN A PLASMA WAVE GUIDE 
M. F. Gorbatenko, V. I. Kurilko 


As was shown in (Ref. 1), when a surface wave is propagated along a 
plane boundary of half-space occupied by a plasma having a temperature 
which is different from zero, Cherenkov absorption of the wave energy by 
thermal plasma electrons takes place. This leads to damping of this wave 
even when there are no collisions. In contrast to the damping of a longi- 
tudinal wave in an unconfined plasma, in this case for small thermal 
velocities the damping coefficient is proportional to the thermal velocity 
of plasma electrons. 


We investigated this phenomenon for the case of a plasma wave guide 
produced by a plasma layer having a finite thickness (21, |x| < 1). 
Since the electrons moving at a thermal velocity are successively reflected 
from both walls of the wave guide, it was not known previously that in 
this case the absorption investigated in (Ref. 1) will not decrease con- 
siderably. 


The initial system of equations consists of the kinetic equation for 
a high frequency addition to the distribution function and of a Maxwell 
equation: 


OF peas > a 
at. oyf — = Ey; o =0 


rot H = S744 a 
gs (1) 
+ 1 OH 
rot E =— => ar 
— 4-0 > 
J=—e f vfdv 
It is assumed that the equilibrium distribution function is a Max- [84 


z A n vpo? 
well distribution: f, = —* exp = 
T U. 


5 | (the z axis passes along the 
Ur 


T 


direction in which a surface wave is propagated). 
Let us write the solution for system (1) in the following form 
f(x, z, t) = 9 (x) exp [i (kaz — of )], E(x, z, t) = 
= E (x) exp [i (kz — of]. 
We then obtain the following equation for the distribution functions of 


particles moving in positive (+) and negative (~) directions along the 
x axis 


75 


as 


—iwtọt +, - Ele z (vE: +v,E,) = (2) 


We shall employ the following geet: as the boundary conditions for 


the distribution function $* (x) in the planes x = +1 of a plasma layer 
gt (x = — l) =p (x = — h), 
g(x =l) = pet (x = 2. 


The solution for equation (2) which satisfies the boundary conditions 
has the following form 


(0 <p< 1) (3) 


mv,v 2A 


ot(x) = —- 2i l f (v-Ez + 0,E,) exp ie w x 29 |e + dx! 4. 
L 
te | et ntaen [Ew aaja feg tae 
1 


x exp a a (x —x— 20) aa (4) 


=x) = Sar (f (v-Ez2 — vsEx) exp [o (x’ —x— das 


t x 
—p f (v-Ez + v,E,) exp k eta] ae — p? | (v.E.— O,Ex) X 
“1 * 4 
x exp [a Beeria, 


where 


Die LO 
0 m » 
No -- plasma electron density. 


By substituting system (4) in the expression for the current, and 
by substituting the latter in the Maxwell equation, we obtain a system of 
two integro-differential equations for the E, and E, fields. In the 


general case, this system is cumbersome, and it is difficult to obtain an 
analytical solution for it. However, in the case of slight thermal 
scatter (vr >+ 0), a solution may be found for the system. In actuality, 
in this case, as may be readily shown, current corrections due to the 
electron thermal motion decrease with an increase in the distance from 
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the surface of the plasma layer proportionally to 


x 


api Re|t—at(R=z, n=); 


Since there are significant eee scatter phenomena only in a layer 


having a thickness on the order of YI close to the plasma surface, they 
w 

may be taken into account by means of effective conditions which the 

fields determined according to the Maxwell equation for "a cold" plasma 

must satisfy, 


The effective boundary conditions follow from the Maxwell equations, 

in which the currents are expressed by the fields by means of system (4). 

For purposes of definition, let us investigate a wave in which the longi- 
tudinal field component is symmetrical with respect to the plane 

=0 [E > (-x) = E (x), Hy (-x) = -H , (x) J (Ref. 2). As would be expected, 


the currents Jy. and J, -- according to system (4) -- retain the symmetry 


of the corresponding fields. Therefore, it is sufficient to obtain the 
boundary conditions for one of the plasma layer surfaces; they will be 
satisfied on the second surface due to the field symmetry. Let us find 
the boundary conditions which must be satisfied by the tangential field 
components on the x = 1 surface. Integrating equations containing J, and 


Jz, over an infinitely thin layer et << € << 6, l, where £ -- layer /86 
thickenss, ô -- depth of skin-layer), we obtain 
i 4al 
-H =% =E Ja (n) dr; 
(5) 


E! — Et = ily aa KKA (n) dn. 


When the integrals in equations (5) are computed, it is advantageous first 
to perform integration over n: 


1 
cal J; = ifi a mal ge — a f Ex (n') X 
0 


w 
g l—e 


sd a oo 
vee ( Edt QR 
—2 È Es) dn’ \ — ip? zt ge z|- E, (x’) x 
r ) J f 
I+e (Ed 
xema esha) (ae +2 { Ex (0) da’ | —ip (5x 
0 
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1 
x exp (— 8) { È Es (a) exp i E+ m| dr! — 
—l 
1 n 
— f Es er)exp [z (+ ané ; 
—1 


ara s r 
< =if Be exp r |- a sc ae =x 
1 


(6) 


ò 
D a= ip? k exp e] feex 


1 
x ole (1l Sean bin (Bex (=ef f E: (q) x 
=1 
x exp = e+ mjer - { Ez(n') exp E a+ m) ar), 
—l 


where only terms making the largest contribution in the case of Vr > 0 


are retained. Performing integration over n' in equations (6), we find 


aik 2iR 
one ; EERU) e = +pet 
AY J, = fep (i as | a i 
20 0 e & — pet 
tal 
2(E— E} æ e f+eé 
q Ee) af Bie, (1) exp (— £”) &8dé — aR iR 
3 ò e, E — pe È 
2R zR (7) 
= g 2e È 
4n y To Ei)? 3 e PE 
tS ae exp (— 2) tdt LER IR 
we I R J Xp ( a s — pte € 


[= iR iR}]2 


-2 fE: (1) exp (— £?) Edé 


e 


wN 
œ 
N 


Since Je and J in equations (7) are proportional to vp, when the integrals 
v. 


are computed with respect to € = = in these equations, it is sufficient 
T 


to confine oneself to the zero-th approximation with respect to Ve (i.e., to 
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strive to zero in the vp intervals). Taking the fact into account that 
Imw > 0, we finally obtain 


__ 9 (1 -- PEL); 


I= Ba VaoR lesl >I ies 
2 e i 2 

. _ eo (EZ— Ez) | oU +P) 

t= “arom tyro =) 


Thus, the effective boundary conditions taking into account the 
thermal motion of plasma electrons have the following form 


klo (1 — p) ; 
GN = see: er 
la (1 + p) io? (1 +p) (9) 


ʻi 
a(i =E) aa t pag he mO 
Substituting the solution of the Maxwell equations for a plasma 
Cal < 1) in a vacuum ({n| > 1) in these boundary conditions 


Et (n) = Achzln; Ez (a) =A shrin; Hin) = — ‘Bs A shtin, 
where 
2 
T? = k? — e3k?; eg = 1 — 5; k=; 
Eel) =B —il +e a Beat aly ot ik =i, 
s(n) e » Ez (n) = — «Be , Hi n) E Be , (10) 
x? =k? — k?, 


we obtain the dispersion equation 


iw 


a. ~ i 
Vint [a Fe ea (h pl (11) 


E3 [f -} ea” th | = 


In the case of 1 + ~, this equation changes into an equation for half-space 


(Ref. 1). 


Thus, in the case of a plasma layer having a finite thickness, the 


surface wave damping is proportional to vp. However, the wave damping in- 


creases with a decrease in the layer thickness, because the wave phase 
velocity decreases: 
iw 20, 


Ak, = — —— 
3 V Tkzot3lo3 


L= 
[o + p) Ri, P -+ za]; les 1>> l, €3 < 0; (12) 


I 
Hy = Bt ari AP KI, 


79 


REFERENCES 


1. Gorbatenko, M. F., Kurilko, V. I. Zhurnal Teoreticheskoy Fiziki, 
34, 6, 1964. 


2. Vaynshteyn, L. A. Electromagnetic Waves. "Soviet Radio" (Elektro- 
magnitnyye volny. '"Sovetskoye radio"). Moscow, 1955. 


x 
(oe) 
Ke) 


SINGULARITIES OF AN ELECTROMAGNETIC FIELD IN A 
NONUNIFORM, MAGNETOACTIVE PLASMA 


V. V. Dolgopolov 


As is well known, in an isotropic, nonuniform transmittant medium there 
is a sharp increase in the electric field strength of an electromagnetic 
wave at the point where the dielectric constant vanishes. This increase is 
caused by plasma resonance. This phenomenon, which has been called "infla- 
tion" of the field has been studied in several articles. It is pointed 
out in (Ref. 1) that the "inflation" must occur in a magnetoactive medium 
at the point where the refractive index for the wave becomes infinite. The 
phenomenon of field "inflation" for normal wave incidence on the layer was 
studied in (Ref. 2) for the case when a constant magnetic field was parallel 
to the plasma boundary. 


This article investigates the behavior of the electric field of an 
"inflated" wave close to the "inflation" point in a magnetoactive plasma 
for the case of an arbitrary angle between the direction of the constant 
magnetic field and the normal to the plasma layer surface, and for the case 
of oblique wave incidence on the layer. A solution is found over a wider 
region for certain special cases for a wavelength which is a little less 
than the distance at which the layer parameters change significantly. 


We shall assume that the parameters of the plasma magnetoactive layer 
change along the x axis, and the vector of the constant magnetic field lies 
in the xOz plane, forming the angle = with the z axis. The tensor of the 
dielectric constant ¢;, can then be written in the following form 


e, cos? O + e,sin?Q@ iegcosQO (e — £) cos O sin © 
— iscosð8 , £ ie, sin O 
(e3 — £ı)cos O sin O — iez sin © e; sin?@ + e; cos? O 


(1) 


Eik = 
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The quantities £], €2, €3 will be real, if dissipative processes are not 
taken into account. 


Assuming that the dependence of the electromagnetic field strength 
on time, y and z, is determined by the factor ei(-wt + kyy + k,z), we 


obtain the following system of equations from the Maxwell equation for 
the electric field components 
w? . ðE 0E; w? 
(r? = = eu] Ey+ i(ty rd + kazi) — z (e12Ey + e132) = 0 I 
., GE 2 OE, w? 
iky 30 — F 2E — g + (x2 = Sen) Ey — 


xX 
— (koki + Fees)Ee=01, (2) 


., 3E o? PE. 2 
ike Se — Fe Ee— g — (hy ot Sealer 


w? 
i (r -3 e9) Ee =0 | 


where k2 = k? + k2, 
y Zz 


/90 


The following inequality is fulfilled in the vicinity of the "inflation" 


point for an “inflated" field 


Æ 
Ox 


>a |ë. 
This enables us to represent the field E in this region in the form of a 
series 


-Rop B04..., 


where 
[£o|~ee\[é| lE. 
We shall always select the origin along the x axis, so that x = 0 at the 


inflation point. We obtain the following from the two lower equations 
of system (2) 


ðE!) DELP : 
E% = EO =0, a = ihkyE®) + Cy, ios ik,E) + Cz; 
PED 4 GED a EO PED OEM. a o 
ge S hyg T E aE o Ge Me ge Saas 


where cy and C, are the integration constants. 
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Substituting the expressions obtained in the first equation (which 
is differentiated with respect to x) of system (2) (we shall disregard 


components containing gE 12 and Ila ) and taking into account (1), we [91 
ox ox 
can write the equation for E(°); 
x 


2 (en EW) T i2k.e;;EO + const = 0. (3) 


It can be seen from this equation that "inflation" can only occur where 
€q] vanishes. In the given approximation, it must be assumed that £13 


is constant, and is limited to the linear component in Taylor expansion 
for £]. Then the equation for EC) assumes the following form 


aE k £13 EO const 0 
: (4) 


1+ i124 == 
ox +( 5 a) * ae 


where + ĝen 


We can write the solution of equation (4) in the form of a sum of 
particular solution and the general solution of the corresponding homo- 
geneous equation. The particular solution is a constant quantity, and is of 
no interest. The general solution of the homogeneous equation for the 
field E, close to the "inflation" point yields the following expression 


E, =A ee nen) (5) 


where A is a quantity which is constant with respect to x; 


kar 
= 24. 


Eir 
If o #0, we have 
E = — k A ,—icin (kx) C.: A —is ln (kx) 
yY Y o + u Ez = — kT e + Cs, 
where Cı and Cy are integration constants. 


In the case of o = 0, we have 


E, = ikyA ln (kx); Ez = ik,A ln (kx). 
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Thus, if k, = 0 or €,, = 0, the nature of the field "inflation" È 
is the same as in an isotropic medium. 


If kz #0 and €,3 +0, the closer the "inflation" point, the more 
intensely does the field È oscillate on both sides of it. Thus, the 
amplitude of the field component Ba increases as 4, and the amplitudes [92 
of the components Ey and E, remain finite. It can be seen from equation 
(3) that if den 0, at the point where ej] = 0 "inflation" does not 


ox 


occur in the case of k2£13 #0. 


Let us investigate the case when the medium parameters change slowly 
and the method of geometric optics may be employed when solving the system 
of equations (2) outside of the vicinity of the "inflation" point 


ifa wdx, (6) 


Substituting expression (6) in the system (2), we obtain the following 
equation of the fourth power for x(x) 


ez1xt + 2kaegax? + a (x) x? + B(x) t 7 (x) = 0. (7) 


We shall not give the expressions for the coefficients a, 8, y, due to 
their cumbersome nature. 


If k, and €,3 are different from zero, we have the following from 


equation (7) for the wave vector component x(x) of an "inflated" wave 
close to the "inflation" point 


€ 
LS — 2k: Ž. 
n 


This result coincides with formula (5). The wave oscillates in space on 
both sides of the "inflation" point. 


1. If the magnetic field is parallel to the plasma layer surface 
(€ = 0), €13 and B(x) vanish and equation (7) becomes a biquadratic equa- 
tion. We than have 


= 2_ utip 07 °? 2 — (Gt) 
x? =— k? — e k? 2a + 


w? w w? 
$ a sy [e — ca) (k —ž x gl +4 a Rests. 


(8) 
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Taking into account the components of subsequent order in the expansion 
with respect to the small parameter 3~/x? in the ares of equations (2), 


we obtain the expression for the vector components A(x): 


KEA 
xz g i7 dx 
g Sy SELI ae’ (9) 
Ame ge cope" 1 A =C e , 
where C is the integration constant corresponding to one of the four 193 


functions x(x); 
Y = ky (EHEH 47t) + ka (CF +O) — 2 (ETET + n°); 
t* =k, as + i gent) tł =k: (K? + ix caky) 
a= K(e—Fa)—a i Ke =t pet — Sey 


Z=— kk [23 ek (P+ F (ega) 34 


If we may disregard dissipation, in the region where x is real, Y and Z 
are real, and the exponential factors in equations (9) characterize only 
the wave phase. 


As may be seen from (8), when €] strives to zero (in the case of 
E = 0 €11] = £1) one value of x? goes to a finite limit and the other be- 


comes infinite as 1: 
€ 
1 


a 2 2 BS | 2 w? 2 

al =k teim g) keti). 

1 y wo? 2 1 p2 2 £ z c 
zt z E3 


,=0 e,+0 


In the case of x? = x2,, the solutions of (9) are valid for the point 
where €; = 0. They are not "inflated". The remaining two solutions of 
system (2) for the vicinity of the point where £] = 0 cannot be represented 


in the form of formula (9). When trying to determine them, one may assume 
that £2 and £3 are constant, and for £] one may restrict oneself to the 


linear components in the Taylor expansion: 


pes , 
E = ex. 


After several cumbersome computations, we obtained the following 
equation from the system of equations (2) for the field component E, 
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, OE, OE, aE 


a w? OE w? í 
ex ga + 4 a — (= s3 + kea) ae + (= Es — k?) (ex aa t 


ðE A 2 2 i 
+ 2e; = = [= e? (= e3 — r?) — Rk? ke, + = keea] Ex =0. 


Disregarding the last two components, we may write 


where c 
2 
PE, ae at Hie 
= oe ’ p = 7 . 


If we substitute u(x) = —~3v(E), x = &%, equation (10) becomes a Bessel 
equation 


d?v 1 dv 9 
ate: a (+ p) =o 


and, employing the recurrent formulas for the Bessel function, we obtain 
„E, = AMV bx) + BK, (2V px) (11) 
V px 


where A and B are the integration constants; I1,(2V px) and K} (2Ypx) —- 
the McDonald and Bessel functions. In the region under consideration, 


the components of the field E are related by the following relationships 


CE . OFF @E OE. 
a Thy zi oat ike, 
x Ox Ox Ox 


which enable us to determine Ey and Ez: 


gi = gia [An (2V7) — BK, (2Vpz)}. (12) 


Ss 
N 


Taking the asymptotic form of the Bessel function and requiring that ex- 
pressions (11) and (12) change into expression (6) for large values of |x|, 
we may find the relationship between the coefficients A and B and the co- 
efficients C included in (9). In the case of p > 0, the coefficients C 
are related to A and B in the following way to the right of the "inflation' 
point (x is imaginary in this case): 
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1 wej A 1 wre, \ "3 iA —: 
C} = — (5 3) =~, CL== i A oe 
(—i) TC P j/: \2nc p 
where C = Cy in the case of x imaginary negative; C = C_ in the case of 


x imaginary positive. 


2. If the magnetic field is perpendicular to the layer surface 195 


(E = ree the system of equations (2) may be reduced to a system of two 


second order equations 


da w? . wo? 
wt (æ e1 — k?) g + ieat = 0 
Oe, 
520% >; (13) 
d? 3 d w? - 
es T — ° ae + (a % — k?) ap — iene) =0 | 
T Egs — k J 


where 1 i i 
Ey = pa (kup + hp) Es = qe (bop — kug); Ex =~. 


z ts — k? 


Outside of the vicinity of the "inflation" point, and representing $¢ 
and y similarly to (6) in the following form, we have: 


x 
i§ xdx 


9, p = Ay, p (x)e 9 , 


We obtain the following from the system of equations (13) in the approxi- 
mation of geometric optics 


o? k E 1 hp) ekg eo Ee. es TOL, T 
WP = ey aei +2) +N ktl — e) + 45 ele (= es — k?) ; 


A =i T eC Q's; Ay = K?CQ'+, (14) 
where zi 
-T €g — k? “ 


2 
e : 2 2 è 
v? appo | K? =x +k — a fp 


OT fee Fa). 
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C is an integration constant corresponding to one of the four 
functions x(x). 

As may be seen from (14), close to the point €3 = 0 (in the case of 
= a5 €11 = €3) one solution for x* is finite, and the other solution has 
the following form 


a=k, 
E3 


In the first case, solutions of (14) are valid at the point £3 = 0. 


They correspond to a "non-inflated" wave. In the solution of system (13), /96 
we shall assume that ¢€, and €9 are constant, and we shall expand £€3 in 


Taylor series in the vicinity of the point £3 = 0, and shall confine our- 
selves to the linear component 


We may find the following from the first equation of system (13) 


E E (15) 


dx? wre, © dx?’ 


Substituting this expression in the second equation of system (13), 
disregarding the small components, and making the substitution v = 


=- s4 » we obtain 


1 
rg ae 


where 


ke, 


p=— >. 
3 


When € = 2Y px is used as the independent variable, this equation 
is reduced to a Bessel equation 


dy , 1 do = 
ee qe te =0. 
Utilizing (15), we find the following from the latter equation 
p = C,H (2V px) + CoH (2V px), (16) 


where C}, Cy are the integration constants; HG) (E), u (3) (€) -- Hankel 
functions. 
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As may be seen from (15), close to the "inflation" point 9$|%|e3y| << 
<< |p| -- i.e., the contribution from p to the expression for the fields 
in this region is negligibly small. If we make the stipulation that (16) 
change into an expression like (14) (outside of the vicinity of the "in- 
flation" point), we may find the relationship between the constants C2, C2 


in (16) and the constants C in formula (14). 


Just as in the case of an isotropic layer, the presence of "inflation" 
points in a magnetoactive, tionuniform, transmittant medium leads to the 
fact that the thermal radiation intensity of this medium is on the order /97 
of the radiation intensity of an absolute black body, if the wavelength 
being studied is comparable to the distance from the "inflation" point 
to the layer boundary. When an electromagnetic wave, whose length is 
comparable to the plasma layer thickness, falls on this layer, the amount 
of energy absorbed by the plasma in the vicinity of the "inflation" point 
per unit of time (per 1 cm^ of layer surface) is 


S~ iF? (17) 
n > 


where E is the electric field amplitude in a vacuun. 


However, (17) is valid until the phenomena caused by the oscillation 
nonlinearity and spatial dispersion in the vicinity of the "inflation" 
point exceed the phenomena caused by particle collisions in this 
region, i.e., 


loyal lw; KV Ls», 


where v is the electron velocity caused by the wave field; v -- the effec- 


F > š Ww 1 “3 
tive frequency of electron collisions with ions; K Pigs or od effective 
wave vector in the vicinity of the "inflation point; £ -- layer thickness; 


m -- electron mass; T - temperature. This leads to a limitation for the 
field E and the temperature T, which can be fulfilled by (17): 


Esl H: TS aL?esngi\ 
wc “9 mo , 
where Hg -- the constant magnetic field strength; L -- Coulomb logarithm; 


Ng -- electron density in the plasma; e -- electron charge. The maximum 
energy obtained by a charged particle per unit of time is 


v8lH2 


W ~ Toten, ' 
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In the case of ng = 10!5 cm-3, T = 106°K, Hg = 105G, w ~ 5.1020 
sec-!, w % 100 ev/sec, E < 1 v/cm. 


Thus, it is not advantageous to employ "inflated" fields for plasma 
heating. 
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EXCITATION OF A MAGNETOHYDRODYNAMIC WAVE 
GUIDE IN A COAXIAL LINE 


S. S. Kalmykova, V. I. Kurilko 


As is known (Ref. 1), magnetohydrodynamic waves have low frequencies 
(w << y ES 1.57108 -Ł, with respect to Avac >> 102 cm). At these fre- 


quencies, a coaxial line represents the most reasonable method for trans- 
mitting energy from an oscillator to a wave guide. Therefore, the problem 
of excitation of an axially symmetrical E-wave in a magnetohydrodynamic 
wave guide by a TEM wave of a coaxial line is of great interest. If the 
plasma temperature may be disregarded, and its conductivity is great enough, 
under these conditions, the plasma may be characterized by the dielectric 
constant tensor without spatial dispersion. The problem under considera- 
tion is then a special case of the problem regarding the matching between 
an anisotropic dielectric wave guide and a coaxial line. Equations for 
determining the stray field were obtained in (Ref. 2, 3): 


w+ te zezena Cp a TOT 
O+ TET S ZORO | Zw Rn] rr a) 


S k Z(t) Z(t) dt’ 
=~ ia Y ZE) + Za UN (E k (1) 


; Imk >0. 
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The diffracted magnetic field in the space between the wave guide 
(radius a) and the casing (radius b) is expressed by means of the in- 
determinate function $(t) = ¢t(t) - $7(t) = ot(t) - ot(-t) (the indices /99 
(+) designate the analyticity in the upper or the lower half plane of 
the complex variable t): 


+o 
H(z) = {¥- tae (t) edt npa z < 0; (2) 
te, Ph 1 1 
= + 
Ho- (3 Rw morro Orr O- 


x 3} exp (itz) dt for z>0, (3) 
where 
; 1 
+A — xt (—) = g Ot aaa 


_ 2 | bola), BY) Io(Ba) 
20 =F na LOS na E 


An (F) = Ia (0r) Ko (vb) — (—1)" K, (vr) Io (vb), n =0,1; 
o = (£? — R22; i B= [ze = eel” 


If the plasma density is large (we a2 >> c?), equation (1) can be 
solved according to the iteration method. 


The reflection coefficient of a coaxial wave in the first approxima- 


tion, as a function of the parameter —77>, is 
3 


R=- zmt |! zn] a1 m 


ha Sahn (b/a)? >i. 


When the amplitude of three-dimensional waves excited in a magnetohydrodynamic 
wave guide is computed, a distinction must be drawn between the four regions 
of the plasma parameters and the wave parameters. The magnetic field 
amplitudes have the following values corresponding to these regions (the 
amplitude of the coaxial wave. equals unity): 


0233 


Ta + In-?(b/a))-', —1€4 Ki 


2 
= mea a 
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22 792 2)2 
Ta S a Mikeli; 


~ 32 tn (bay L@ eS 
3 wy 2234 
eee EN a "1 Keg: 
Soe Ue Zl Dae tile ») 100 
2 nae 


[_ = >t, 

* 226, In (6/a) a> 
Thus, the effective excitation of magnetohydrodynamic waves by a 

coaxial line is only observed for strong magnetic fields (ej - 1 << 1). 


Harmonics with large numbers corresponding to the limiting case i >> 


>> a*/52 are always only slightly excited. For weak magnetic fields, the 
main portion of the coaxial wave power is dispersed into excitation of 
the surface plasma wave. 
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THEORY OF MAGNETOHYDRODYNAMIC WAVE SCATTERING 
AT THE END OF A WAVE GUIDE 


V. I. Kurilko 
The study of magnetohydrodynamic waves is of great interest in 
solving several problems of plasma physics, such as controlled thermo- 


nuclear synthesis, magnetohydrodynamic oscillators, etc. (Ref. 1). A 
great many articles have been recently published which investigated the 
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propagation of magnetohydrodynamic waves in uniform, unconfined magneto- 
hydrodynamic wave guides. However, practically every wave guide is con- 
fined. Therefore, it becomes necessary to investigate the phenomena re- /101 
lated to the scattering of magnetohydrodynamic waves at the end of a 

wave guide (for example, reflection of one of the eigen waves of such a 

wave guide, and its transformation into other waves). In addition, the 

theory of electromagnetic wave scattering includes the excitation of a 
confined wave guide. 


This article investigates these phenomena for a semi-infinite 
plasma wave guide. In the general case, this problem may be reduced to 
a system of two coupled integral equations or (in the presence of an 
infinite casing) to two infinite, coupled systems of algebraic equations, 
whose solution may only be found by numerical methods. Therefore, let 
us investigate the case when the end of the wave guide is covered by a 
conducting diaphragm. As will be shown below, the problem of determining 
the Fourier component of the scattered field may be reduced to an integral 
Fredholm equation of the second type, whose solution may be obtained by 
a numerical method, and when there is a small parameter -- in analytical 
form. 


We shall assume that the plasma conductivity is infinite, and its 
temperature equals zero. In this case, the electrodynamic properties 
of the plasma wave guide, as is known, may be characterized by the di- 


2 2 
w (A) 2 
electric constants ej = 1 + —, €= 1- zo w2 = Sten = Ho 
wy ay we 0 m mc 
eo 


= 7 >> o) » which only depend on frequency. Since the frequency is 
c 
fixed in our problem, the plasma wave guide may be regarded as a speciai 
case of an anisotropic dielectric wave guide. Therefore, let us first 
investigate the more general problem of electromagnetic wave scattering 
by the jump in the dielectric constants of an anisotropic dielectric wave 
guide, whose uniform sections are separated by the conducting diaphragm. 
The plasma wave guide parameters may be employed to find the analytical 
solutions of the general equation obtained. 


Thus, let us investigate an anisotropic dielectric wave guide 
(r < a, -> < 2 < +”) with a piecewise-uniform tensor dielectric constant 
A ae 
e(z>0, O0<r<aj= (eP, e0); 
^ : vo 
e(z<0, O0<r<a)= (P. e); 


^ 
e(—œ <z < +o, a<r <b)=l. 
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Let us assume that the uniform sections of the wave guide are /102 
separated by an ideally conducting diaphragm (z = 0, 0 < r < a), and 
that the conducting casing is not confined and is located at the distance 
r = b from the wave guide axis. Let us assume that one of the eigen, 
axially symmetrical E-waves of this wave guide, which is characterized 
by the wave number h,, falls on the nonuniform section from the right 
wave guide. Let us determine the amplitudes for the eigen waves of 
both wave guides which are excited due to scattering of this wave. We 
may write the solution for the fields in the following form 


nE +o ’ 
He = $ H (tyA, (t) edt + exp (—ihmz) r=a+0, 
Et= (11, Wd, (etdi + ZL exp (ihng) | Z2 <?< 
z k 0 m J <+ o; (1) 


An (t) = In (va) K, (vb) — (— 1)"K, (va) Io (vb), n =0,1; 
v (t) = (0 — Ry k =>, Im o > 0; 


+o 
H= $ hs (£) J, (Bsa) eff2dt + 2A, cos hmz 


r =a — 0, 


+œ 
E= ea (£) hs (t) J, (Bsa) e#dt + 2A,Z' cos hm, s=1—z>0, 
4-00 s=2—2<0, 
pics (h (£) J; (Bsa) ef? + 2A Am (sin Amz) . (2) 
r : kels) S 1 s'em in keD 


Žo i! 


where 


l B (O Jo (8,2) i (9) Ha 
Zs (£) = A 7 Jy (Bsa) i Bs Ea [5 (eR? m £) , 


A, =l; A, =0. 


The fields thus selected satisfy the Maxwell equation. We can deter- 
mine the remaining, unknown Fourier amplitudes H, h, of the desired fields 


from the boundary conditions on the lateral surfaces of the wave guide 103 
and on the conducting diaphragm. The boundary conditions on the wave 
guide -- vacuum surface have the following form 
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HP =He, EW = Re for r=a, z>0; 
: H® = He, E®9 =E; forr=a,z<0. 3) 


Substituting (1) and (2) in the boundary conditions (3) and employing 
the results derived in (Ref. 2), we may express the unknown functions 
H and hg by the boundary values on the contour Im t = 0 of the functions 


which are analytical in the upper (+) and the lower (-) half-planes of 
the complex variable t: 


1 _— 
hi (H) Jy Ga) = py {2 09+ (t)— ot 0-3 ; 
1 7 
mAn = p ZOTO- O gee (4) 


D(t) 
ıp —z! 
H(t) A (t) = D5 fz (E) pt (t) — p+ (£) -A 


—_ 1 Oink — zZ — Zl 
=p (tx (t)— E 0 — FE 


= 


where 
D: (=Z) — Z0), ZO =F EH. 


The latter equation (4) represents a boundary problem for determining 
the unknown functions qt, yr, —& and x”. The relationships lacking between 
these functions may be determined from the boundary conditions on the con- 
ducting diaphragm 


Es(z =0; 0<r<a)=0. (5) 
Substituting E$ in the boundary conditions (5), we obtain 
hs (—t) = As (t). (6) 
We obtain the following by means of the latter equation and the Sommerfeld 


condition for the finiteness of the magnetic field and integrability of 
the electric field close to the diaphragm edge: 


OS —¢t (i); TOS— pt (-O, (7) 


where /104 
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— — — - Zm 
TO = O+ a at TOSTO + aga 


1 s 
O=O Ery O=O a TER 


Im hnm > 0. 


(Equations (6) are only necessary for fulfilling the boundary conditions 
(5). We shall assume that these equations are sufficient, although this 
may only be proven for a rectangular wedge [Ref. 3]). 


Thus, the boundary problem for determining the function assumes the 
following form 


| Qty Z(t) — Z. 
Daa ZO (4 Fe. OT Fal 


(8) 
1 


at Mh, 22 (t)—Z) 
=pig|-2 er 0 +4 — ar ar i 


This type of boundary problem for one special case (ey El, E> > œ) 


was first studied in (Ref. 4). It was shown that in the presence of a 
casing, when the coefficients have singularities of only the pole type, 

it can be reduced to an infinite system of algebraic equations. However, 
by employing the formulas of Sokhotskiy -- Plemel', it is more advantageous 
to reduce the problem (8) to an integral equation. By combining and sub- 


tracting equation (8) with its mirror image at the point t = 0, we obtain 
(Ref. 5) 


we 
bi) at! zoff AO e 
a ia +h2oyeu o Ier ) Zø E i 
s=l 
l l (9) 
= Z(t) Wai— 
=— ot. aa PO Pn = ft) Z(t) 91 (£) = $1 (£). 


The index of the latter equation equals zero (Ref. 5). Therefore, it is /105 
is equivalent to one integral Fredholm equation of the second type for 
vice) = vice) - yt(-t): 
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+a 2 +- 
Z(t") de’ QV z5 (¢") 1 1 pog 


1 
nO +a f Ft LD, C | ZO Z7@)| Sr 
— o s=1 —o 


(10) 
ie dt’ 
maf ©. ” 


~ 00 


According to the general theory of singular equations with a kernel 
of the Cauchy type (Ref. 5), the solution of (9) which vanishes in the 
case of t > æ exists and is real. Due to the equivalence between (9) and 
(10), the same holds for the solution of equation (10). It can be solved 
numerically in the general case, and in the presence of parameters -- it 
can be solved analytically, even if there is no casing (b > œ) and the 
coefficients in equation (8) have singularities of the point branching 


type. 
In the case of the magnetohydrodynamic wave guide, which we are in- 


vestigating, the ratio between the wave guide radius and the wave length 
in a vacuum serves as one of the large parameters. Even for Hp ~% 10' gauss, 


Ao = a ~% 10? cm. Therefore, in the case of a < 5 cm the ratio a/d (A -- 


wavelength in a vacuum) cannot exceed 107? (with allowance for the require- 
ment that A >> Ag), so that even In (A/a) is large (In (A/a) 2 4). In 


addition, in a significant number of important cases the linear plasma 


. ô Co 
density na? is great, so that the ratio a bee is small. 


Assuming that the inequality 6 << a << à is fulfilled, we can signifi- 
cantly simplify equation (10). Disregarding terms on the order of /106 
3 7 P AR 
(6/a)In(A/a), x) In(\/a) and taking the fact into account that €; = 1, 


we obtain 


ay Pt Ye ZOU HH ee 


fo $o 
le veya’ ZOA yea hm Z(t) Zm 
Port z f: + ZH (10") 


(the index «<1 > for y is omitted from this point on). Equation (10'), as 
may be readily seen, is equivalent to the scalar boundary problem 


An 1 a 
yt) = zofa OSa aali + zal ey 


The solution of (11) has the following form 
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tin OHA y4 dt hZ) z} qe — kh 
y(t) =a X ofS : a(t? a8) 1 +705 Xt ey)? (12) 
c ' 


where 


+ __ Ko (va) y— pay. + a 1 
Xt =X: X*O=— 


and the singularity at the point t" = t passes around the contour C from 
below.. 


2 
Within an accuracy of a term on the order of (2) In à << 1, we 


have 


; T xa 1 
X+ (£) = (ta + kay op i eae FFA X 
ka 


x In [in (ž-imż)]}. 


where y is the Euler constant, and in the case of In à >> 1 we have 


(12') 


nha a Gs : 1 In(. + t/k) 
x D=|3 in] (ta + karfi — 7: Inà/a }. 


Expansion in powers of 


-1 
In x is usually employed in antenna 


theory (Ref. 6), and numerical integration may be employed to determine 
expression (12') for S4-(t) which are not too small. 


In the same approximation (6 << a << A), the boundary problem for de- 
termining the function ¢t(t) has the following form 


He fists Mn 1 Zh 
t= (Oo ae l+7H . (13) 


The integral in the solution of (12) may be readily computed by closing 
the integration contour in the lower halfplane in the first term, and 
in the upper halfplane in the second term, by the complex variable /107 
t. The integral in the solution of (13) may be computed similarly for 
¢+(t). Substituting the expressions thus obtained for yt(t) and gt(t) into 
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the formula for the amplitude of the Fourier-field in a vacuum 


6 OHO = pm? MOOI So, (14) 


by means of (1) we may find the expression for the amplitude for a wave 
with the number 1, which is excited in a wave guide (z > 0) during the 


scattering of a wave with the number m: 


l 1/2 
Gi [aane et Mant 2zina | et 
== : G 


Mm H 2 
nm t1 \ PXA i 
+ ae + 2Zinttm a+ 
nm (Bm — u) my +1 xF nê — n? 
l mn+1\ xE z} — z} (15) 
Err T E + mgr — 
Ny (Am — 1) fm tl Xa ni — nn 
ziz ia 
ee ee nS, 
where Dux Gq Dalim mi Zi iy =e (h, mË Xt n= Xt (ty, m); 


g = gt (h); 4; = yt (h) . 


The latter expression may be significantly simplified in the most in- 
teresting cases of small and large retardations ngm’ It thus appears that 
the sum in the parentheses,in (15) is represented in the form of a power 
series of the logarithms for the ratios A/a and à/ô. Retaining the old 
terms in this series, we obtain 


2 
2ln — 
Rigg oe ai 14“ “BAK 
2in2? — —. 
Maloy FPP by (16) 


Ji Oi m) =0; n? =e; 


ay 2 inçe, ~1 
Rn = an] Ped (leper ae 


2 —_-_ 
In a(e, we 1) 
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EER ea Te (17)  /108 


B7y3, 


em À 
a K In: K a (18) 


pees E eT L 
Ym byi m n at n by, 7 bv, 


l = 1 [a « 
m VE (itm) Ea ? 


jas + paan 
Ing in by, 


2 2 
Kim <p (19) 


It can be seen from these formulas that the coefficients of the inter-trans- 
formation of the magnetohydrodynamic waves Rin decrease with an increase 


in the wave number and a decrease in the magnetic field strength. If the 
plasma wave guide is surrounded by a conducting casing r = b, -~ < z < +>, 
it is impossible to study magnetohydrodynamic waves in the case of à >> b, 
and only the transformation of one oscillation into another occurs at the 
end of the wave guide. The transformation coefficients for this case may 
determined by substituting Xt (t) from the solution of the corresponding 


problem X + (t) = [vatn] X - (t) and the wave numbers hg m -- from the 
solution of the dispersion equation, in expression(15). 


i (h? — k? 
Paea E ht o; 
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DETERMINATION OF PLASMA TEMPERATURE AND DENSITY 
DISTRIBUTION BY REFRACTION AND DAMPING OF A BEAM 


V. L. Sizonenko, K. N. Stepanov 
Let us investigate the passage of a beam through a nonuniform plasma. 
We shall show that, by changing the angle of incidence or the frequency of 
the microwave signal, we can determine the plasma density distribution 
by the beam refraction, and can determine the plasma electron temperature 
distribution by the damping of the wave energy along the beam. 


In order to determine the beam trajectory in a plasma, let us employ 
the Fermi principle 


è f kds =0, (1) 
a 
where k is the wave vector; ds -- an element of length along the beam tra- 
jectory. In an isotropic plasma, we have 
=1Y =o, (2) 
2 
where w is the wave frequency; wy (x) = Arale) -- plasma Langmuir fre- 


quency; n(r) -- plasma electron density. Modulation is performed in 
formula (1) for fixed ends (8x; |a = 6x4 |b =0. Taking the fact into account 


that ds = Y g;,dx;dx,, from expression (1), with allowance for (2), we may 
obtain the equation which determines the beam trajectory for given wo (£) 110 
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df Voi Ot) 1 Va wk i Vo 
Alen ot — 0b Ft) 4 Vo 0 “Gy, as ds zV o w = 0. (3) 


By employing equation (3), let us investigate the inverse problem: 
Let us find wo(t), i.e., the plasma density distribution based on the beam 
refraction, as a function of the angle of incidence for a fixed frequency, 
or as a function of the frequency for a fixed angle of incidence. 


Planar Problem 
If the density of a plasma filling the halfspace x > 0 depends only 
on one coordinate x, and it increases monotonically with an increase in x, 


then the beam trajectory in the case under consideration is flat (Figure 1). 
Equation (3) may be written as follows 


g 
ATE = 0; ds = 


= V dx® + dz?. 


(4) 


Integrating equation (4), we obtain 
Vo — oF =a = const, (5) 


where a = w sin y, y is the angle of incidence. 


The position of the point at which the beam leaves the plasma z = & 
can be determined from the plasma density distribution; it depends on the 
angle of incidence and frequency. If the dependence £() is known from 
experiments, we may obtain the density distribution n(x) in the case of 
x < x*, where wg (x*) = w (this problem is similar to the problem of deter- 
mining the potential energy by the specific dependence of the oscillation 
period on energy [Ref. 1]}. In actuality, it follows from equation (5) 
that 


dz _ sinp 
dx we g 
cos? p — z (6) 
We thus h ie i 
e us have Ky) =2 sinddx = /111 
2 25 
cos“ bp — z 
7 (7) 
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where x = x9 is the rotation point of the beam, determined from the 
following condition 


cos? ) = 2n 
2 
wo 
Introducing a new variable u = ~; instead of x, from formula (7) we 
u) 
obtain cosy 
nm 


sin ye du 
LW) =2\ Vospua' 


0 (8) 
It thus follows that 
1 1 cos 
L(Y) d cos? Lae 
cos? b ðu 
— 2 \ dcos? = 
sinp V y= a TF ? V lost} — u) q — cos*t) 
ó 0 ô 
1 
= + 2a | ude = + 2nx(1). 
0 
2 
WO 


Assuming that y = “52 we finally find that 
w 


2 

20 

oF 
l _ Lpdcosy d cos? } 


: (9) 
gem a — cos*p 


Equation (9) determines the function x(wo), i.e., the dependence of density 
on the coordinate in implicit form, according to the specific dependence 
&(y). We may employ this equation in the case of a cylindrical plasma, 
when the plasma is probed by a beam in the plane $ = const passing through /112 
the cylinder axis. Thus, wọ = wo(r) and x = r = R, where R is the plasma 
radius (wọ(R) = 0). By performing similar measurements for different 
values of the aximuthal angle $, we may obtain the density distribution 

in the case when density depends on the angle $ (but does not depend on z). 


The plasma density distribution may also be determined by probing the 
plasma with microwave signals having a different frequency. For a given 
angle of incidence, the position of the exit point of the beam z = £ de- 
pends on the frequency w. Knowing the function 2(w) from experiment, we 
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P 


Figure 1 


may readily determine the function wo (x) according to equation (7). 
follows from equation (7) that 


1 
__ Le) de oo ede 
a = 
tgy tev V 1i — oi m 2 qx (y). (10) 
Ges a ( Di Vocale] 


— we obtain the expression determining the dependence 


Assuming that y= Goa 


wo(x) in implicit form 


ere 


[o)do da , (11) 


Zi ig Ft) eyo wZ/cos? p — a? 


We may also obtain formula (11) from the expression given in (Ref. 3) 
[see also (Ref. 4)] for the real layer height x(wọ) according to the specific 


effective height Xp (w) for the case of normal wave incidence on the layer, 


with allowance for the following relationship (Ref. 4) 


1(%) = 2tg bx,(w cos ¢). 


Cylindrical Problem 


If the beam passes through a plasma cylinder in a plane which is per- 
pendicular to the cylinder axis (Figure 2), equation (3) may be represented 


in the following form 


d 3d dr? + rdo? 
a (PV 0? 00 g) = 05 ds = Va Frag (12) 
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¥ 


Figure 2 


The plasma density distribution n(r) may be found from the specific 
experimental dependence of the angle at which the beam leaves the plasma 
$9 on the angle of incidence y * (this problem is similar to the problem 
of determining the potential energy from the specific dependence 


£113 _ 


of the scattering cross section on the scattering angle (Ref. 1, 2). Inte- 


grating equation (12), we obtain 


2 
dr r? v R? , 
E E (13) 


It thus follows that 


Pol?) = 2R sing a 
m TE Boe Ea 


where r = rg is the rotation point of a beam which is determined from the 
following condition 


o (ro) 


z + moi y= l. 


w 


Assuming that 


* This problem was solved in (Ref. 5). 
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2 
R @, 2 
't=lnž; eale la; 
we can write equation (14) in the following form 


sin? 4 
Jú du 
Pol) = 2sing f aa 5) 


from which we find that 


-A t sin? @ 


Or 
go (p) d sing of sin dd si du t“ 
-oe = in = n - 
Js sin? = 7 J $ y Yy (sin? m Pe m (— sin? p) TT (1) 
w2\ a 
Assuming that w=p=(1— Ze. we may find the dependence r = r(wg): 
j 1 £ 
r=Rexp| _1( eo @dsing 
n a . 


sin? ¢ — u (16) 
Va 


Formula (16) determines the function r = r(wo), i.e., the distribution /114 
of density over the radius in implicit form, according to the specific 
dependence of the beam exit angle on the angle of incidence. 


Performing measurements in different planes z = const, we may thus 
obtain the density distribution along the cylinder axis in the case of 
axial symmetry. For a fixed angle of incidence y, equation (15) determines 
the beam exit angle ¢) as a function of frequency. A knowledge of the 
function ¢9(w) enables us to determine the density distribution. We thus 
assume 


w (r) 
Q (r) = - 7 , 
y: a sin? p 


t = arc sin (F sing). (17) 


Equation (14) may then be written as 
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Po (9) 02 
2% Vor — g’ 
We thus find 
r E 2%% ao 
Po (0) dw = 0Q Besa as PE 
Seca = eo (pases — 7M — FO. 


Since &(0) = Elona = ~, assuming that y = 2, we can write 
e 
1 
Eoso 5 | eee. 
0 


We thus obtain the expression determining the dependence Q(r) in implicit 
form 


jee Rimni 
1 {2o (0) de (18) 
singy + > == 
| ja] 


Temperature Determination from Beam Damping 


Due to collisions, the energy of a beam leaving the plasma is 
eT times less than the incident energy, where t = f x ds is the 


optical plasma thickness; x -~ damping coefficient. If the wave fre- /115 
quency is considerably greater than the frequency of collisions of elec- 
trons with ions and neutral particles v(r), we then have 


o? 
ee as 


2 
cw? / 1—0 2 
W 


Knowing x(r), we may find the plasma temperature distribution from the 


specific dependence of v on temperature. 


The damping t depends on the angle of incidence y. We may find the de- 
pendence of the damping coefficient on the radius from the quantity t(y), 
which is measured experimentally, and from the specific density distribution. 


Let us first investigate the planar problem (in the case of a cylin- 
drical plasma, the beam trajectories lie in the plane ¢ = const). For 
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t(~), we may write the following expression 
Xe 


swej y EF t 


0 


Substituting gz from formula (5) in this, we obtain 


costp 


(H) = afao Vo pdx = 2 | x je (19) 
0 


Let us introduce the function f(u) according to the equation 


(20) 
d 7— ad 
T xu) yT uE 
We then have 
cos? ẹ 
f of du 
ðu 
=2 n 
(4) 7 y costy—u (21) 
Solving this equation and equation (8), we find 
A òr ; 
aja) 1. | acosg cos Y 
d m) Vicos ` 
v 
2 
WO 
Assuming that A = w2? we obtain 
y= tay Ee ee (22) 
a(x) = = \ 7- es 
at e costy 


In explicit form, this formula determines the dependence of the damping 
coefficient on coordinate x. 


If the beam trajectories lie in the z = const plane, we then have 
R. -c 
2- /dẹ\? 
“ =2f “inf ite (Z) dr. 
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Substituting the value of ag from formula (13), we obtain 


dr 
R 
wo? — og 
(9) =2 \ x(r) E T ET.: R ag dr. (23) 


Assuming that 


r? o? di 2 P 
"e ee T) “ag 1-3-4 


we may write equation (23) in the following form 


df du 
(Y) =2 sio 
; ; „yesin ọ (24) 


Solving it in the same way as equation (19), we find 


1 
yoo afa NY 
Vsin? pr 


m 


Assuming that y = u(r), we obtain the following expression 


1 
ar 
(y= 2 (2 | asap TY 
x | dr Vsintg—u a ( =) (25) 
Vu u = Ri Le 


Expression (25) determines the dependence of the damping coefficient 
coefficient on the radius according to the specific dependence of T on jy. 
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SECTION III 


PLASMA NONLINEAR OSCILLATIONS AND WAVE INTERACTION 118 


EXCHANGE OF ENERGY BETWEEN HIGH FREQUENCY AND LOW FREQUENCY 
OSCILLATIONS IN A PLASMA 


V. D. Fedorchenko, V. I. Muratov, B. N. Rutkevich 


Oscillations with frequencies v on the order of 100 ke (Ref. 1), 
as well as oscillations at the electron plasma frequency (wpe % 50 Mc), 
occur in a plasma produced by an electron bundle in a longitudinal magnetic 
field at a pressure of 10-6 — 10-5 mm Hg. It has been found that there 
is a relationship between the low frequency and high frequency oscillations 
(Ref. 1 - 3). High frequency oscillations cause an increase in the low 
frequency oscillation amplitude. On the other hand, they undergo ampli- 
tude modulation at a low frequency, which leads to the occurrence of com- 
bined frequencies w + v. The amplitudes of oscillations having the fre- 
quencies w + v and w - v are not the same, which -- just as with the 
Landsberg-Mandel'shten-Raman effect -- can indicate the direction and 
effectiveness of energy transfer from oscillations of one frequency to 
oscillations of another frequency. 


As has been shown previously (Ref. 1), low frequency oscillations /119 
represent ion motion (which is transverse with respect to the bundle) in 
the field of the space charge of electrons contained by a magnetic field. 
The frequency of these oscillations may be readily computed from the 
following simple model. Let us assume that electrons and ions fill one 
and the same cylindrical region in the equilibrium state. The simultaneous 
shift of all ions by the quantity i> and of electrons by the quantity Že» 


leads to polarization, which can be determined as follows in the case of 
a small shift 


ESE r) T 
where n is the plasma density. 
We can write the following expression for electrons and ions 
= + 2 (rh: (2) 
na É (3) 
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or, with allowance for (1), 


ae Bo toy abs (4) 
ri = — obi (fi — re) + = [rB] 
i 
= ee am (5) 
Te = woe (i — re) — = [r.B], 
i e 
where 
2 ne. 
Oe gym,” (6) 
w2 — ne 
Oe gm (7) 


Disregarding the electron shift during orbital motion, we shall only 
allow for their drift in crossed electric and magnetic fields, and 
we shall designate the position of the electron guiding center by the 
vector ae According to expression (5), we then have 


> 


Ti— Te = < [rB] = S [rB]. (8) 
e’0e Moy 


It follows from formulas (4) and (8) 120 


~ ect +> 
n= — p eB] + a [r:B}, 


and thus 
~> e > — e — e? ss > 1 B32 % 
tı = ——[r.B) + —[7,B] = ERA 
i im, | 2B) + m, [r:B] mal, (ir. B] B] moe, Te (9) 
or R a 
le = Fr = — Qz’, 
Hi 
(10) 
where 
oH 
SS (11) 


-- the ratio between the ion cyclotron frequency and the plasma frequency. 


Substituting the expression obtained for te in formula (4), we obtain 
ri — ia [rB] + 03, (1 + 257) r =0. (12) 
Equation (12) has a ty ~ elvt type of solution, and the frequency 
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must satisfy the following equation 
Q? + Qul — 1—27 = > 


(13) 
i.e., 
es Yita (14) 
In the case of 97 X l, v X woi vy. 
If Qy >> 1, equation (13) has two roots: 
Q, = —Qy (15) 
2% = z- (16) 
If y << 1 Q= tg 
(17) 


The oscillation frequency v decreases with an increase in the magnetic 
field strength, which coincides with the solutions of (16) and (17), if 
it is assumed that the frequency woi depends slightly on the magnetic field 
strength. If woi does not depend on B -- i.e., Opi is constant -- it may 121 


be readily determined according to any pair of measured v and B, which 
would ‘then enable us to compile a graph showing the dependence of v on B 
according to formula (14). The dependence of ion oscillation frequency on 
the magnetic field strength is shown in Figure 1. The curves were re- 
corded for different current values in an electron bundle (1 - I = 10 ma, 

2 - I = 30 ma). The dashed curves are drawn through the points computed 
according to formula (14), under the assumption that the plasma density 

does not change when there is a change in the magnetic field strength. 

The computed points do not lie on the experimental curves (particularly 

in the region of small fields), which is no doubt related to the variability 


of Woi: Nevertheless, it may be stated that the nature of the dependence 


is correctly imparted by our simple model. 
Modulation of high frequency oscillations may be due to low frequency 
oscillations of the plasma density. Let us assume, for example, that the 


> 
plasma is located in an external field Felt which is transverse with 
respect to the bundle, and the frequency w is so great that only allowance 
for the electron component shift has any meaning. The external field 


Feit gives rise to a shift and the occurrence of a polarized field 
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> 
Me Tae The equation for electron motion in these fields may be written 


in the following form 


` 


7 =—£(Ee +E) Er. (18) 


We obtain the following expression for forced electron oscillations 


© Belot 
Me (19) 


ab, +o (oy, —@) 7 


> 
r= 


If it is assumed that the plasma density undergoes oscillations at the fre- /122 


quency v << w and the quantity ioe: has the form wa (1 + a cos vt), while 


a << 1, then expression (19) changes into the sum of the oscillations with 
the principal frequencies (w) and the combined frequencies (w + v): 


E aa ee, O T 
5 3 6 2 m, wp, Oo. ig 
mpe] b+ o, (°He— o) h: 4- w2, (°He— | 
ie i(w—y) t (20) 
a e E e 
+ 7 Š m r im. . 


-E E O 
1+ -7 ("He — w) 
Oe 
The oscillations have resonance close to the electron cyclotron frequency. 


This model is inadequate for determining oscillation intensity at the 
combined frequencies. 


The experiments were performed on a hollow electron bundle in a longi- 
tudinal magnetic field with a strength ranging between 200 - 2000 oersted. 
The bundle length was 50cm; diameter -- 2 cm; energy -- 250 v; and the 
current -- 20 - 40 ma. The bundle was located in a metallic tube having 
a diameter of 9 cm. Pressure in the chamber was several units of 1076 mm Hg. 
The interaction between the outer, high frequency field and the low fre- 
quency ion oscillations was studied. The interaction was observed in 
three cases: When the frequency of the outer signal (1) did not coincide 
with any of the plasma eigen frequencies, (2) coincided with the electron 
cyclotron frequency, and (3) coincided with the plasma electron frequency. 


In the nonresonance case, the outer field was transverse with respect 
to the bundle. It was produced between two conductors having a diameter 
of 0.2 cm and a length which was close to the bundle length. The conductors 
were parallel to the bundle axis at a distance of 6.5 cm from each other. 
Figure 2 shows the diagram of the experimental apparatus employed to /123 
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Figure 1 


study the interaction between ion oscillations and the outer electric 
field, whose frequency does not coincide with any of the plasma eigen fre- 


quencies (l -- solenoid; 2 -- electron gun; 3 -- probe; 4 -- hollow elec- 
tron bundle; 5 -~ collector; 6 -- conductors between which a high fre- 
quency field is produced; 7 -- coils connecting the generator and the 


measurement circuit. The resonance curve for the circuit producing the 
high frequency field is shown in the upper right). The capacitor is a 
section of the resonance circuit weakly connected to the GS-23 generator 
and the C4-8 spectrum analyzer. The resonance frequency of the circuit 
was 13.88 Mc. The resonance barely shifted when the bundle was switched 
on. The resonance width was quite large (Q = 50), so that the combined 
frequencies did not go beyond it. The connection between the circuit, the 
generator, and the measurement circuit was selected so that the resonance 
curve was symmetrical. This was important for comparing the intensities 
of combined oscillations with the frequencies w + v and w - v. 


Figure 3 presents typical spectra of oscillations produced when the 
resonance frequency was used (the spectra were obtained under the following 
conditions: Pressure p = 3.4°10-© mm Hg, electron bundle current I = 30 ma, 
anode voltage U =.250 v, effective variable outer field strength v, = 64 v, 
magnetic field strength H: a -- 600 oersted; b -- 780 oersted; c -- 920 
oersted). Lateral frequencies spaced at the frequency of ion oscillations 
may be seen, in addition to the frequency employed. Employing the termin- 
ology used in the theory of combined scattering, we shall call the lateral /124 
lines red (w - v) and violet (w + v) "companions". The relative height 
of the companions (with respect to the carrier height) was 1 - 2%, and it 
increased with an increase in the amplitude of ion oscillations and the 
amplitude of the outer signal. 


The heights of the red and violet companions, generally speaking, 
were different, and this difference depended on the magnetic field 
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strength. As may be seen from Figure 3, when the magnetic field strength 
was 600 oersteds, the violet companion was higher, and for 920 oersteds 
the red companion was higher. At a field strength of 780 oersteds, the 
companions were the same. In order to explain the difference in the 
heights, we should point out that the increase in the ion oscillation 
amplitude when a high frequency outer signal was employed was always 
greater for large magnetic fields. For small magnetic fields and the 

same intensity of the outer signal, the amplitude increase in ion oscilla- 
tions disappeared. Thus, an increase in the magnetic field strength can 
improve the conditions for transferring energy from high frequency oscilla- 
tions to ion oscillations. 


Comparing this with the data given in Figure 3, we may arrive at the 
conclusion that there is a relationship between the companion heights and 
the direction of energy transfer. The predominance of the red companion 
corresponds to the transfer of energy from high frequency oscillations to 
low frequency oscillations; the predominance of the violet companion 
corresponds to the energy transfer in the opposite direction. The direc- 
tion of energy transfer no doubt depends on the relationship between the 
amplitudes of the ion oscillations and the outer signal. In actuality, 
an increase in the magnetic field strength decreases the ion oscillation 
amplitude, which leads to a more effective energy transfer from the high 
frequency to the low frequency for a given amplitude of the outer signal. 


The same result may be achieved in another way: by changing the /125 


outer signal amplitude for a constant magnetic field. With an increase 
in the outer signal, the red companion becomes higher as compared with 
the violet. 


Let us turn to an experiment in which the outer field frequency 
coincides with the electron cyclotron frequency. Figure 4 shows the 
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Figure 3 


diagram of the experimental apparatus for studying the interaction be- 
tween ion oscillations and the outer field at the electron cyclotron fre- 
quency (1 -~ solenoid; 2 -- electron gun; 3 -- probe; 4 -- hollow electron 
bundle; 5 -- collector; 6 -~ volumetric resonator; 7 -- connection with 
the oscillation source and measurement circuit). The outer field was 
produced in the volumetric resonator which encompassed almost all of the 
bundle. The mode H,, was excited in the resonator at a frequency of 

2265 Mc. When the bundle was switched on, the oscillation level in the 
resonator sharply decreased, when the magnetic field strength reached a 
value corresponding to electron cyclotron resonance (Figure 5). The 
width of the resonance absorption curve was primarily determined by the 
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nonuniformity of the magnetic field over the length of the system. Compan- 
ions appeared close to resonance which coincided with formula (20). The 
violet companion was higher than the red companion, which was apparently 
related to the small amplitude of the signal supplied. 


Let us examine the case when the outer field frequency coincides 

with the electron plasma frequency. The outer signal is supplied to a 
grid located at the bundle origin. Measurements are performed by a probe. 
Figure 6 shows the spectra for several frequencies close to the electron 
plasma frequency (p = 4-10-© mm Hg, I = 40 ma, U = 250 v, H = 1000 oersted, 
voltage on the grid u, = 0.1 v. Oscillation frequencies of voltage on the 
grid: a -- 39 Me; b -- 40 Mc; c -- 41 Mc; d -- 42 Me; e -- 43 Mc). It 
can be seen that the interaction is resonant in nature. The heights of 
the companions are large (the total altitude of the carrier is shown in 
the photographs), which points to the effectiveness of the interaction be- 


tween ion oscillations and the outer signal at the electron plasma fre- /126 
quency. 


The amplitude of the signal supplied is 0.1 v os which explains the 


* The effective voltages are always employed. 
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Figure 6 


predominance of the violet companion. When the amplitude of the outer 
high frequency field increases (for a magnetic field strength of 1000 
oersted), the red companion begins to predominate. A similar relationship 
between the companion heights and the amplitude of the outer field is ob- 
served for H = 7000 oersted. With an increase in the red companion, the 
altitude of the principal line (carrier) decreases. If H = 700 oersted, 
a decrease in even the principal line may be observed when the outer 
signal is intensified (from ug = 0.5 v tou = 0.95 v). For a compara~ 
tively small magnetic field (H = 400 oersted) and a significant outer 
signal amplitude, there is very strong interaction which is accompanied 
by the appearance of many combined frequencies. 


These data point to the effective transfer of energy from electron 
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plasma oscillations to ion oscillations. This phenomenon may probably 
be employed to increase the energy of the plasma ion component. 


REFERENCES 


1. Fedorchenko, V. D., Rutkevich, B. N., Muratov, V. I., Chernyy, B. M. 
Zhurnal Teoreticheskoy Fiziki, 32, 958, 1962. 


2. Fedorchenko, V. D., Muratov, V. I., Rutkevich, B. N. Zhurnal Teoreti- 
cheskoy Fiziki, 34, 458, 1964. 


3. Fedorchenko, V. D., Muratov, V. I., Rutkevich, B. N. Zhurnal Teoreti- 
cheskoy Fiziki, 34, 463, 1964. 


DISSIPATION OF PLASMA OSCILLATIONS EXCITED IN A 
CURRENT-CARRYING PLASMA 


Ye. A. Sukhomlin, V. A. Suprunenko, N. I. Reva, V. T. Tolok 


Several experimental and theoretical investigations have studied the 
development of bunched instabilities in a current-carrying plasma for 
large electric field strengths (Ref. 1 - 5). It has been shown that, as 
only the mean energy of the ordered electron drift is larger than their 
thermal energy, intense longitudinal electrostatic oscillations develop 127 
in a plasma. Their energy reaches the initial energy level of electron 
drift usually after several tens of plasma oscillation periods. 


Computations of the multi-flux motion of electrons in a current- 
carrying plasma, which were performed by O. Buneman (Ref. 6), J. Dawson 
(Ref. 7), and Ya. B. Faynberg (Ref. 8), have shown that very intense 
"thermalization" of the plasma oscillation energy occurs, if this energy 
is considerably greater than the electron thermal energy. Thermalization 
occurs due to nonlinear phenomena leading to the transformation of longi- 
tudinal oscillations into transverse oscillations and to their rapid phase 
mixing. 


This process takes place until the energy of the ordered oscillations 
equals the electron thermal energy. It may be assumed that ion heating 
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will occur due to "collective" friction of electrons on ions in the case 
of bunched instabilities. 


Thus, an investigation of bunched instabilities in a high-current 
gas discharge opens up new possibilities for effective plasma heating. 


The studies (Ref. 9, 10) performed detailed investigations of the 
excitation conditions of bunched instabilities in a current-carrying 
plasma, as well as the plasma characteristics in the presence of these 
instabilities. The occurrence of an anomalously high discharge resistance 
and intense microwave plasma radiation was discovered. 


This article investigates heating and containment of a plasma in a 
strong magnetic field, under conditions when bunched instabilities excited 
by "escaping" electrons develop in the plasma. The experiments were per- 
formed on an apparatus representing a rectilinear tube made of alundum 
with a diameter of 10 cm and a Length of 25 pi which was usually filled 
with hydrogen at a pressure of 5-10-73 - 10-* mm Hg. Aluminum electrodes 
were placed at the two ends of the tube; a battery of capacitors having 
an over-all capacitance of 15 microfarads was discharged between the 
electrodes. The battery was charged to a voltage of 30-40 kv. The 
discharge current through the gas amounted to 100 ka with a period of 
9 microseconds. In order to eliminate hydromagnetic phenomena, the dis- 
charge was performed in a strong longitudinal magnetic field (on the order 
of 1.2 tl), at which the Shafranov condition of stability would be ful- 
filled (Ref. 13). In order that the plasma did not touch the walls, a 
diaphragm with an opening which was 80 mm in diameter was placed between 
the electrodes. 


During the first half-period in the discharge,a highly ionized plasma/128 
filament, which was separated from the wall and which had a diameter of 
80 mm, was produced; no macroscale hydromagnetic instabilities were 

tame in s plasma filament. The plasma density changed between 

- 1013 cm-3. The construction of the apparatus and the experimental 
ee were En in detail in (Ref. 9). X-ray and microwave radia- 
tion from the discharge, the current of "escaping" electrons, the over-all 
discharge current, and the voltage between the electrodes were studied 
experimentally. 


Figure 1 shows the following oscillograms: a -- microwave radiation 
from the plasma; b -- over-all discharge current; c -- current of "escaping" 
electrons; d -- voltage between the electrodes of the discharge tube re- 


duced to a single time scale. The oartliograns were recorded at an initial 
hydrogen pressure in the chamber of 2. 10-2 mm Hg, a magnetic field strength 
of 0.64 tl, and a charge voltage of 34 kv. 


Characteristic, inter-correlated oscillations are observed at high 
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electric field strengths; these oscillations are due to the development 

of bunched instabilities. During the initial period, the field strength 

in the discharge center increases as the voltage wave penetrates the 

plasma. With an increase in the electric field strength in the discharge, 
accelerating processes begin to develop. A current of "escaping" electrons 
first appears due to the "tail" of the Maxwell distribution. However, as soon 
as the electric field strength in the plasma begins to exceed the criti- 

cal value, all of the electrons acquire a drift velocity which is greater 

than the thermal velocity, and bunched instabilities develop in the plasma, 

to which the occurrence of epithermal microwave radiation corresponds. 

In this case, a large portion of the directional drift energy of the 

electrons is transmitted to excitation of oscillations, and the current 

of "escaping" electrons sharply decreases. This leads to an increase /129 
in the effective plasma resistance and to a dip on the oscillogram for 

the over-all discharge current. 


The amount of energy contributed by the outer source to the buildup 
of longitudinal electron oscillations may be computed from the additional 
current at the moment an instability develops. For the case shown in 
Figure 1, this energy amounts to 10 kv per particle. The energy of these 
oscillations considerably exceeds the initial thermal energy equalling 
30 electron volts, which leads to effective thermalization of plasma 
oscillations due to nonlinear phenomena. As a result, the random electron 
energy will equal, in order of magnitude, the energy of plasma oscilla- 
tions. Due to thermalization, intense X-ray radiation occurs as a result 
of energetic electrons falling on the target. 
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Figure 2 


Since the transverse component of electron energy has increased con- 
siderably, and the velocity of directed drift has remained the same as 
previously, the condition for excitation of plasma bunched instabilities 
has been disturbed and the plasma has returned to the initial unexcited 
state. Figure 1 clearly illustrates two cycles of such oscillations with 
a period on the order of 1.5 microseconds. This period is probably deter- 
mined by the time required for electrons to change into a state of "escape 
According to the computations of Dreicer (Ref. 11), this period equals, 
in order of magnitude, the time between two Coulomb collisions for an 
electric field strength which is greater than the critical strength. 
Assuming that 30% of the energy of plasma oscillations is "thermalized" 
(Te = 3-103 electronvolts) (Ref. 12), we find that for a plasma density 
of 7-10!* cm-3 the time between two Coulomb collisions is 4.5 microseconds 


The period in which the "heating" cycles are repeated depends on the 
initial gas pressure in the chamber. It increases considerably with a 
decrease in the plasma density. The effective electron temperature of the 
plasma must thus increase, since the total energy transmitted into the 
buildup of plasma oscillations from the outer source changes very little. 


Thus, it would be expected that intense electron heating occurs due 
to the development of bunched instabilities in the discharge. Direct 
measurements of the. electron temperature are of great interest. 


The effective electron temperature was determined by the absorption 
of electron braking radiation in thin beryllium foils located in front 
of a scintillation crystal on the wall within the vacuum chamber. The 
plasma electrons falling on the foil-target are braked in the very thin 
surface layer. Their energy is transformed into braking X-ray radiation 
which, after partial absorption in the foil, falls on the crystal causing 
a flash of light. In order to determine the radiation hardness, without 
disturbing the vacuum it is possible to place beryllium foils having 
different thicknesses in front of the crystal. The light from the crystal 
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Figure 3 


is supplied to the photomultiplier by means of the wave guide. The signal 
from the photomultiplier through the cathode follower is supplied to the 
oscillograph amplifier. 


Figure 2 shows the dependence of the photomultiplier (PM) signal 
intensity on the thickness of the absorber foil. This dependence may 
be employed to determine the electron temperature for a specific form 
of the electron energy distribution function. Maxwell distribution, 
rectangular distribution with the width T,, and Drayvesten distribution 


lead to similar temperature values. This enables us to employ the curve 
shown in Figure 2 for a rough estimate of the electron temperature when 
the electron energy distribution function is not known precisely. 


The curves in Figure 2 were compiled under the assumption of Maxwell 
distribution for three temperatures: 1 -- 1 kev; 2 -- 2 kev; 3 - 3 kev. 
It can be seen that the experimental points correspond to a plasma elec- 
tron temperature on the order to 2 kev. 


Figure 3 presents the following oscillograms: a -- current of 
"escaping" electrons; b -- X-ray radiation from the plasma due to braking 
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of thermalized electrons in the foil. X-ray radiation arises simultaneously 
with the current, and continues for a long period of time after the 

braking of "escaping" electrons due to the energy of transverse motion. 

The maximum quanta energy of this radiation is about 15 kev. This points 

to the effective transformation of the energy of electron longitudinal os- 
cillations into the energy of transverse motion. In the absence of 

bunched instabilities, the electron temperature in the discharge is 30 
electronvolts. 7 


124 


Figure 4 presents the following oscillograms: a -- current of /131 
"escaping" electrons; b -- X-ray radiation; c -- microwave radiation; 
d -- light from the discharge center; and e -- over-all discharge current. 
X-ray radiation arises simultaneously with intense epithermal microwave 
radiation at a frequency close to 


CRIES M 
v=V tR N a 
t 


(wo -- plasma electron frequency; wy -~ electron cyclotron frequency). 
The power of this radiation is approximately four orders of magnitude 
greater than the power of thermal radiation from the plasma at an elec- 
tron temperature of 10% electronvolts. It was shown in (Ref. 10) that 


the radiation is related to longitudinal electron oscillations in the 
discharge. 


All of these statements indicate that X-ray radiation from discharge /132 
is related to the thermalization of bunched instability energy. 


The period of pair interaction for electrons having an energy on 
the order of 2 kev is considerably greater than the time of the process 
being studied. Therefore, we may assume that nonlinear phenomena during 
collective oscillations with a large amplitude play the main role in 
thermalization of longitudinal oscillations. In our case, the Larmor 
radius of hot electrons is a little less than the discharge chamber 
diameter (less than 1 mm). Therefore, X-ray radiation is primarily 
caused by electrons diffusing toward the walls across the magnetic field. 
When the X-ray radiation reaches a maximum (see Figure 4), the lumines- 
cence in the discharge center sharply increases; this is determined pri- 
marily by admixtures dislodged from the walls by hot electrons. There is 
a simultaneous strong increase in the current on the boundary electro- 
static probe, which is executed in such a way that the current upon it is 
determined by the resistance of the discharge plasma across the magnetic 
field. Thus, there is a rapid cooling of the electrons, even 5 - 6 micro- 
seconds after the X-ray radiation has terminated. 


In the absence of pair collisions, strong diffusion and great con- /133 
ductivity across the magnetic field would not be expected. Anomalous 
diffusion can occur due to nonlinear phenomena for a large longitudinal 
oscillation amplitude. However, in the experiments described, the anoma- 
lous diffusion across the magnetic field may be explained by an increase 
in the gaskinetic plasma pressure as compared with the magnetic pressure, 
due to intense electron heating. The decrease in the X-ray radiation in- 
tensity with an increase in the magnetic field strength also points to 
this conclusion. 
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DAMPING OF INITIAL PERTURBATION AND STEADY FLUCTUATIONS 
IN A COLLISIONLESS PLASMA 


A. I. Akhiyezer, I. A. Akhiyezer, R. V. Polovin 


When the perturbation of the distribution function is a nonanalyti- 
cal function of velocity, the damping of plasma waves when there are no 
collisions may be different from that determined by the classical formula 
of Landau. If the perturbation of the distribution function has a /134 


-like singularity, then perturbations of the macroscopic quantities are 
not damped with the passage of time. If the perturbation of the distribu- 
tion function has a discontinuity of the n-th derivative, then the perturba- 


tion is damped according to the law t-(ntl) (and not according to the 
exponential law, as in the case of Landau damping). 


Let us investigate the mechanism for establishing fluctuations 
(which are not dependent on the initial perturbations) of the macroscopic 
quantities in a non-equilibrium plasma. We shall show that these fluctua- 
tions are established due to the "survival" of a singular component in 
the expression for the distribution function perturbation. 


The dependence on time of the perturbation of the k-th component of 
the Fourier potential ¢,(t) in an unconfined plasma is determined by the 


following expression (Ref. 1) 


ge (t) = (2ni)t “i paperidp, (1) 
where Ane? N (ks), 
eee E (2) 
D (k, p) = 1 — maT (w) r: 
N (ks p) = § elo) iee as 
w=", F,(w) = §Fo(v)dv., galw) = fge(v)dva, vi =v —k $? 
Fo(v) is the unperturbed distribution function; g,(v) -- the Fourier com- 


ponent of the distribution function perturbation in the case of t = 0; 
integration is performed in formula (1) along the line Re p = o lying to 
the right of all the singularities of the function op: 


Formula (1) makes it possible to determine the behavior of the poten- 
tial ¢,(t) with an increase in t. As is known, the asymptotic behavior 
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of the function (t) for large t is determined by the nature of the 
singularities of the Laplace transform dp: The function bp was 


determined above for only large values of Re p (in this region it has 

no singularities). In order to study its singularities, we must first 
determine this function in the entire complex variable plane -- i.e., 

we must analytically continue the determination of (2) to decreasing 
values of Re p. The analytical continuation of p is determined according 


to the previous formula (2) along the imaginary axis p. 


For purely imaginary values of p, the denominators in the integrals /135 
which determine $,, vanish in the case of w = ip/k. Therefore, for 


analytical continuation of Pax in the region Re p < 0, it is necessary 
to deform the integration path in the integrals (3) and (4), so that it 
passes around the pole w = ip/k from below. Deformation of the path 


assumes, in its turn, the possibility of analytical continuation of the 
functions Fo (w) and 8k (w) determined initially only for real w in the 


region of complex values of w. 


Thus, a clarification of the singularitites for the function cr» 
which determines the nature of the asymptotic behavior of $,(t) for large 
t, requires a knowledge of the analytical properties of the functions 
Fo(w) and g;,(w). 


Let us confine ourselves to investigating the functions Fo(w) per- 
mitting analytical continuation in the region of complex values. The 
function D(k, p), determined in the case of Re p > 0 by relationship (3), 
may be continued analytically in the region Re p < 0, by determining it 
everywhere as 


_ 1 _ trie | Fo (w) dw 
D(k, p)=1 me | (5) 


where integration is performed along the real axis w with passage around 
the pole from below in the case of w = ip/k. 


We have found the denominator of expression (2) for $¢. -- i.e., the 


function D(k, p) -- over the entire plane of the complex variable p. Let 

us now calculate the analytical properties of the numerator for this expres- 
sion, i.e., the function N(k, p). Formula (4) determines it in the case 

of Re p > 0. As has been pointed out, the function N(k, p) has no singu- 
larities in this region. The position and nature of the singularities of 
this function are determined by the properties of the function g,(w) in 


the case of Re px 0. 


If the function g,(w) has singularities (which may be integrated) for 
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real w, then the function N(k, p) will have singularities for purely 
imaginary p. In particular, such a situation is observed if the func- 
tion g,(w) has a -like singularity, a discontinuity, or a break, and 
also if any of its derivatives has a break (in these cases, the func- 
tion g,(w), generally speaking, does not permit analytical continuation 


on the real axis. 


If the function 81, (w) has no singularities on the real axis and per- 


mits analytical continuation in the region of complex values of w, then 
the function N(k, p), and consequently the function bps will have no 


singularities on the imaginary axis p. However, generally speaking, it /136 
may have a singularity in the case of Re p < 0 at the points p = -ikw,, 


where w,. is any singularity of the function 8p (w) lying in the lower 
halfplane of the complex variable w. 


Let us elaborate further on the nature of the asymptotic behavior of 
k(t) for noninteger functions g,(w). In this case, singularities of the 


function N(k, p) are added to the singularity by determining the roots of 


the dispersion equation D(k, p) = 0. The distribution of these singulari- 
ties depends only on the form of the function Bk (W) ~- i.e., on the nature 


of the initial perturbation -- and does not depend on the plasma properties 
(on the function Fg(w).) As has been indicated, one significant property 


of the singularities for the function of N is the fact that they may all 
lie only in the left halfplane p. Therefore, if only one of the roots 
Py = -iwr-Yr of the dispersion equation D(k, p) = O lies in the right 


halfplane p, yp < 0 (which corresponds to the possibility of an oscilla- 
tion increase), then the nature of the initial perturbation has no signifi- 
cant influence on the asymptotic behavior of $ (t) in the case t t >œ, 

If N(k, p) has singularities at the points p = py = -y,-iw, (n = 1, 
2, 3, ...), the contribution made by these singularities to the asymptotic 
behavior of $,(t) in the case of t > œ may be written as D, n exp {—Int — iont}, 

n 

where o, represents certain constants. Adding this sum to the contribu- 
tion from zeros D(k, p), we find the asymptotic expression for $,(t) in 


the general case of noninteger functions g,(w) (which have no singulari- 
ties for real w) 


“ge (£) ~ Vol? exp {—yt — iv,t} + $ an Exp {—Ynt — font}, (6) 


where g(r) is the residue or at zero of the function D(k, p) (the point 
cr 
P = Py = -Yz;~iw,). 
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Thus, for large t the potential k(t) represents superposition of 
the eigen plasma oscillations, whose complex frequecies w;-iy, are deter- 


mined by the plasma properties [the right sum in (6)], and oscillations 
whose complex frequencies w)-ty, are determined by the form of the initial 


perturbation g,(w) [the second sum in (6)]. The eigen oscillations may 


be both damped and intensified. The oscillations whose frequencies are 
determined by the form of the function g,(w) may be only nonincreasing 


(i.e., damped or oscillating oscillations with constant amplitude). 
Let us give two examples of oscillations whose frequency and damping 


decrement are determined by the initial perturbation, and do not depend 
on the plasma properties. 


As the first example, let us investigate oscillations produced in /137 
the case of 


SoH 


BO) = ay ae (7) 


where 9, Wo, w are certain constants. In this case we have 


$ TK 
N (k, p= o Ta ° 


The function N(k, p) has a singularity in the case of p = ~ikwg-kw), 


which introduces the following contribution to the asymptotic behavior 
of $),(t) in the case of t > œ 


or (£) ~ go exp [—kwt — ikw}. (8) 


Thus, the frequency and damping decrement of oscillations produced in 
the case of initial perturbation such as (7) equal kwọ and kwj, respec- 


tively. In the case of wı > 0, the damping disappears. We should point 
out that the function 8p (w) acquires an ô-like singularity on the real 


axis, g,(w) > 7g95(w-wo). 


Let us study oscillations produced in the case of the discontinuous 
functions 8, (w) Let us set, for example, 


Eo (—Wm<w< Ww), 9 
gw) = ( (wl > w). P 


Thus, the function 
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_ Zop PH ikw 


has branch points on the imaginary axis p, p = + ikw. The contribu- 
tion made by the singularities of the function N(k, p) to the perturba- 
tion of the potential (t) has the following form 


ge (£) © gy TELD, (10) 


The 6-like singularity on the real axis on the function g,(w) leads to 
nondamped oscillations of the potential k(t). The discontinuity of the 


function g,(w) -- i.e., the ô-like singularity of its first derivative -- 
leads to potential oscillations which are damped as t=}. It can readily 
be shown that the discontinuity of the n-th derivative of the function 

8, (w) -- i.e., the 6-like singularity of its (n + 1)-th derivative -- 
leads to asymptotic behavior of a t~("+l)exp{ikwgt} like potential, where 
Wo is the discontinuity point. 


Let us determine the manner in which the fluctuations of macroscopic /138 
quantities, which do not depend on the initial conditions, are established 
in a collisionless plasma with an arbitrary (not necessarily equilibrium) 
distribution function (Ref. 2). In order to do this, we should note that 
the averaged product of the distribution function fluctuations for particles 
of the a-th type can be represented as follows at corresponding periods 
of time 


(BEV) BE (V1) ) = dad (k + KNB (v — VY FE (V) + Yaar (VV K) ay 


where the first component (Ref. 3) describes the "correlation of the particle 
with itself", and the second component is related to the interaction between 
particles (and is determined by the previous history of the system). It 

is important that the first component contain 6(v - v'), while the second 
component is a smooth function of velocity. 


In order to obtain the correlation function of the potential, we 
must express the potential ¢,(t) by gf(v) by means of relationships (1) - 


(4), and we must then perform averaging by means of formula (11). It may 
be readily seen that the presence of the 6-like component in formula (11) 
leads to a nondamped (and oscillating according to the law expfikwt}) com- 
ponent in the expression for the potential correlator. The remaining 
(smooth) components in formula (11) are damped according to the law 
exp{-y,t}, where Yk is the customary Landau damping decrement (i.e., the 


imaginary part of the root w, of the dispersion equation e(w,, k) = 0). 
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Thus, in the time ‘eo the potential fluctuation distribution (and 


the distribution of all other macroscopic quantities), which does not 
depend on the initial perturbation, is established after the outer per- 
turbation is shut off in the plasma. 


According to relationships (1) - (4), (11), the correlation function 
of the potential has the following form 


(g(t) ¢u(0)) =8(k Ka) f do] (w,k)[-?exp { — iot) x 
x (12) 


x Me? f dvFo(v)3(w — ky). 


This relationship was obtained by Rostoker (Ref. 4) by employing a 
different method. Our derivation presents a clearer explanation of the 
mechanism, and makes it possible to determine the time required to estab- 
lish fluctuations, which do not depend on the initial perturbation, in a 
nonequilibrium plasma. 
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CHARGED PARTICLE INTERACTION WITH A TURBULENT PLASMA 
I. A. Akhiyezer 


This article computes the energy lost (or acquired) per unit of time 
by a charged particle when moving through a turbulent plasma. The depen- 
dence of the particle energy change on the magnitude and direction of its 
velocity is established in explicit form. It is shown that the turbu- 
lence spectrum elements have no influence on this dependence. 


The strength of charged particle interaction with the plasma is deter- 
mined by the level of plasma fluctuations. In particular, the particle 
energy losses per unit of time P are related to the charge density corre- 
lator <p qu by the following relationship 


where ez, W, v are, respectively, the particle charge, mass, and velocity. 
If the plasma consists of cold ions and hot electrons moving at the mean 
velocity u with respect to yhe ions, the “barge: density correlator in 

the "sound region" (q(Tį/M) 12 << w << q(T,/m) /2, aq << 1) may be repre- 


sented in the following form 


(Pg =I *(ag)*{|T (9, qu) — $o] (o — qs) + 


L (2) 
T [r (9, —qu) — pw] (w + gs)} 
where Te, T is the Pepea and m, M -- the mass, respectively, of 
electrons aad ions; s = (T,/M) Hae speed of sound; T(q, qu) -- effective 
temperature of sound waves; a= T a! 2(4re 2n)=}/2 —— Debye radius. Substi- /140 
tuting (2) in (1) and assuming Ton T >> uv? > we obtain 
—YVs 
P= care (q'dq i dg feos 8 — (3 == 1) sin 8 cos JE T (q, 7), (3) 


where 9 is the angle between v and u 


s? 1/3 r 
n= qu = gu [È cost + (1—5) sin 0 cosg}. 


In the wave vector region in which the sound waves are damped (qu < 
< qs), the effective temperature is (Ref. 1 - 3) T,Q - qu/qs)~*. Close 


to the boundary of the stability region (qu qs) the function T increases 
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sharply. In the case of qu > qs, the linear theory predicts an exponen- 
tial increase in T with time. This increase stops due to nonlinear phe- 
nomena, and a stationary fluctuation distribution is established which 
is characterized by a very high effective temperature, i.e., the state 
of stationary turbulence (Ref. 4). 


It follows from this behavior of the function T(q,n) that 9T(q,n)/an 
has a sharp maximum for a certain value of n close to qs, n = no % qs. 
Noting that the energy losses are determined by the derivative dT/dn, 


and not by the function T itself, we can thus express P by a small number 
of parameters which characterize 9T/dn in the case of n X ng, without 


including more detailed properties of the turbulence spectrun. 


Close to n X no, we have 


fara a =E (RNa +o —2)}, (4) 


where Aj 2 equals unity in order of magnitude, and T* is a large quantity 
equalling T(q, n) in order of magnitude in the case of n > qs. 


Substituting (4) in formula (3) and assuming that 6, < © < 6_, where 


cos 04 = (uv)! {s? + (v? TA s*)'/2(u? —s*)'}, (5) 


we obtain 


27% — s? 
P= ee a — (cos 84 — cos 8)-'/: (cos 8 — cos 8_)-":, (6) 


where a = aJ(aq)3(AjA2)~!dq. This relationship determines the dependence 
of P on v in explicit form. In the case of 0 < 69(@9 = arc cos s*/uv), 
the particle energy decreases and in the case of 6 > 69 it increases. 


In the case of |6 ~ 0,|< T,/T*, relationship (6) ceases to be valid. L141 
In the case of |0 - 6,| << T,/T*, we have 


s 2 1 . 
P=+ ERT Bs CA (u? Ta s?) (v? =s $?)—"% ; 
ap(uvT, sin 65) hs 


(7) 


where a, =z af (ag) Xr Az "dg (the signs «+> correspond to 0 % 0,). 
Thus, P is proportional to T* in the region 0, < 90 < 8_, except for the 
boundary of this region where P ~% (T*) /2, 


In the case of v > u, the critical value of the angle 6, strives to 
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zero, and relationship (7) does not hold. For 6, << T,/1T*, and employing 
formulas (3), (4), we obtain 


(ezT *)? 


P = aol, 


as, aa = a f(ag)? dz? dq. (8) 


In this case, the energy losses are particularly large (they are propor- 
tional to (T*2). 


If 6 < 6; or 6 > 6_, then the expression for the energy loss does 
not contain the large parameter T*/T,. For le - o, | << 1, nevertheless, 
P is proportional to Je - e, |7 /2, and consequently it is large. 


In order to have critical values of the angles 0,, it is necessary 
that both u and v exceed s. If u > s, then 0, > 6, = arc cos s/v. In 
the case of 0 X 9c, the energy losses are propotional to crx) 1/2 

(ezs)*(T*T,)'/s I ee 
P= yO to = gaflar Ay dg. (9) 

When expressions (6) - (9) were derived, it was assumed that the 

difference v - s was not too small. If 1 - s?/v? << Te/T*, the energy 


losses will be at a maximum in the case of vu = s? (in this case P is 
e by formula (8) and sharply decreases with an increase in 
vu - s‘/. 


In conclusion, we would like to point out that the dependence of P 
on the angle @ holds, even if the assumption that 93T/3n is small in the 
region n > no is not fulfilled. The contribution made by the quantity 


dT/dn with n > no in the expression for P can only change the function P 
somewhat in the case of 6, < 6 < 6_, without changing P in the case of 
6 % 6. Consequently, the nature of the dependence of P on the angles 


is not changed in the case of 6; < © < Ə_. In particular, the function 
P, which is positive in the case of 0 = 6; and negative in the case 
of 6 = 8_, must vanish for a certain value of the angle 9 = 6, Oy < 

< 69 < 6_ (thus 69 can differ somewhat from arc cos s*/uv). 
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THEORY OF NONLINEAR MOTIONS OF A NONEQUILIBRIUM PLASMA 
I. A. Akhiyezer 


As is well known, low frequency oscillations with a linear law of 
dispersion -- so-called ion sound ~- are possible in a collisionless 
plasma consisting of hot electrons and cold ions (Ref. 1, 2). It is 
interesting to study nonlinear motions of a collisionless plasma consisting 
of hot electrons and cold ions, and primarily simple waves. The study of 
simple waves not only makes it possible to trace the development of an 
initial perturbation, but it is also of interest as an independent investi- 
gation, since only the region of simple waves can (when there are no dis- 
continuities) be contiguous to an unperturbed plasma [see (Ref. 3)]. 


A. A. Vedenov, Ye. P. Belikhov, and R. Z. Sagdeyev (Ref. 4) have 
studied simple waves in a two-temperature plasma on the basis of an ios- 
thermal hydrodynamic model. This article investigates simple waves in a 
nonequilibrium plasma on the basis of a kinetic equation, without employing 
a special model. A system of equations has been obtained for the moments 
(introduced in a specific way) of the electron distribution function, which 
enabled us to determine the direction of change for quantities characterizing 
the plasma in a sound wave, and to trace the development of a perturbation 
having finite amplitude’. 


` Equations Describing a Simple Wave /143 


The system of equations describing the motion of a collisionless 


* Yu. L. Klimontovich and V. P. Silin (Ref. 5) have investigated the 
problem of a hydrodynamic description of a two-temperature plasma 
without collisions. Nonlinear motions of such a plasma were studied 
in (Ref. 4, 6). 
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plasma consisting of hot electrons and cold ions has the following form 


ð G) e nð z 
(+e +ee2)rao, 


PERTE: = 0; on + div (nu) = (1) 
div E = 4; ne ((Fdv — n); rot E = 0, 


where F(v) is the electron distribution function; n and u are, respectively, 
the ion density and hydrodynamic velocity; E -- electric field; m, M -- 
electron and ion masses, respectively. (It is assumed that the electron 
mean energy considerably exceeds the ion mean energy.) Being interested 

in sound oscillations whose phase velocity is small as compared with the 
mean thermal velocity of electrons, we do not have to take the term dF/dt 
into account in the first of the equations (1). Confining ourselves to 
one-dimensional plasma motions and making allowance for the fact that the 
charge spatial distribution is small in a soundwave 


jn —f Edy | ~n (alh) <n 
(A -- the length at which the quantities characterizing the plasma change 


significantly; a -- Debye radius), we can reduce the system of equations 
(1) to the following form 


OF M du OF . dn duo. 
oe mae T a ta =% (2) 
n = § Fdv, 


where d/dt = 3/dt + ux (the x axis is selected in the direction of wave 
propagation; the subscript x for the velocity component u, is omitted from 
this point on). 


Let us introduce the "moments of the distribution function": 


Dj (x, n= (soy mf s Adv, j=0, 1... j 


Employing system (2), we obtain 


oD; M d o d 
a t+5 a Pits = 0; wt nt — 0; Di =n. (4) 


In order to study nonlinear plasma motions, system (4) is more advan- /144 
tageous than the initial system of equations (2), since it includes terms 
which are only dependent on x and t, while equations (2) also include the 
electron velocity distributions. In this sense, equations (4) are similar 
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to equations of hydrodynamics, although -~ in contrast to hydrodynamics 
which operate with a finite number of quantities -~ they include an 
infinite number of “hydrodynamic quantities" n, Dj, u. 


We are interested in simple waves, i.e., those plasma motions for 
which the perturbations of all quantities characterizing the plasma are 
propagated at the same velocity -- in other words, for which each of the 


functions X [X = u, n, Dj» F(v)] satisfies the following equation 


(at V(x, t) 5) X =0. 


In the case of simple waves, as is well known, all of the quantities X 

can be represented in the form of a function of one of them (for example, 
n), which in its turn is a function of x, t. System (4) thus changes into 
a system of customary differential equations for the functions D; (n), u(n), 


and the phase velocity V(n) is determined from the solvability condition 
of this system. After simple transformations, we obtain 


du Vs, 4aDi Dian, 


antn? Gn ~ D,’ (5) 
mn \'4 
V =u + EVs, v; = (2) J 


where £ = +1 (e = -1), if the wave is propagated in the positive (negative) 
direction of the x axis. 


The determination of the electron distribution function in the case of 
a simple wave may also be reduced to solving the customary differential 
equation. Rewriting the kinetic equation (2) in the following form 


ƏF (v) M (V~u)? Fv) 
me a do TO 


n u,0v, 


and introducing the notation 


F (o3; va 2) =F (v; n) 


[ve = (vy> v,)], we obtain 
F (oss vi n) =F (2 + B(2); w), (6) 
where the function B(n) satisfies equation /145 
ap, 2M Vs _ o 
dn m n | (7) 
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We should point out that Landau damping of sound waves was not 
taken into account when the initial equations (2) were derived. There- 
fore, it is necessary that the wave amplitude An not be too small, 
An/n >> (m/M)°/2, in order that equations (2), and consequently relation- 
ships (4) - (7) may be valid. The role of nonlinear phenomena in the 
development of the perturbation is much greater in fulfilling this condi- 
tion, than is the role of sound damping. 


Development of a Perturbation_Having a Finite Amplitude 


The system of equations for "hydrodynamic" quantities (5), together 
with relationships (6) and (7), enables us to study the direction of the 
change in quantities characterizing the plasma (including the electron 
distribution function) and to trace the development of a perturbation 
having a finite amplitude. 


First of all, let us determine the manner in which the electron dis- 
tribution function changes in a simple wave. It follows from (7) that 8 
decreases in a contraction wave, and increases in a rarefaction wave. 
Therefore, for values of v for which 3 F/V 2 Vx < 0, the number of electrons 


having a velocity ín the (v, v + dv) range increases in the contraction 
wave, and decreases in the rarefaction wave. Conversely, at values of v 
for which 3 F/v,8v, > 0, the number of electrons with velocities in the 


(v, v + dv) range increases in a rarefaction wave, and decreases in a con- 
traction wave. In particular, if the initial electron velocity distribu- 
tion has a spike encompassing a small velocity region, along with a maximum 


for v = 0, the spike shifts to the region of larger (smaller) values of 


RA as the contraction wave (rarefaction) moves. 


Let us dwell in somewhat greater detail on the case of Maxwell distri- 
bution. Employing system (5), we may state that in this case Vg, Dj/n do 


not depend on n, and consequently are motion integrals. The electron tem- 
perature and the distribution function F/n, which is normalized to one 
particle, do not change during wave propagation. Thus, in the case of 

a Maxwell velocity distribution of electrons it is valid to describe a 
two-temperature plasma by means of isothermic hydrodynamics. 


In order to determine the manner in which the form of the sound wave /146 
changes, it is necessary to compute the derivative dV/dn [see (Ref. 7)]. 
Employing the system of equations (5), and assuming, for purposes of 
definition, that € = 1, we obtain 


av V nD 
T= sa ce) 
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Depending on the electron distribution function, dV/dn may be positive, 
negative, equal to zero, or alternating [positive for single values of 
the parameter $ and negative for other values of this parameter, see 
formula (6)]. (We should point out that the derivative dV/dn is always 
positive both in customary and in magnetic hydrodynamics.) 


If dV/dn > 0 for all values of the parameter 8, then (just as in 
customary hydrodynamics) points with a large density move at a large velo- 
city. Therefore, discontinuities arise at the contraction sections.** 
Self-similar waves are rarefaction waves. In particular, this possibility 
exists for a Maxwell electron velocity distribution and for a distribution 
in the form of a step F v © (vj(n) - v7), © (x) = 4a + sign x). 


If for all values of 8, dV/dn = 0, all the points move at the same 
velocity during wave propagation. Therefore, the wave profile is not 
deformed and no discontinuities arise. Employing equations (5) and (7), 
we may state that the velocity of two-temperature sound and the guantity 
pg /2 change in an inverse proportion to density, Vn = const, Bn“ = const. 


The case dV/dn = 0 is realized, in particular, for a Cauchy distribution 
Fa {vő (n) + v2}, 


If dV/dn < 0 (independently of the value for the parameter 8), then 
points with a large density move at a low velocity. Therefore, discon- 
tinuities arise in the rarefaction sections. Self-similar waves are con- 
tractions waves. This possibility is realized, in particular, if the dis- 
tribution function is the superposition of two Cauchy distributions 


F ~v, (n) (0? + of (n)}7? va (n) fo? + 03 (n)}-?. 


We should note that in this case the velocity V, increases in the rarefac- 
tion wave, and decreases in the contraction wave. 


Finally, let us discuss the case when dV/dn may be both positive and 
negative, depending on the value of the parameter B. For purposes of 
definition, we shall assume that dV/dn > 0 in the case of 8 > 8; and 
dv/dn < 0 in the case of B < B1, where 8; is a certain critical value of 
the parameter 8. When a contraction wave moves in such a plasma, the 


žk We have employed the term "discontinuity" to designate the narrow 
regions in which the gradients of the quantities characterizing the 
plasma become so large that the initial equations (2) are not 
applicable. In the case of a/A >1, sound dispersion must be taken 
into account. With a further increase in the gradients, multi-flux 
flows [see (Ref. 4)] or shockwaves may arise in these regions. 
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"apex" of the wave (points with the density n > n}, where nı is deter- 
mined from the equation ß(nı) = 81) lags behind the "base" of the wave 
(points with n < nı). Therefore, the density at the point of the dis- 
continuity which develops from such a wave cannot exceed n}. A discon- 
tinuity may also arise during the motion of a rarefaction wave. Thus, 
the density cannot be less than n, at the discontinuity point. If, on 
the other hand, dV/dn < 0 in the case of 8 > 82 and dV/dn > O in the 
case of B < B2, then -- as may be readily confirmed -- the density at 
the discontinuity point cannot exceed ny when a discontinuity develops 
from a rarefaction wave, and cannot be less than nz when a discontinuity 
develops from a contraction wave (the critical density ng is determined 
from the equation B(n2) = Bo). 


Both of the above possibilities may be realized, in particular, if 
the electron velocity distribution is a superposition of two Maxwell dis- 
tributions, "hot" and "cold", 


F(v; n) = v, (n) exp { — 51: | + v2 (a) exp { -mh TiD Ts. 


For small density values (n < (a,Tj*/a:TY*)"”"), according to equations (6), 

(7), Y= 417) ve = GgnT/T (a,,a,  -- constants). It may be readily confirmed 

that in this case dV/dn > 0 in the case of n <n, and dV/dn < 0 in the case 
2a, TYAS 


of n > nj, where n= (2 7) . For large density values (2 > aT Y*lasTY*) 
2 1° 


vy, = ayn! yo = aon (ai, a2 -- constants). In this case, dV/dn < 0 for 


n < ng and dV/dn > O in the case of n > n2, where nọ = 2 (T Ta)". 
aa 


Let us investigate the motion of a two-temperature plasma arising 
during its uniform contraction or expansion [similarly to the problem of 
the plunger in hydrodynamics, see the monograph (Ref. 7)]. We shall 
assume that the plasma occupies the halfspace x > Vot, which is uniformly /148 


limited by a moving plane. (Such a boundary may represent, in particular, 
the region of a very strong magnetic field.) As is well known, oniy self- 
similar waves (in the absence of shock waves) can be steady motions of 

a uniformly contracted (or expanding) medium. If dV/dn > 0, a self- 
similar wave (which is in this case a rarefaction wave) arises during 
plasma expansion (Vp < 0). If dV/dn < 0, a self-similar wave (which is 

in this case a contraction wave) arises during plasma contraction (Vp > 0). 


Employing formulas (5), (6), (7), we may relate the change in all 


the quantities X characterizing the plasma in a self-similar wave with the 
"plunger" velocity Vo. Assuming, for purposes of simplicity, that 
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Vo << V „> we obtain 


V; Dı 7’ 
AV, = z (1— $) Vo; (9) 
J 


AF (v) = — vv, E0, 


0 > 
v,0v, 


AX = X k=vut pE X|e- coe 
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NONLINEAR PROCESSES IN A UNIFORM AND ONE-COMPONENT PLASMA 
N. A. Khizhnyak, A. M. Korsunskiy 


Nonlinear solutions of a one-dimensional kinetic equation without a 
collision term, which depend on x and t by the combination § = x = Vot, 
where -- Vo is the constant wave velocity, were compiled in (Ref. 1) 149 
and studied for several cases in (Ref. 1 - 4). The conditions at which 
these solutions may be realized were found in (Ref. 4), and a limiting 


transition to small oscillations was performed. 


General Theory of Nonlinear Waves 


This article investigates the more general nonlinear solutions of 
a one-dimensional kinetic equation without collisions 


BD ay EO he (1) 


where ọọ is the potential of the self-consistent electric field with the 


factor ~, determined by the Poisson equation 


ap (2) 


dne? 
Je 7 Z£ (f fdu — n), 
where ngo is the unperturbed density of ions whose mass is assumed to be 


infinitely large. It is assumed that the distribution function f depends 
only on the variables u and ¢. Then the electron density 


(a) = f fdu =F(9), (3) 
the density of the electron flux 
(nu) =JSufdu = ® (ẹ) (4) 
and the energy density 
(nF) = fp E du=Y (o) (5) 


are explicit functions of only the potential ġ. 


We shall show that in this case the main plasma characteristics may 
be compiled within the framework of a hydrodynamic approximation, and 
that the electron velocity distribution function may be found relatively 
simply. From the equation of continuity 
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2 are : ar ae a nu) =0 


we find that the potential ġġ must satisfy the following equation 


C) 
where Vo = 2 is a certain velocity of longitudinal waves in the plasma 


which depends on the potential ¢. It can be readily seen that the specific /150 
solutions depending on x - Vot, where Vg = const, are special cases of 


equation (6) in the case of Vo(¢) = const. 


From this point on, it will be assumed that the function Vọo(¢) is 


given, and it may be employed to formulate the solutions both for the 
equations of the hydrodynamic approximation and for the kinetic equation. 


The equations of the hydrodynamic approximation 
ð ð 
E E 1e (7) 


may be reduced to equation (6) and to the equation for hydrodynamic velo- 
city v(x, t). In actuality, since 


( nu } 
(na)? 


it follows from expressions (3) and (4) that v = v($¢). Therefore, the 
Navier-Stokes equation (Ref. 7) has the following form 


y= 


APRE a = 1 dp æ 


and is identical to equation (6), if only v(¢) is determined by the equa- 
tion 


’ 1 , 
v-o t1+7 Py 


U 


= V, (9). (8) 
? 


The electron plasma density <n> may be found from the specific wave 
velocity Vo(¢) and the hydrodynamic velocity v(¢) from equation (8), by 
means of the following relationship 

e 


a CSAS de 
u (9) — Vo (P) 
F(9) =ne % i , 
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which follows from the equation of continuity. Thus, the dependence 
of the main hydrodynamic plasma parameters on the potential ¢ is always 
formulated in quadratures for a certain Vo($). 


The potential of the self-consistent field $ is determined according 
to equations (6) and (2). The general solution is found in the following 
form from equation (6) according to a certain function Vo ($) 


9 = (Cc), 


where c = c(x, t) is the equation of its characteristic. The characteris- 
tics of a quasilinear equation in partial derivatives (6) may be com- 
piled according to the well known method (Ref. 5). 


We shall regard t, x and ọọ as functions of a certain parameter s. /151 
Then the parametric equation of characteristics can be written as 
follows 


d . d 
d l; z = V, (P); z 70 


from which it follows that ¢ equals $9 and does not depend on s. There- 
fore, x = Vo(¢)s + xgand t = s + tg,where xg, to and ¢9 are the values of 
t, x and ọọ on a certain line s = 0 through which the characteristics pass. 
In particular, if primary interest is directed toward the development of 
moving waves, which is caused by the nonlinearity of the medium, xo, to 
and $9 must characterize a given initial moving wave. 


Let us assume that at the initial stage of the process there is a 
given moving wave with a constant phase velocity wo. We then have 


t= 1; x= Wot, Po = Po (1), 


where $9(t) characterizes the initial dependence of the field potential 
on time. In this case, we have 


t=s +s; x = Va (Po (*)) SHWT P = P(t) 
and after excluding s we find t from the following equation 
x = Vo (g C) (E — 1) + ws, 
which determines the characteristics of the quasilinear equation (6) 
T = T(X, t). Using the characteristics from equation (2), let us deter- 


mine the field potential in the plasma ¢(x, t) and all of the hydrodynamic 
parameters of the mediwm. 
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We can find the electron velocity distribution function from the 
kinetic equation (1). Actually, sinee equation (1) is transformed into 
the following form, together with equation (6), 


lu— Vo (oN 5-H =O 


the equation of characteristics 
dp 
da =U + Vo (9) 


enables us to find the particle velocity distribution function in a 
general form from the potential ọ and the kinetic velocity u. Since the 
potential ¢ may be compiled within the framework of the simpler hydrody- 
namic approximation, the specific dependence of the distribution function 
on the coordinates and time has been established. 


Propagation of Waves in Media with Linear Dispersion /152 


We shall assume that the given law for the dependence V)(¢) of the 


nonlinear wave phase velocity on potential determines the law of the 
medium dispersion. We shall study nonlinear, nonstationary waves in a 
medium with linear dispersion 


Vo (y) = Vy + 1? 


where Vg and y are certain constant parameters. 


We shall disregard the plasma pressure, so that equation (8) assumes 
the following form 


d i 
a =— oFV, +49. 


The dependence of the hydrodynamic velocity v of the plasma on the poten- 
tial can be determined according to the following relationship 


9 = ce” + eH, 
where x = v - Vo; c -- integration constant. [In the case of v = Vo,. dis- 
persion disappears, i.e., in the sense that Vo is the limiting plasma 


velocity at which the electric field potential vanishes. We then have 


p= h [0X + 1) — er]. (9) 
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Since 
a ae 
dp v — Vo (P)? 


the electron plasma density may be found according to formula (9) 


x dọ 


? = 
4 dọ dx 
f pov, AF J —rer dX 
n = Nne., = N£ i 
Since C 


f z f 

dx aa | diy) e™ 
-m dX = = In 
[x — yel k 1 — e” 1—nt ’ 


%e Xe 


e™ 
F(?) = "07 


where xg is the value v = Vp at which the electron density equals the ion 
density. 


Consequently, the Poisson equation can be written as follows 


and x may be found from relationship (9). 


Let us determine the characteristic of equation (6) for the linear 
dispersion law 


x = [Vo + T90 (*)] (£ — 1) + wot. 


Let us assume that at an initial moment of time the field is switched on 
whose potential increases linearly with time: 


T 
Po (1) = Po y 


We then have T ; 
t+ (F Vo—m) —#}e +55 © — Vol) =0, 


from which it follows that 
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alri (w — o 


1 tiy Frie ge T ww, — vali T w Vet). 


Consequently, the characteristic is real only up to a certain value Xpounqg 
which depends on time 


x =V + alee E (w — vo]. 


bound 


If X>X ound» the signal which is switched on in the case of t = 0 has still 
not reached the point under consideration. If X<Xþound» the characteristic 


is determined in more than one way (it has two values). The requisite 
value of the characteristic is determined by additional considerations. 


The propagation rate of the front boundary is 


dX bound _ 
Vpound™ dt Yot te 4 9 Ply 


Consequently, for y > 0 the front velocity increases with the time, while 
in the case of y < 0 it decreases. 


The velocity of a point with a constant given potential, as may be 
seen from (9), is determined by the wave phase velocity Vo($), i.e., by 
the specific value of the potential ọ. We may formulate the potential and 
the electron velocity distribution function by the hydrodynamic parameters 
which are found. 


Electron Velocity Distribution Function in the 
Case of a Square Dispersion Law 


Let us find the electron velocity distribution function in the case of 
media with a square dispersion law 


Va (p) = Vo £V —2¢. 


The equation of characteristics 


d — 
a E — X + V —29, 
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where xı = u - Vo; u -- electron kinetic velocity in the plasma, can be 
integrated by substituting ¢ = xit. 


We finally find 
Tuo a 
| var (10) 
x Se eer ONE ek Set ey ~ 


for a positive sign of the root, and 


ni| -2 y —29 1 o 
RT |u V U-V} 
29 =e 
2y/ 2 _ le 
|- V =e (11) 
x 


V5+1 


2W wate 
faa 


VY 5-1 


= Cz = const 


for a negative sign of the root. 


The solutions of (10) and (11) enable us to compile the distribution /155 
functions in the case of ¢ > 0 which change into a Maxwell distribution. 
On the other hand, all the distributions changing into Maxwell distributions 
in the case of 9 > 0, in regions with a non-zero electric field, have the 
following form 


f =AexpC, for u> Vw 
f= AexpC: for u <V- 


Consequently, electrons have velocities lying outside of the boundaries 


2Y —2% 
Vs—1 


2V —29 
V5+1 


u>V,+ 


and 
u < V — 
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The electrons whose velocities are included within 


2V —29 2 V —29 
V, — S= V -=-—, 
are damped by a wave, and are not included in the distribution functions 
of (10) and (11). The relationships obtained enable us to study the 


collisionless transition of electrons from the region of trapped particles 
into other plasma electrons and the associated energy distribution. 
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INDUCED SCATTERING OF LANGMUIR OSCILLATIONS IN A PLASMA 
LOCATED IN A STRONG MAGNETIC FIELD 


V. D. Shapiro, V. I. Shevchenko 


This article investigates the nonlinear interaction of harmonics in 
the long wave spectral region of Langmuir oscillations (kyp,<<1). The 
linear damping of these oscillations, which is caused by the interaction 
with resonance particles, is negligibly small. It is assumed that the 
plasma is located in a rather strong magnetic field, so that the plasma 
particle oscillations are possible only in the direction of the magnetic 


field which is parallel to oz. 


The initial system of equations for the distribution functions of 
electrons and ions and the electric field has the following form 


of 3 oft. 
art ilka — wa) + aire ap Bede i, 
kz "2 a 7 k —q2 “~z (1) 
T een 
RÈ. = ta Yie, {fait m 
> ikr N 
= sy ay re. ft) 
E=7hE_¢ k Heres pag i + c.c. (3) 
k k 


(fọ is the background distribution function whose change with time can be 
disregarded, due to the small number of resonance particles (i — kipe < 1) 


and the small dissipation of oscillation energy during scattering). The 

notation in equations (1) - (3) is standard; summation in equation (2) is 
performed for plasma ions and electrons. We obtained the following rela- 
tionship for damping frequency and decrement, disregarding the nonlinear 

terms in the first equation, from equations (1) and (2): 


— 3 42 A x be OF. 
o = vod + zk a3.) cos 8; T> =7' ae cos 0 ay. (4) 


(0 is the angle between the direction of oscillation propagation and the /157 

magnetic field). The nonlinear interaction of harmonics leads to a change 
in the oscillation spectrum due to processes of wave decay and scattering 

by Plasma particles. The laws of conservation must be fulfilled in the 


* coe = = complex ‘conjugate. 
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case of two-plasma decays 


Assuming that ky and k, lie in one plane, we obtain the following condi- 
tion for the spectrum (4) from (5): 
k, cos 6, + ka cos 8; 


= cos 6, + cos bz. (6) 


Decay is possible if kı; > 0, koz < 0, which corresponds to the sign "--" 


in condition (6). In several cases, the spectrum of Langmuir oscillations 
in a strong magnetic field is a nondecay spectrum, particularly if it is 
close to a one-dimensional spectrum 0, ™ 65%. 


Nonlinear wave scattering is caused by the interaction of plasma parti- 
cles with the beats of different frequency. This process becomes signifi-~ 
o,—o 
cant if the condition —* "=< Ue. is fulfilled. The beats which are 
l2 ~~ *22 
the cause of wave scattering cannot arise due to decay, since the following 
condition is fulfilled for the waves formed during the decay 


w, — t, o_, 

pH =~ => oe 

k =h > > K 
1z 2z z P = ka| 


Therefore, in this case the processes of wave decay and scattering are in- 
dependent. The transformation of the oscillation spectrum due to the non- 
linear scattering process will be subsequently investigated. 


In solving the nonlinear equation (1), we shall employ the method of /158_ 
perturbations. Substituting £$, which is found from the linear theory, 
k 


in the nonlinear terms of this equation, we may employ (2) to obtain the 
following formulas for the electric field amplitude and the distribution 


function in the second approximation 


* Ion oscillation branches are not examined. Actually, in the approxima- 
tion under consideration, decay of a Langmuir wave into a Langmuir 
ion-sound wave is possible. However, the energy primarily remains in 
high-frequency oscillations (Ref. 1). 
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Bee Be Hay ato. — os) 


cual 3 a hid 
go zí y= i DY Rog a A kaa a Redy 
k m? R? — 
2 a E e ee +o, *) 
q k—q q 
fo 
1 af 9%, = 
s kv, —o, e dv, \ 9,%,— 4, se 
k—q q q 
` E, — E, exp —i Or + @_, — oO, 4 
E ines > E eg al S : 
k a? Rv, eg a ay 
> k—q 4 
4 (7) 
afo 
A R E pe E E E. . 
m? 0v, \ 9,0, — 4 Mm, k? ov, = 
£ E > + O_» ‘) 
k—q q 
sne f 1 ð du, do. 
m kv, — 9 —w, Ov, \ qV — 


In the third approximation, the system of equations for determining 
the electric field amplitude and the distribution function has the 


following form 
of : i a Ca 2t ea >? 1 3 > 
e ERAS E ie e AS (e9 E + 
q 
(2) 1) ; 
HEE Se (On ee sla: j 


RE = 4r a 5 v, 
ikE~ aie f F du 


2 2 

where rl), g2 are determined from formulas (7). We obtain the kinetic /159 
k k 

equation for waves from these equations by simple computations: 


3E 


k 
-a = 1-8 


ys H (È, q, *) E „E, xX (9) 
k kz e k—qz q—xz2z 


qa 


+ + o aed i 


POR bar er a 
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where Y is the decrement of Landau linear damping: 


te, +to,-e, 


> = k a 
A(k,9,%) =— 5: Ant hE 
€ ee +o, „Ho; k 
k—q q-% x ) 
4net 1 
a ( » af | ras ug oe So 
a k—q q—xz x 
oo 
1 ð ðv, > 
x duz f 9222 —%, .—®, dv,\%,0,—_, du + 
q-* x z 
qz l drez 1 
tg me m? kpr — 2a oe, o=o, a: 
ae a> ! ons #) b =f q-* + 
q—« x a 


ð 1 dfo 
vz 92%, — §_, e (kz — q2) v — o, _ } Ov, 
q x k—q 


of 


4nes 1 dv = 
oD es eee ee a 
m 9292 —- 9, 1 Oo, Ov, %z0z — O_, 
B q x i 


Multiplying (9) by E* , averaging over time (only the terms with 
kz 
> > > > 
= k and x = q - k, which change slowly with time, remain in the right- 
hand side) and combining the equation obtained with the complex conjugate 


/160 


one, we can write 


(11) 
ar |g f= 25 [E gf + Peele lal 
ie = Zim ( wii tn Gd j= 
q 
2 es a 
k ðv 
=F age im ( ai (te — 9) | aaae * 
{en ( i) 
k 
1 — 4)? 


Miao an i(k 95 won do +l 
re) ee a rae 
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du] ). (12) 


In computing the integrals included in equation (12), we assume that 
the following conditions are fulfilled: w, >> kvp, which corresponds to 
k 


weak wave absorption, and w, w, which is necessary in order that a sig- 
k q 

nificant number of plasma particles may participate in the interaction 
with beats of different frequency. Let us investigate the case when the 
contribution made by ions to the plasma polarizability at a different 

oi 
(w =, o)? 

m k q 
pore Or is fulfilled. Then, assuming that the distribution func- 161 
i" “De 


frequency is negligibly small -- i.e., the condition 


tion of the plasma particles £3 is a Maxwell distribution, we obtain the 


following from equation (12) after very cumbersome computations 


k — Èn 
220e ETIA *% 
(13) 


2 (cos 8 — cos 8") m ge (cos 8 — cos 6’)? 
x len ee exp See oom en’ 
3 (k? — q?) X? De cos QT, (k — q)? cos? 8 


27 T e 
u 6’) = - E 0k?q? 
q" (k, 9, q, 0) 2 z gaa COS% 


The increment FAK, 0,4,0") which determines the induced scattering rate 
of plasma oscillations in k-space differs considerably from zero only in 


v 
the narrow range of angles: ev — 815 5" in the case of 6 #0 and 


p-e y E in the case of 6 %0. 


Thus, during scattering the spectrum which is initially one-dimensional 
remains close to a one-dimensional spectrum. In the case of 6" - are, it 


coincides with the increment obtained previously in a one-dimensional model, 
within an accuracy of the factor cos 6 (Ref. 2, 3). It increases somewhat 


2 1 o7 ; 
R up to Thha z" TEA (0 =0), and 


with an increase in 0 - 0' = 
y 3ye 


then rapidly decreases to 0. 
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In the real case of a confined plasma the spectrum of the values 0 


v 
is always discreet, and for a sufficiently small ratio =, only one 
© 


possible value 6' = 6 can enter the angular range 6' ~- 6 in which bh 
is large. The scattering of plasma oscillations in a strong magnetic 
field then takes place, in fact, in the same way as in a one-dimensional 
model, in contrast to scattering of oscillations in a plasma which is not 
located in a magnetic field. In the latter case, scattering at large 
angles is possible, and, if the angle between the wave vectors of two 


V 
interacting waves considerably exceeds Te, y*2increases by a factor of 
Vo 
TAR as compared with yn in a one-dimensional model (Ref. 4, 5). The 
e 
contribution made by ions to the plasma dielectric constant is quite sig- 
2 cos? 0 
nificant when the condition A z>1 is fulfilled. We should note 
Gea) 
that this condition is fulfilled most readily in the case of 0 % 6, when 
the difference “z - w, is at a minimum: Op = 0 = 2 oe cos 0 (k? — q°)à 2, and 162 
q 

w2, cos? 8 m 4 l m, i m, 

To,—o.\? m, pat” Assuming that the conditions 7° me! H X 

Car Tae aes 

xaz «1, are fulfilled, we obtain the following expression from (12) in 

De 
the case of 6 = 6' for yet 
4 7 e [MA 23 g?> cos 8 
Utp 0g. 09) — 4 Le BON at oe, SCORE , (14) 
T (k0, a, 27 2 meo, lk—al ktoe 


i.e., in this case yx" increases somewhat as compared with (13). In the 
PE 
“or Me 


v. e 
case of 00 (ee 6 — cos 6’ ~ ale i ~X _1 and when the condi- 
(a A 

1 


m, 
my RPA? 
De 


€1 is fulfilled, yt is determined in this region by formula 


We should recall that this investigation pertains to a case when the 
plasma is located in a magnetic field which is so strong that the plasma 
particle oscillations are only possible in the direction of the magnetic 
field. The condition of "magnetization" of the electron component, as is 
customary, has the form woe << wy,, i.e., it is fulfilled for field 


strengths which are not too large. In this case, when the ions make a 
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significant contribution to the dielectric constant of the plasma, it 
must be required that the ion plasma component be "magnetized". The 
corresponding condition is harder: 

my 


hd2,,. 


|; T =| < end es He D Moe = 
k 4 


e 


Let us turn to certain general characteristics of the nonlinear 
change in the wave spectrum. Since the increment of nonlinear scattering 


Yk, 854,09") in the one-dimensional case 6 = 6" changes sign when k? q 


is substituted [see formulas (13), (14)], the change in the total oscilla- 
tion energy during scattering in the approximation under consideration 


eqnele-2ete d NJE Pe y yh, 09,0) /E.P [Ef =0 
Tor ca 
eal 


= (15) 


Employing formulas (13), (14), we may also readily see that the non- /163 
linear interaction of harmonics leads to a transfer of energy along the 
spectrum to smaller k: 

a 2 Lik —a)yi lk, 9, 9, 0) |En] En 2 <0. 
rr a rd ea a ms e a 
k ka 


(16) 


The total oscillation energy during scattering by plasma particles 
changes in the subsequent series with respect to Kea » Since in this 
e 
approximation an addition to y? appears, which is symmetrical with respect 


> 

to the keq substitution. Assuming, for purposes of simplicity, that 
the oscillation spectrum is one-dimensional, we obtain the equation for 
the change in the total energy in the spectrum: 


KAJ 


! 3 
a ee = Ee ee ar i — 
ade; |= D 5 ae, d 
k k 


RY FB Dee aT 


4 2° mê 


(17) 


> 


kg 


The second term in this equation which describes the oscillation energy 
change during scattering becomes more significant than the first tern, 
which characterizes the oscillation energy change as a result of their 
interaction with resonance particles, when the following condition is 
fulfilled 
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BIE, P 3 zr 
1 73e 2k 
I E 2 D 18 
4xN,1 ~ 60 € a (15) 
e 


i.e., for relatively small oscillation amplitudes [E, |? << NoT, if the 


VTe 
parameter Fs is fairly small. 


The possible dissipation of oscillation energy when they undergo 
nonlinear scattering by plasma particles was pointed out in (Ref. 6, 7). 


However, for Langmuir oscillations this phenomenon is Kn? times less 


than the change in the field intensity in separate harmonics in the spec- 
trum due to energy transfer. 
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NONLINEAR THEORY OF LOW FREQUENCY OSCILLATIONS EXCITED 
BY AN ION BUNDLE IN A PLASMA 


D. G. Lominadze, V. I. Shevchenko 


As is well known, in the case of the interaction between a bundle 
of rapid electrons and a plasma, the energy lost by the bundle during 
relaxation changes into thermal energy of electrons in the plasma and in 
the bundle, and into energy of high frequency Langmuir oscillations 
(Ref. 1). When investigating the possibility of heating the plasma ion 
component during bunched instabilities, it is of interest to investigate 
the excitation of low frequency oscillations by the bundle, in the non- 
linear approximation. A previous article by D. G. Lominadze (in collabora- 
tion with K. N. Stepanov) investigated the linear theory of low frequency 
oscillation excitation in a plasma located in a magnetic field by an ion 
bundle. In a strongly non-isothermic plasma (T, >> Ti)» during the passage 


of an ion bundle, longitudinal longwave k Yia << ') oscillations may be 


Ha 


excited, whose frequency is determined by the following relationship 


1 1 2 292 20” 
og (0 + ef) ty (8 ea) aes si (1) 
where ; wh; g Amn, yp _ Te, 
o= 1? "0 T TM’? “De T inng’ 
L+ pag, 
Hie is the Larmor frequency of ions and electrons, respectively; 0 -- /165 


the angle between Hy and the direction of oscillation propagation. If 
T; ~ Tg, the ion bundle can excite longitudinal shortwave oscillations 


LY Ti 
—— ~ 1} in the plasma with the frequencies (Ref. 2) 


Hi 
Ia (e 
Wn = Noy (1 + Pa) Pa= aet ; , 
Tì m — WY nH) 
ee ee (2) 
NAA, Bees 
2 
kivey 
where wy = 3 T, (44) is the Bessel function of the imaginary argument. 
w 
Hi 


The process by which bunched instabilities develop may be 
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qualitatively divided into two stages (Ref. 1). In the initial (hydro- 
kv, 
z 


dynamic) stage, the bundle remains monoenergetic << 1, vp, — 


Ky 


k 
thermal velocity of bundle ions, and instability develops very rapidly; 


Ys} 
s~ (2) "(io n, -- densities of plasma and bundle, n, << ng, w -- 
2 


excitable frequency). For rather large oscillation amplitudes, the 


k, u 
thermal energy in the bundle is so large ( J ~1} that the bundle 
T 


relaxation may be investigated in the quasilinear approximation (quasi- 
linear stage). The time required for the development of this instability 


i he ord OL 
Tquasi ÍS on the order of a a 


This article investigates the development of instability at the first 


stage, the most unstable oscillation branches are found which produce the 


dynamics of the instability development, and the change in the macroscopic 


parameters of the bundle and the plasma is determined (thermal energy, di- 


rected velocity) in the case of instability. In a strong magnetic field 
(wyi >> Wg) it is possible to trace the development of instability at the 


quasilinear stage and to determine the state at which the bundle and the 
plasma arrive as a result of the quasilinear relaxation process. 


Excitation of Longwave, Low Frequency Oscillations 
(Hydrodynamic Stage) 


We shall assume that at this stage the following condition is ful- 
filled 


Rv.) <i; kdu l 
Mae an (3) 


k k 


(ŝu is the change in the bundle velocity at the initial stage), at which 
the dispersion equation of longwave oscillations has the form 


(4) 


When investigating the dispersion equation (4), we can examine two cases. 
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l1. wy << üs (weak magnetic field). During Cherenkov excitation 
(w = kz up), oscillations with the frequency 


2 
ow, (: pian ) r 


2 
os 


have the largest increasing amplitude increment. The maximum value of 
the corresponding increment is 


3 ‘V/s i 
ngg rs ° 


The formation of instabilities in this case is possible in the case of 
v 
cos 8 < a, if Uy ? Vg» and for any cos 0, if Uy < Vs (7s ya} 


2. Hi >> We (strong magnetic field). In this case, ion-sound waves 
with the frequency w = W, cos 6 have the largest increment. The increasing 


increment of these waves is 


vs y 
e = Ta F Ww, COS 6. (6) 
v 
Excitation of oscillations is possible if — > 1. 
ug 
Relationships (4) - (6) were obtained by disregarding damping by /167 


plasma electrons, which is valid if the following conditions are fulfilled 


ez» 


no Uy 
If condition (3) is fulfilled, all of the bundle particles are in resonance 


with the wave, and its hydrodynamic description is possible by means of 
the moments of the velocity distribution function 


u" = ~ ç ofsdo; Ty, = +f (v: — ut) (0 — ug) Fedo. 


We can obtain the equations describing the change in these quantities with 
time from the kinetic equation of the Fokker Planck type, which is derived 
in (Ref. 3) 

af, a (o: n) 

Of v, \ $k Dog) * (7) 
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B = e, i pertain, respectively, to plasma electrons and ions and bundle 


ions; a;;, -- diffusion coefficients in velocity space: 
5 g ba) a2) 
ai= DJE Ai RE 2o [La + Thi; 
‘ k? kyo, © nd? (Ny) (8) 
aiz = S21 = s; maler | z onp 2 eon —T,): 
alla EJO Ty + Y I3 (2) E-a + Tal. 


n=l 
Here we have 


ku 
k aes, 
r 


= oe + mon) +e ae 2 oy 
k 


r 
The frequency wr and the increment 1 ate determined by the disper- /168 


sion relationship of the linear theory (4). 


We obtain the system of equations for the change in the "bundle" 
parameters due to the development of instability from formula (7) 


7 To 1- 
(fauo — a ui) +1, (Fao — v + emi) + © : 
aT yy 2 We jk. Wied (9) 
fe > ai (too PEE? 
dT 
a= Lie; 
k 
CNN i 
(o — of on)? F $ (Fzto— o CE ent) + qe 
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The change with time of jz |? is described by the following equation 
k 


ajep (10) 


mae La 


We may obtain the equation for the change in the plasma ion para- 
meters from formulas (9), assuming that ub = 0. The condition for the 


<< 1 is 


kAv. 
applicability of a hydrodynamic description of plasma ions i 


less hard in the case of ny << ng than equation (3), and this description 
is applicable throughout the entire development of instability. 


Integrating equations (9) with respect to t, we obtain the equations 
determining the change in the energy of directed bundle motion and in the 
thermal energy of ions of the bundle and the plasma. 


In the case of a weak magnetic field (Ons << wy), we have /169 
n,Mu,su = = $ 3DE; K 
k 


poor S BATEZ 


2 u 


we fs 
nT,ı ~ = (a) (i -4i alezi <a> 


no (11) 
1 aT e RELA; rf 
2 0 ae 3 u 
è aaz EET 
nol u= g ~ 5 ae 8x | zl 
In the opposite case (wg >> Wa), we have 
1 1 
n,Mu,du = — aaa T D| E> 2. 
o — 
E k 
AS | 
nati = (2) " SB D] Ef" (12) 
—- k 


nT u =. z >| E> 2, 


The change in the electron distribution function in the case of 
instability is determined according to the following equation 


qt e 
of ð |e ke k alo 
2-36 . oar | EJE; "Paka . (13) 


The phase velocities of the excited oscillations change between vp; << 

<< vg << Vre: Therefore, the main portion of plasma electrons is in 170 
resonance with the waves v >> Vg. The change in the electron distribu- 

tion function for these v is described by the equation 


of o kz k f 
Sf aloe: Sa BIG Biot 30, ° (14) 


Thus, for a change in the electron thermal energy in the case of insta- 
bility, we obtain the following in the case of a weak magnetic field 


Eng ist. m= Def (15) 


and in the case of a strong magnetic field 


FT = —25- a Des | (16) 


For plasma electrons which are in resonance with excited oscillations, 
whose velocities lie in the narrow range ax £ Vir: we obtain the following 
Te 


expression from equation (13) 


2 
f JE, et (e; — ho) 5. (17) 


a Atel 


T, 
Changing to the variablet=2(1/ 7 <t< 1) in formula (17) and 
s e 


representing fje in the following form 
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he = ete (1- AC ): (18) 


2 x 
where ọ © Dheo = 5, p(s.) is the portion of the electron distribution 


function depending on v,, we obtain the following equation for the change 
in $(&, t) 


2 
ap 1 e % 2 sin8 ig L. el, 
eee aia mite 3) (as) 
The interaction of resonance electrons with oscillations leads to /171 


the occurrence of a plateau in the electron distribution function in the 
region of excited oscillation phase velocities. The time required to 
establish the plateau may be determined from formula (19): 


= on” m2" 55 cos 4 


_ (48. 1 æ STE Ept go (20) 
D 0 


(A& ~- the dimensionless width of the plateau). 


The change in the energy of resonance electrons is 


na) . m7. | Ez |< x E. K (21) 
k 


aa) = 

M, \m k 

i.e., it is considerably less than the change in the total energy of plasma 
electrons. 


We may determine the balance of energy during the development of 
instability from formulas (11) and (15), (12) and (16): The energy of 
directed motion, which is lost by a bundle, changes into the thermal 
energy of particles in the bundle and the plasma and into energy of 
electrostatic oscillations. Since conditions (3) must be fulfilled in 
the initial (hydrodynamic) stage, we may readily estimate the maximum 
energy of oscillations excited at this stage: 


DO “mes (22) 
k 


At this oscillation energy, the time required to establish a plateau in 
the electron distribution function, as follows from formula (20), is 
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1 
a?’ (23) 


i.e., it is a little less than the time required for the development of in- 
stability in the initial, quasilinear stages. 


Thus, the plasma electron distribution function changes most rapidly 
in the resonance region: A plateau appears in this region in the velocity 
range in which oscillations are excited at a given moment of time. A 
change in the spectral density of oscillations and ion parameters occurs 
much more slowly. Under these conditions, the electrons have no signifi- 
cant influence on the dynamics of the instability development. 


Excitation of Low Frequency Oscillations /172 
(Quasilinear Stage) 


Further development of instability leads to a still greater increase 
in the thermal scatter in the bundle, and its distribution function be- 
comes so diffused that the quasilinear approximation is applicable. 


In the case of w, >> wyy, it follows from the expressions for the 
diffusion coefficients (8) that a,j; % a,,, if kôv ~ w. The problem is 


thus essentially a three-dimensional problem. Let us investigate the 
opposite case, Wy; >> We, since in this case a,j; << azz at the quasi- 


linear stage, i.e., only longitudinal diffusion is significant. 


The initial system of equations for the quasilinear approximation 
has the following form 


2 


og__ 9 ok Og 
F ERr in Esa de — kava) 5 ; (24) 
k 
OJE_ |? 
| z] = zo," , 2898 E}? (25) 
ot eM k*dv ft ra 
-E 
"=i 


Here g(t, v) is the distribution function of bundle ions integrated with 
respect to vj. 


Substituting = E from equation (25) in formula (24), changing 
k 
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> 
from summation over k to integration, and integrating with respect to k 
due to the 6-function, we obtain 


E .}2 
Og to nm a ali. or Ve 
Ot Det ne” MA oS f sina, Leh . Sa 
La 
#(1—3) 


Integrating with respect to v and t, we may write the following 
equation 


C sin 049 (1E (9) 2 2M m, 1 
f {Ez | -|E Op = 2a B.S x 


(27) 


U 


5 u? Ys o 
xo ( 2a) f lgt, v) — g0, o'av’, 


where vı is the lower boundary of the instability region caused by damping /173 
by plasma ions. We may employ this equation to determine the oscillation 
energy at the quasilinear stage: 


Og 


Yad Epp = ae (o(t— 3) term. v)— es 
k 7 A 


0: 


— g (0, v')]dv’dv, 


where g(%, v) is the distribution function at the end of the quasilinear 
stage -- the plateau, whose height is determined from the following condi- 
tion 

ni 


glo, v)(vs— 04) =n; g (oo, 0) = 5: (29) 


vı and vy (the upper boundary of the instability region) are determined by 
the following relationships 


&(, vı) = gi(%); g(co, və) = go (V2); (30) 
go(v) -- the distribution function of the ion bundle at the beginning of 
the quasilinear stage; g; (v) -- the distribution function of plasma ions 


integrated with respect to vj. We obtain the following from (30) 


oT. Se Ys 
y= y n” NED X lo; ve = uo(1 aa] ). (31) 
ny 


Performing integration with respect to v in (28), we may write 
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2. 
w= g Ming TE “ay (32) 


The change in the energy of the plasma ion thermal motion may be deter- 
mined from (11): 


1 2 
ne 1 kz 1 1 
= w2,— DE pk cet 
2 NÒT ri "2% 2| Ez | Das 6 Mno}. (33) 
k k 
The increase in the energy of the plasma electron thermal motion is 
1 l l 1 vf 
FNT pe => VIE.P = =a Mn, 2 34 
g otl ge ar 2| ra AB, T 10 ig (34) 
k 


The finite bundle velocity established at the end of the quasilinear stage /174 


1s 
“Us 


u” = \ vg(œ, v)dv = 2 =F, 
U: 
i.e., the energy loss of the ordered bundle motion 
be = — 3nMuj. (35) 


The thermal bundle energy acquired during the development of instability 
is 


U: 


nòT 2 Mu2 
E n= fe. yo 25%) zier (36) 


v: 


We may employ formulas (32) - (36) to verify the fact that the law of con- 
servation of energy is fulfilled 


PSE eyes 37 
Bet |E ty nT ie + peel e+ ET 1 = 0. (37) 


k 


Thus, in contrast to the excitation of high frequency oscillations, 
the excitation of low frequency oscillations leads to the transfer of a 


considerable portion ve of the bundle energy to plasma ions, and atso 


leads to significant heating of the ion component. 


The quasilinear theory disregards the nonlinear phenomena of oscilla- 
tion scattering by plasma particles. As was shown in (Ref. 1), in the case 
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of the excitation of Langmuir oscillations by an electron bundle, non- 
linear interaction of harmonics is usually insignificant in the quasi- 
linear stage. However, there is a parameter region in which these phenomena 
determined the dynamics of the instability development (Ref. 4). 


We shall continue to study the role of the nonlinear interaction of 
harmonics in the case (which we are considering) of excitation of ion- 
sound oscillations by an ion bundle. 


Excitation of Shortwave, Low Frequency Oscillations /175 
(Hydrodynamic Stage) 


When allowance is made for the finite, Larmor radius of plasma ions, 
excitation of longitudinal oscillations by the harmonics Miss is possible. 


One important feature of these oscillations is excitation in the case of 
T; ~ Te and propagation almost perpendicularly to the magnetic field. 


Therefore, it is more likely that low frequency oscillations are the 
reason for anomalous plasma diffusion perpendicularly to the magnetic 


k 
field, rather than ion-sound oscillations, for which £ << 1. 


Zz 


Let us investigate the manner in which the ion bundle and the plasma 
parameters change at the initial excitation stage of low frequency oscilla- 
tions. 


The dispersion equation for shortwave oscillations for conditions (3) 
has the following form 


2 2 
Tio a of n (Hi) 201 29 or ly 
1+ T~ i 2; acron PST AT a (© — ko)? — o fi (38) 
n=— o 
x sin? @ = 0. 


Waves with frequencies of wp = nwy;y (1 + Yn) have the largest increment. 
Their increasing increment in the case of n = 1 is 


where 
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The change in the ordered velocity and thermal energy of bundle ions 
may be described by means of (9), where w, and y, are determined from 
k k 
formulas (38) - (39). Making use of the fact that y„ has a maximum with 
k 
respect to k, we may obtain the following in the case of instability at 
the frequency w = wy; (1 + y1) 


2 
` £ 1 
mMudu=— | R=) ge ad 7 | ef 
R alk) j, 
M k=ko 
AA (40) 
Tu=. a l 
Mèi a No ory, k? TT "aÍ dk) EP; 
2 Fp? 
1 oo; | #249 fn NY te l p 
Shel ip G SPRES 2 
2 on pu Mu ae en au J lager dk 
k=k, 
(Ko -- the value of k at which Yy has a maximum). Let us employ the {176 


k 
approximate formulas given in (Ref. 2), and we shall assume that a(k) 


2 2 
reaches a maximum in the case of Paor 1.5: 


Hi 
0,22 k 
a(k) = mra poo ae (41) 
i 10 0 
I+ 7) 


Assuming that a << py, << 1, we obtain the following expression for the 
0 


plasma ion diffusion coefficients in the case of instability at the fre- 
quency w = wg; (1 + y1) 
ku 
wi L +) 
gis k i nt 


kr (e = er)” FE’ 
k 


k 


1 (42) 
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We obtain the following by means of (42) from formula (7): /177, 
2 2 
1 wi [Rz 1 1 1 aa 
DNT = > >: = oa \ dei E}: 
z "o HC E) Ta aap § | zl 


(43) 


2 a 
oy 1 1 l aA 
nT ,,=(—p- —~ Bi” F | dk] EP 
m ole) jy 
We may determine the change in the plasma electron energy by formula (14) 


1 2 

1 06 l 1 ay 

gT ie “le i Gaps) dk Ef i 
Le 


Comparing expressions (40), (43) and (44), we can see that the 
energy of the bundle ordered motion changes primarily into energy of the 
transverse thermal motion of plasma ions, and may lead to significant 
ion diffusion perpendicularly to the magnetic field. Relationships (40), 
(43) and (44) are valid as long as the conditions of a monoenergetic bundle 
are fulfilled (3). 


The maximum energy of low frequency fields, obtained at this stage, 


is 
a | E-)~ = (i Ty (E) “natu, ass 
k 
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NONLINEAR PHENOMENA IN A PLASMA WAVE GUIDE 
(ION CYCLOTRON RESONANCE’ AT A DIFFERENCE 
FREQUENCY) 

B. I. Ivanov 


Ion cyclotron resonance (ICR) has been extensively studied both /178 


theoretically and experimentally [see, for example, the summary in 

(Ref. 1)]. Several works have appeared recently (Ref. 2 - 4) which 
examined the problems of the nonlinear theory of ion and electron cyclo- 
tron resonance. As is known, the non-linearity criterion has the 
following form (Ref. 5, 6) 


eEgh Uo —! 
2runc? Py (1 S) =l 


v 
(Eg -- strength of the wave field; à ~- wave length; Bg ~~ retardation; 


vo -- ordered plasma velocity). Formation of nonlinear phenomena is 
facilitated during resonance (Ref. 6), since in this case the non-linearity 
parameter contains the additional factor ~ w (w-w,)~°. Thus, it is 

possible that nonlinear phenomena may occur in the case of ICR, because 

for this case the occurrence of large strengths of the wave field, small 
phase velocities, and large wavelengths is characteristic. In this case, 
nonlinear phenomena may play a significant role during heating (nonlinear 
damping, nonlinear shift of the resonance frequency) and during the intro- 
duction of high frequency energy into the plasma (interaction of frequencies). 


In principle, it is possible to introduce large UHF power into the 
plasma at two frequencies, and then to perform ICR at the difference fre- 
quency. Such a mechanism is also possible during the excitation of low 
frequency oscillations in the plasma-bundle system (Ref. 7). As is known, 
in unstable plasmas, low frequency oscillations, whose origin is 
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sometimes difficult to establish, occur simultaneously with high fre- 
quency oscillations. The occurrence of the high frequency oscillations 
is satisfactorily explained by the theory of plasma-bundle interaction. 


One of the mechanisms leading to the formation of low frequency 
oscillations may be the nonlinear interaction of high frequency oscilla- 
tions [separation of difference frequencies, decay instabilities (Ref. 8, 
9)] with the subsequent transfer of energy from high frequency oscilla- 
tions to low frequency oscillations [parametric amplification (Ref. 10)]. 
This article makes an attempt to provide a model for this mechanism by /179 
which low frequency oscillations are excited close to the ion cyclotron 
frequency. The parameters of the apparatus are the same as in the pre- 
ceding article (Ref. 12), which investigated the nonlinear distortions 
of the signal form and the formation of combined frequencies. We employed 
generators having a small power (1 w) which had a relatively small dis- 
turbing influence on the plasma which was produced independently. In 
order to fulfill the non-linearity condition, it was necessary to operate 
at low frequencies (f ~% 1 Mc) and with low phase velocities (Bẹ ^% 1072), 
which, in its turn, made it necessary to employ a low-density plasma 
(n ~ 10? cm-3) (Ref. 12). On the other hand, fL X fei >> Vin (Yin © 108 p 


-- collision frequency) represents the necessary condition for observing 
the ICR at the difference frequency. In view of these considerations, 


the main ("beat") frequencies and the difference frequency were of 
order of magnitude one: fı % f2 v fL% 1 Mc. 


In order to observe the weak ICR signal, a sensitive system of a 
balanced, high frequency bridge was employed (measures were taken to re- 
duce the noise level). Figure 1 shows the diagram of the apparatus (1 -- /180 
current regulator; 2 -- voltage regulator; 3 - high frequency generators; 

4 -- amplifiers; 5 -- phase inverter; 6 -- phase rotators; 7 -- AVC unit; 

8 -- two-ray oscillograph; 9 -- heterodyne receiver; 10 -- heterodyne fre- 
quency meter; 11 -~- self-excited oscillator; 12 -- main anode; 13 -- quartz 
tube; 14 -- water cone; 15 -- auxilliary anode; 16 -- cathode; 17 -- mag- 
netic field recorder; 18 -- palladium filters). The plasma wave guide 
consisted of the following, which were distributed coaxially: A plasma 
core with a diameter of 1 cm, a lead tube with a diameter of 3 cm, a 

copper casing with a diameter of 23 cm with the cross section along the 
generatrix, and a solenoid. The total length was about 180 cm. The 

quartz discharge tube was evacuated from both sides to a vacuum of 

~ 107" n/n2, after which hydrogen was introduced from both sides through 
the palladium filters. When the entire length of the tube was continuously 
evacuated, it was possible to obtain a constant pressure. Before the 
measurements, the tube was treated to preliminary processing with pro- 
longed, high frequency discharge (wavelength } % 2 m, generator power 

P ~ 500 w). The plasma was produced by discharge at a constant current, 
and the discharge current was stabilized. In order to increase the 
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Figure 1 


ənization coefficient, a cathode was employed made of lanthanum hexa- 
oride with indirect direct current heating. Two anodes were em- 

Loyed in order to obtain a stable discharge: The main anode and the 
uxilliary anode. The anode potentials were selected according to the 
lasma noise minimum. 


The main frequencies fı and fo from the generators were excited in 
ne plasma wave guide by short spirals located close to the left end of 
he wave guide. An adiabatic, absorbing water charge (length of about 
0 cm) was located at the right end of the wave guide. The reflection 
coefficient from the right end of the wave guide k equalled 0.1 - 0.3 
Ref. 13). Thus, moving waves with main (f) and fọ) and combined 
nf; + mf2; n and m —- whole numbers) frequencies could be propagated 
a the wave guide. The output signal was employed on two spirals, one 
f which was located in the magnetic field section which could be 
odulated by a commercial frequency. From the receiving spirals, the 
ignal was supplied to the two arms of the high frequency bridge. Each 
rm consisted of an amplifier, a phase rotator, and the AVC unit. The 
atter was used to eliminate relatively slow unbalancing of the bridge 
Taye 7? Tod? The signals from the arm of the bridge were supplied 


ut of phase to the heterodyne receiver adjusted to the resonance fre- 
uency, and after rectification through the low-frequency filter they 
ere supplied to the oscillograph. In the normal position, both arms 
ere almost completely balanced, but the amplitude of the signal with a 
odulating frequency remained somewhat larger. During the modulation 
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Figure 2 


of the magnetic field, at the moment the resonance value was crossed 
(£_ ot na) -- twice during the modulation period -- a signal for the 


bridge unbalance was produced, which could be observed on the oscillo- 
graph. This recording system has good sensitivity and noiseproof qualities 
(subtraction of the signals from the two receiving channels leads to an in- 
crease in the modulation depth of the carrier frequency by the ICR signal, 
and simultaneously eliminates the plasma noise correlations along the wave 
guide length). 


Figure 2 presents oscillograms showing the dependence of the reso- 
nance position on the magnetic field strength (fı = 1.5 Mc, fo = 2.5 Mc, 
f_ = 1.0 Mc; n = 1-109 em-9, p = 7-1072 n/m’, Hy = 7 ka/m, H % 56-62 ka/m). 
The upper line corresponds to the resonance signal which is inverted 
during rectification. The lower line corresponds to the signal coming 
from the generator recorder of the magnetic field. With an increase in 
the constant magnetic field strength (for a fixed difference frequency 
and a constant amplitude of the variable magnetic field), the resonances 
converge, since the resonance condition (f_ % jci) is fulfilled in the 
negative halfperiod of the variable magnetic field. 


Figure 3 shows the dependence of the resonance position on plasma 
density (fı = 1.5 Mc, fg = 2.5 Mc, f| = 1.0 Mc, H = 62 ka/m, Hy = 7 ka/m, 
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Figure 3 


p= 710-2 n/m’, n & (1 - 2) 10? cm-3). With an increase in the plasma 
density, the resonances diverge, i.e., they shift into the region of 
large magnetic field strengths. 


In order to obtain quantitative estimates, the reflection coeffi- 
cients from the wave guide ends, the magnetic field strength, the phase 
velocity, and the plasma density were measured. 


Dynamic measurements of the magnetic field strength were performed 
by the generator recorder (Ref. 11) (see Figure 1). The self-excited 
oscillator circuit was located in the modulated section of the solenoid. 
Carbonyl iron was used as the induction core. Due to the small dimen- 
sions and the small value of yp (~ 10), the recorder disturbed the magnetic 
field to an insignificant extent. With a change in the magnetic field /183 


strength, due to the dependence (H) the circuit inductance and the self- 
excited oscillator frequency changed. The latter was measured by the 
heterodyne frequency meter according to the zero beats, which could be 
recorded simultaneously with the ICR signal by the two-ray oscillograph 
(see Figure 2). This system was calibrated initially by nuclear magnetic 
resonance, 


The phase velocity was measured by the system shown in Figure 4 
(1 -- AVC unit; 2 -- phase meter; 3 ~- amplifier; 4 -- delay line; 
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5 -~ limiting attenuator; 6 -- generator; 7 -- selective microvoltmeter). 
A wave with the frequency f = f_ X foi was excited in the plasma wave 
guide by the generator. A high frequencg signal was applied to the two 
antennae(a non-mobile short spiral and a mobile whip), and was supplied 
to the phase meter input. The compensation method was employed to 
measure 8g. A coaxial delay line, consisting of segments of the high 
frequency cable (length, 2.5; 5; 10 and 20 m), was switched into the 
circuit of the mobile collapsible-whip antenna. These high frequency 
cable segments could be subsequently combined in any combinations. With 
a minimum distance between the antennae, the delay line was completely 
introduced, and the phase meter indicator pointed to zero. As the dis- 
tance between the antennae increased, the delay line decreased so as to 
compensate for the phase shift produced. In spite of the fact that the 
phase meter system was designed so that the phase reading was not depen- 
dent on the signal amplitude, there was a possibility of error for small 
phase shifts and significant changes in the signal amplitude. In order 
to eliminate this possibility, an AVC unit was switched into the circuit 
of both antennae, and also a selective microvoltmeter (for controlling 
the high frequency signal amplitude) and a limiting attenuator (for main- 
taining the amplitude at a definite level) were introduced into the cir- 
cuit of the mobile whip. 


The length of the delay line AL changed linearly as a function of 
the distance between the antennae £. The retardation was determined 
according to the following relationship 

l 


= — ~ 3. —2 
Po TA 3-10 
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(e -- dielectric constant of cable insulation). 


The electron density n was measured by the shift in the eigen 
resonator frequency Af, and the dependence Af(n) was determined by pre- 
liminary calibration with respect to the electron bundle (Ref. 12). The 
relative content of atomic and molecular hydrogen ions in the plasma was 
not measured. 


The following conclusions may be drawn on the basis of the measure- 
ments performed. At the moment that resonance is passed, the amplitude 
of the "difference" wave increased. The resonance frequency increased 
with an increase in the magnetic field strength, and decreased with an 
increase in the plasma density, while f_ ~ fei’ The picture observed 
corresponds qualitatively to the excitation of ion cyclotron waves at 
the difference frequency. 


The experimental data agree qualitatively with the dispersion rela- 
tionship for ion cyclotron waves 


®, 
ot = (iga). 
a 


ci 


The quantitative divergences (stronger experimental dependence of resonance 
frequency on plasma density) do not as yet yield to a satisfactory explana- 
tion. 
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SECTION IV 


EXCITATION OF PLASMA OSCILLATIONS 


RADIATION OF ELECTRONS IN THE PLASMA-MAGNETIC FIELD /186 
BOUNDARY LAYER 


V. V. Dolgopolov, V. I. Pakhomov, K. N. Stepanov 


The cyclotron radiation of electrons in the plasma-magnetic field 
boundary layer can make a significant contribution to the energy balance 
of thermonuclear reactors with a small density, which employ magnetic 
grids to contain the plasma. This problem was examined in (Ref. 1). 


The thickness of the transitional layer between the plasma and the 
magnetic field may comprise several Larmor electron radii 


Ve T eBo 
e= igo 7 Vz 8 = ine} 


The trajectory curvature of electrons having a velocity on the order of 
Ve is also on the order of pẹ» Moving along such a trajectory in a 


vacuum, a non-relativistic electron radiates the following energy per /187 
unit of time 


2,22 
d € OBY, 
r E, T&me. (1) 


The number of emissive electrons per unit of layer area is nope (ny -- 


plasma density). Therefore, the total intensity of cyclotron radiation 
from unit of layer surface -- if it is assumed that all the electrons in 


the layer radiate the same way as in a vacuum —- is 
dw eno ne a 
I ‘as at NoPe aru 3 . 
Since 
2 
Bo ~ als (3) 
eng *T? 
a cim): ’ (4) 


where a ~ 1. Expression (4) coincides, within an accuracy of a coefficient 
on the order of unity, with the result derived by Burhardt (Ref. 1). 
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We may employ expression (1) only if the refractive index n for 
radiated frequencies w % ùp is close to unity. Im the case under con- 
sideration, the radiated frequencies lie in the region of the anomalous 
skin-effect 


. Q c 
Ren~Ima~ zn! (5) 


(2 = V = 
~m  —- plasma frequency}. In a medium with a large, complex 


refractive index (5) the radiation intensity differs greatly from the 
radiation intensity (1). For example, in a dense, non-relativistic 


plasma (2> "t ab) » the intensity of cyclotron radiation of the main fre- 


quency decreases, as compared with radiation in a vacuum, by a factor of 
aT (Ref. 2 - 4). 


Let us determine the intensity of cyclotron radiation, assuming that 
the radiated waves correspond to stable plasma oscillations. In the case 
of radiation equilspr tus the radiation flux falling on the plasma-magnetic 
wT 
AR 
the plasma I(w)and the flux reflected from the plasma Ipj;R (R -- reflec- 
tion coefficient). We thus find that I(w)= Ipj;(1 - R). In the case under /188 _ 


v 
consideration, L- R% Ea < The width of the radiation spectrum, 


field boundary, IkJ = » equals the sum of the flux emanating from 


caused by the Doppler effect during radiation and by the nonuniformity of 
the magnetic field in the boundary layer, is on the order of wp, in the 
case under consideration. Therefore, the total intensity of electron 
cyclotron radiation in the layer is 


ope eT Po end'T? (6) 
4nics anh? 


I~ 1 (o) ~ 


A comparison of (6) and (4) shows that allowance for plasma polariza- 
me 


tion decreases the intensity of cyclotron radiation by a factor of T° 


Due to plasma resonance, radiation of boundary layer electrons strongly 
increases in the region of frequencies w which are less than the maximum 
Langmuir frequency 29, but considerably greater than wp. In the resonance 


region where e) = ff ey, the intensity of electron braking 


radiation greatly increases at the frequency w. 


Let us investigate absorption of waves whose electric vector lies in 
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the plane of incidence (the XY plane), and who impact on the plasma 
layer from a vacuum. The electric field component which is parallel 


to the plasma boundary, Ey w eikyy-iut | satisfies the following equa- 
tion 
aE R22 de dE ot 
a . 2 
axt ene Git (F243) & =0, (7) 
we (es) 


where e= 1S, 7,28) = (*) is the dielectric constant of the plasma; 
& 


» (x) = ees -~ frequency of collisions; A -- Coulomb logarithm (it 
x 


is assumed that v << w). Since Q % up, the wavelength i = £ 
e 


= is on the 
W 


order of the layer thickness. 


At the point where Re e = 0 (region of plasma resonance), the wave 


electric field increases sharply (Ey ~% E, Ine v Eg In Si Ex von ~ Eo, 189 


Eg -- amplitude of incident wave). This leads to the fact that a consider- 
able portion of the incident wave energy is absorbed in the layer having 
the thickness Ax ^% as in the vicinity of the resonance point. Therefore, 


the intensity of braking radiation at the frequencies w < 2 is 


end'T 8 
I~ Ips Q,~ ae (8) 

‘ F oe ee me, 
Consequently, the intensity of cyclotron radiation is ~y times less than 


the braking radiation intensity. 


2 
We may employ equation (7) only in the case of v > vg = B /3 Re 


If v < vo, we must take into account the formation of plasma waves in the 


resonance region. 
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RADIATION OF LOW FREQUENCY WAVES BY IONS AND ELECTRONS 
OF A NON-ISOTHERMIC MAGNETOACTIVE PLASMA 


V. I. Pakhomov 


As is known, the propagation of three normal waves -- Alfvén, rapid 
and slow magnetosound waves -- is possible in a strongly non-isothermic 
(Te << Ty) magnetoactive plasma, in the low frequency region (w < wyi). 


Each of these waves may be excited due to ion cyclotron radiation or due 
to Cherenkov radiation of ions and electrons moving along a spiral in 
this plasma. 


This article determines the expressions for the intensities of /190 


radiation of these wave types by ions and electrons. The emissive and 

absorbant capacity of the plasma are determined in the frequency region 

w £ Wyy- The case of a low-pressure plasma is investigated in detail, 
H 


when the Alfvén velocity va = 4mnoM 


considerably exceeds the speed of 


i T 
sound in the plasma v, 4/-—= 
M 
slow magnetosound wave is considerably greater than the refractive indi- 
ces of the two other waves. Therefore, it is natural to expect a sharp 
increase in the radiation intensity of a slow magnetosound wave. It is 
shown that in a low-pressure plasma the intensity of cyclotron radiation 
of a slow magnetosound wave by an ion for the s-th harmonic is a 

s 
times greater than the radiation intensity of Alfvén waves and rapid 
magnetosound waves. 


- In this case, the refractive index of a 


Cherenkov radiation of low frequency waves by electrons of a non- 
isothermic plasma may make a basic contribution to the over-all plasma 
radiation. It is shown that the ratio of the intensities of ion cyclotron 
radiation and Cherenkov electron radiation in the case of w ¥ Sugi is on 
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v 


2s—1/y 
the order of (=) (=), where vp; and vp, are the mean thermal 


vb % 
velocities, respectively, of ions and electrons, and V is the phase 
velocity. * 


Propagation of Electromagnetic Waves in a 


Non-Isothermic, Magnetoactive Plasma 


The general expressions which are given in the appendix in (Ref. 2) 
may be employed to derive the dielectric constant tensor of a non-iso- 
thermic plasma located in the outer magnetic field. The sum with respect 
to particle types -- i.e., electrons and ions -- must be taken. Let us 
assume that the following condition is fulfilled 


2 ve 


k a 
Pa = ni Kl, as=i, & 
ofja 
_ elon, 
Zu = yrk, aa 
o — loye ; 
= — 1, f= +1, #2,...5 
tle Peo” 


ee I 
ae Gua 


These conditions are fulfilled if it is assumed that the plasma is greatly /191_ 
non-isothermic (Te >> Ty). As a result, for the tensor of the plasma di- 
electric constant we obtain the following expression 

ej (k, ©) = aby + chihi; + dejjphp + 

-+ e (xi [xA]; — x [xA] + F A [XA], 


where 


H. =Ê.: 
h=7> ERS 


n? 
a= oy ee pene a (20) —1 |; 
; 1 (1) 
d = — ita [ru + 9 Coe) te? |sgn o; 
0q(Zoe) . = — monB, 
e= TT a iV x ——, Mu cos ô’ 


* A portion of the results have been published in (Ref. 1). 
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where 


6 is the angle between the direction of wave propagation k and the 
direction of the outer magnetic field H. If the frequency w is close 
to the harmonics of the ion gyrofrequency swg; (s = 2, 3, ...), the 


following expressions must be added to (1) 
a = — č = — id’ = 2is,, 


where 


ro Vz _ vs? (sBnsin 0)?s—2 os 
s 2 2's! Bn | cos 6 | i 


(2) 


w — So jy v. 


Z. = . p= T yea 
s VZkuycosð © c =] Mea’ 


The dispersion equation determining the longitudinal refractive 
index ny = n cos 6 as the function of the transverse refractive index 


n; = n sin @ and the frequency w has the following form in the given 
case [see also (Ref. 3)] 


nô +n} (=s tage) ae ka | (3) 
= gay (nh + 2nt)] + a[n? (n? + nå) —n2n3] = 
where 
pee z 1l 2 5 
A=V/ FE e 
nyn n 2n? n? n? n? +o 2 
ta e Any AH Nh 4s nio fns (Mi 4.0) + 
2|(14-Vu)* +2 2 nf, 
we ee z) A | (ns + na) — 
nn Vu) nè e z 7 
ea e A; ae NAS GT u 
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When there is no damping (A = 0), equation (3) has three solutions 
ny = nyj (j = 1, 2, 3). Assuming that 


My =A, tiny, NK 
we obtain the following expression for Th; 


7 A (n,);) Bln 5 (4) 
nii = ay HIN Hi]. 
Qn y (047 — 245) (0 — 1) 


The damping coefficient of the j-th normal wave is 
O; > 
% = Š n 1i COS X, (5) 


where x is the angle between the directions of the outer magnetic field 
and the energy flux of the corresponding wave. 


The dispersion equation (3) may be solved in two limiting cases: If 
the Alfvén velocity is considerably greater than the speed of sound in 
the plasma (v, >> v,) and if the frequency w is considerably less than 193 


the ion gyrofrequency (w << w;). 


Let us examine the case when Va >? Vee The dispersion equation (3) 


assumes the following form 
_9 2 
nf ni (e pte el) oh [ol G+ ed) - 


— 2n; ril $ nz =~ (ni — na) =0. 


(6) 


We may find one of the solutions for equation (6) by assuming that ny ^% 
nj “ng >> nas As a result, we obtain (Ref. 4) 


n? 7 
ni =m+ os. 


The imaginary part of the refractive index of a slow magnetosound wave (7) 


ny = aE ate mut), (8) 


The first component in the right part of equation (8) takes into account 
Cherenkov absorption of a slow magnetosound wave in an electron gas. The 
second component takes into account cyclotron absorption in an ion gas. 


We may find two other solutions for equation (6), corresponding to 
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Alfvén waves and rapid magnetosound waves, assuming that ny% nj % 
~% na << ny. The radiation of waves corresponding to these two solutions 
was studied in (Ref. 5, 6). 


In another limiting case, when the frequency w is considerably less 
than the ion gyrofrequency wy; (region of strictly magnetohydrodynamic 


waves), the dispersion equation (3) assumes the following form 


(n3 —n¥) [n5 — 0? (a? + nà — n?) — nån? + (9) 
+ 2 (n4 —n',)] sis, 
where = 
A’ = Vt : sae LG — n3) [ni (ny — ni +n) + 
i s 
n? (10) 
+ 2ni n? nA] + “4 [203 n? (n? — nå — n? ) — nt (On — n? + nå). 
The solution of equation (9) 194 
ni =n +0(7) (11) 


corresponds to an Alfvén wave. The imaginary part of the refractive index 
of an Alfvén wave is 


2n® (n3 — n?) T n$ (n$, + n4) + 20n% (n — n4) = 


, amn — ont (n? +a)? (12) 
ny 1= 8M ° SR? nS PUTE EAE A 
i's 


Two other solutions (Ref. 7) of the dispersion equation (9) 


nies => [ny +n? —n tV (inat n) + Anan’ | (13) 


determine the refractive indices of rapid waves and slow magnetosound waves. 
The imaginary parts of the refractive indices (13) determining the damping 
of magnetosound waves (Ref. 7) are 


— 4,2 2 2 42 2.2 
if am 3 (nhes— Ma +m) +204 o anani 
12,3 = 8M g n n? (nf San? 

s” g 2,3 (P02.3 n? 3,2) 


f (14) 
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General Expression for Ion Radiation Intensity in a 
Non-Isothermic Plasma 


Let us employ the general expressions (Ref. 4, 6) for determining 
the strengths of the electric and magnetic fields produced by an ion 
moving along a spiral in a strongly non-isothermic magnetoactive plasma. 
As a result, the components of the electromagnetic field produced by the 
ion in the wave guide zone assume the following form 


E}? = Pj;sin YEL 4 


E? = Picos VEY; ot 
D ry 
HÈ = + Picos UH; 
HË. = —njcos( xF 9) EQ | 
where 
2eu.; n n, ,e77*siR 
Pi = aR a SE TP ae 
sin X| cos x — H o (aii — nie) (hr nie) 
ni 
Wy =(+ k. sinX + kyjcosy) R — wst F sy FESH 
4 dn 
+ + sgn (cos X —4!"}; 
dn‘, 
EO = 0, (n yr — ls (a i 
gs = Oy (Mey, + My j€33 — €11&33) Js — kyr. (ni AYN — 
0 
— €33) ieiz — Uy Ej1lE23 + Ug Nija (ieiz H leog ctg JA 
(16) 


EQ = + cosX [ 0, [n2 (iee — iens ctg 0) — erresa]; + 
Svi 
kiro 


2,2 2 2 2 
2 | (nini —njess — M1r + e1833 + €23) + Uini xX 
2 ; g 
x (nj — s1) + v aenzens | | — Sin x fo. [ni js (iei — 
r i , su, 2 
— ieza ctg 0) + eqylees] Js + bor. (ninina — niina + 12823) — 
Lo 


2,9 2 2 
—%y (nini į — nje — 11 Fry + £12823 ctg JA ; 


HÌ = — v, COS Î (fe12€33 + €112€23 tg 0) J; = 


sU 

2 2 ' 20: 

— f ; (ness — £1133 — £23 + €11E23 tg 8) cos § + Vie} sin s] Tex 
1 0 


In expressions (16), vy and vj, are, respectively, the ion velocity 
components which are perpendicular and parallel to the direction of the 
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Fe 
outer magnetic field; rg = = -- Larmor radius of an ion; the argument of the 


L 
Bessel function J, and their derivatives J', equals kjrg. The radiated 


frequencies are 


og = Son: + Ry (ws) 04- (17) 


In the case of s = 0 (Cherenkov radiation), equation (17) assumes the 
following form 


Barale = 1. (18) 


In expressions (16) nj and k; are the positive solutions of the PERR {196 _ 
tion for the saddle point 


dnra) _ denka) 
deo sae ee ee (19) 


The negative solutions correspond to saddle points which do not lie on 

the integration path, and therefore cannot make a contribution to the 
integral. The upper sign in equation (19), as well as in (16), corres- 
ponds to a cylindrically diverging wave; the lower sign corresponds to 

a cylindrically converging wave. When there are several positive solu- - 
tions for equation (19), it is necessary to take the sum . of all the solu- 
tions. Each solution for equation (19) determines the relationship between 
the angles x and 6, i.e., the wave toward which the phase velocity is 
directed at the angle 0, and the group velocity at the angle x to the 
direction of the outer magnetic field. 


We may find the intensity of cyclotron (s #0) and Cherenkov radia- 
tion of an ion per unit of solid angle w, j PY employing the following 
expression 


wy == ERAS — EDHO, 


where the bar designates averaging over time. As a result, we obtain 


2 
Lw —2z R 
U 3] 
Wj = = pert . ari tj sf? s. (20) 
? 


oe sin X pula, 
cos SH: ee ee a v? (n? 5-74)" (4; — 2° ,)? 


where 


Us = EQH + cos (x F 6) (FOV. 


189 


In the case of s = 2, 3, ..., formula (20) is valid for the cases 
of slow (vq Vri) and rapid (vj >> Vri) ions. In order to determine the 


radiation intensity at the main harmonics (s = 1), expression (20) may 
be only employed in the case of rapid ions. In the case of s = 1, the 
results given in (Ref. 4) must be employed to determine the radiation 
intensity of slow ions. For Cherenkov radiation (s = 0), formula (20) 
is also valid only in the case of rapid ions. Cherenkov radiation of 
slow ions is greatly absorbed. 


Expression (20) for ions with a velocity on the order of the mean 197 
thermal velocity of plasma ions has the following form (in the case of 
s=2, 3, ...) 


eog? Diaa 
sil yen, eD 


2 ig \Qs-2, og 
where 9 ee ees E I E Ty Se 
2” (61)? sin x H o? (a2 —n2 Dia 


dn | 
cosx -2 
X [[F n} sin (x F (an, — ey; sin 8 + iere sin 8 — iez cos 8) + (e11 — 


L 
— ier?) (€33 cos X F issa sin X) + e2, cos X — €334 cos X ] cos 8 [(e;; —. 


— Fe 12)(€33 — ieza ctg 8) — ness + e3,] + cos (X F 8) [C1 — i12) x 
X (£33 — N? ) — Ain gi (€33 COSO + ieza sin DPL 


Averaging expression (21) over the ion distribution, which is 
assumed to be a Maxwell distribution, we may determine the contribution 
made by ion cyclotron radiation to the emissive power of the plasma 
at the w % suwy, frequencies 


A BslePwne py BETI =z (22) 
TO aheng BEF" 


where ng is the plasma density. 
| 


Cyclotron Radiation of a Slow Magnetosound Wave in a 


Low-Pressure Plasma 


It was noted above that the dispersion equation (3) may be solved 
if the magnetic pressure py is considerably greater than the gasokinetic 
pressure of electrons p, = nol, (this is equivalent to the condition 


Va >> Vg). Let us investigate the cyclotron radiation of a slow 
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magnetosound wave corresponding to the solution (7) of the dispersion 
equation. 


Employing the general expression (20) for the ion radiation inten- 
sity and taking the condition ng << ns into account for the intensity 
of the radiation of a slow magnetosound wave by an ion at. the s-th har- 
monics per unit of solid angle, we obtain the following expression 


2 
ð _ oz Ue 2tst® 
‘sl = “Zc? sl 2 (23) 


where / 19 8 


: n? sin 6 
Uas an 
v2n® sin X| cos x 1—8 
DETS] dn? : 


don, | (i sin(x—O)[o sind x 
dw )| 


‘ R r S > z 
X (ie12 — ites ctg b) Js + (Ps sin 8 — a, cos a) ZAR 1 (Eggfy2 + ` 
A ds 


+ e1123 tg 0) Js + G = E33 COS Î — V g&u Sin e) A -+ cos (X — 0) x 
x [os (Enn? + esa? ,) Ji — (dere — ieog ctg 8) pa tgð +v :) x 
tO ‘ 
. 2 
x nyt Je] }. 


In the case under consideration, the refractive index is determined 
by formula (7), and the damping coefficient is 


2 = 2 
es cos X V3 a Ss 3 (24) 
mme (VE Etim). 


For ions whose velocity is on the order of the mean thermal velocity 
of plasma ions, expression (23) may be simplified (in the case of s = 2, 
33 ess 


onb —~2x (25 ) 
Ws) = One Use ai, 
where 
U s (s8,n 1) nin sin @— 6) 
A T E E TE 
2° (s!)?0 (u — I) a4, sin Xj cosX an? 
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In this case, the equation for the jsaddle point (19) yields the following 
dependence between the angles x and 6: 


tg) = (u— D tex. (26) 


The contribution resulting from magneto-braking radiation of a slow 
magnetosound wave by ions to. the emissive power of the plasma, for fre- 
` quencies close to suwy., is 
q Hi? e wsos ($n) n 3 nSsin 6—X za 


Cr) "7295 s! venga 3 sin Xj cos X|- 


Qs1 = 


(27) 


| v., \2s—2 3 
In order of magnitude, w ~oa) (a): where wọ. is the total /199 
Us ‘ 


Us 


intensity of ion radiation in a vacuum, It was shown in (Ref. 6) that 
the intensities of cyclotron radiation of an Alfvén wave and of a rapid 
magnetosound wave (j = 2, 3) are 


Comparing these results, we find that 
ji va 2541. 
wss \Uy >L 


Thus, in a greatly non-isothermic plasma having a low density, the 
intensity of cyclotron radiation of a slow magnetosound wave for the 


j U ,\2s-+-I 
s-th harmonics is, in order of magnitude, (4) ` times greater than the 
B s 


radiation intensity of an Alfvén wave and a rapid magnetosound wave. 


Cherenkov Radiation of Ma netohydrodynamic Waves by Ions 


In the case of w << wg», the dispersion equation (3) has three solu- 


tions corresponding to magnetohydrodynamic. waves. The refractive indices 
of these waves are on the order of ing (or n AD? and are determined by ex- 


pressions (11), (12). It follows from the eon rie of Cherenkov radia- 
tion Bi ni = 1 that frequencies w << Wyy May exist during radiation 


Bi ~oy re 
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of rapid ions, whose velocity is 


Let us investigate Cherenkov radiation of ions in the low-frequency 
region w << wyy (u >> 1). In order to determine the Cherenkov radiation 


intensity, let us employ the general expressions (20), assuming that 
= 0 and u >> 1. In addition, let us assume that the arguments of the 
Bessel function are small 


Bis ni <1. (28) 


Taking into account (28) for the Cherenkov radiation intensity of 
an Alfvén wave, we obtain the following 


eoo n? sin’ 8 [n38 oi — n?) + Qn? |? 


Woi = = - stone -= 989 8, 
8x0 dîn day e 
1 u2 cos 8 sin X Bee cae : E f(n? — n? 2)? (n3 — n? 3)" (29) 
as 
The damping coefficient of an Alfvén wave is 
2n§ (nh — n?) + ni (nh + nt) + 2an? (nt — 
fom 2 (30) 
msy E o EA 
a= XY M un? nö 
i's 


Equation (19) for the saddle point assumes the following form in this case 


(31) 


nî uns 34 27,2 2 

1a gx + nalns — na) =0. 
It follows from equation (30) that the Alfvén wave is radiated within the 
limits of a narrow cone {x ve along the direction of the outer magnetic 
field. 


The Cherenkov radiation intensity of rapid and slow magnetosound 
waves per unit of solid angle, when condition (28) is fulfilled, is 


2, 3 2 
eo in 0 cos XFO —n 
wors = 3 ars s (X F 0) [Bi ("42,3 s) +2] eT 2*02,3R , 
Bru dn dn (32) 

29 12.3 12. 3 
u cos? 6 sin X cosX—5— - aa 


($e — nis) 


The damping coefficients of magnetosound waves are determined by 
the expression 


— 4f p2 2 2 2 2,2 
mm n(n) 2,3 — Ma tni) +274 anani 


W, 
%23 = -cosy V ga (33) 


ngn? 2,3 (42,3 — 44 3,2) 
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dn js 
The figure presents graphs of the functions =i and nija for 
L 


the case ny > 0 and nA > ng. The points at which the curve dnij inter- 

dny 
sects the line y = + tgx correspond to the solution of equation (19). 201 
As may be seen from the figure, a slow magnetosound wave (j = 2) repre- 


sents a cylindrically converging wave. There are thus two waves for 
each angle x < Xmax» In the case of x > Xmax» a slow magnetosound wave 


is not radiated. The angle x = Xone is determined by the following 


equations 


The rapid magnetosound wave (j = 3) represents a cylindrically 
diverging wave. For any angle x, one such wave is radiated. Expres- 
sion (32) for a slow magnetosound wave is not valid in the limiting 
case, when nj > © and ni > nk + nê. In this case, we must employ the 


general expression (20) for u >> 1 and s = 0. (We should point out that 
Cherenkov radiation in a direction which is perpendicular to the direc- 
tion of the outer magnetic field may be absent since the condition 

Bij nį = 1 thus contradicts the assumption w << üg.) If ny < ng, then 


the substitution ną Z ng must be made in the figure. In the case of 


ny, < 0, the functions dn ij change their sign. 


dny 


Electron Cherenkov Radiation 


It is of interest to study Cherenkov radiation of low frequency 
waves (w % wy,) by electrons of a non-isothermic plasma, because under 


certain conditions it may make the main contribution to the over-all 
plasma radiation. In particular, this is related to the fact that the 
intensity of ion cyclotron radiation sharply decreases as one recedes 
from the center of the line w = Sugi O e E 


We obtained the expressions for the intensity of Cherenkov radia- /202 
tion, by an electron moving in the plasma along a spiral, of each of 
three normal waves. We determined the contribution made by Cherenkov 
electron radiation to the emissive power of a non-isothermic plasma at 
the w v wy; frequencies. 


We shall employ the general expressions given in (Ref. 2) for s = 0 
in order to determine the components of an electromagnetic field 
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Figure 1 


produced by an electron in the wave zone. As a result, we obtain 


gi’ 


E; = P;sin (x + a7) E 


Ey = Pj cos (x J af] Ey a 
Hy = + Pjnjcos (x + a7) Ho: : 
Hy = — nicos (x F 8) Ep; 
where 
W=Rsinx [ki +kij(ki)ctgyl F ot a 
RTA dn, ;\ 7 
= sgn (cos x In? ) ; 
Eşi == 04 (mien + ni j€93 — €11€3a) Jo v5 [nina (iere — ier, ctg 8) + 
$ + ey if23] Jo 
Ex; == + cosX fv [ni (ie;2 — iez ctg 8) + ieee ys] Jo— vı inina (np = Pp? 


~ 11) -H £z2€23] Jo} — sin X {vı [n,n (dey2 — isz; ctg 0) + &1 48230 + 
+0, [njn,7 —nje,, (1 + cos? 6) + efi + etala]; 
Hi = — U, COSÌ (ie,2€35 + £11l823 tg 0) Jo +v 1,sin 8 len(a = 
— £1) + &,28, ctg 0 — efl Jy: 


In (34), vj and vy are, respectively, the electron velocity components 
which are parallel and perpendicular to the direction of the outer 
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magnetic field. The argument of the Bessel functions and their deriva- 


v 
tives equalskj ro, where rg = L is the Larmor electron radius. The 
e 
remaining notation was presented above. 


Employing furmulas (34) and (35), we may obtain the following ex- 
pression for the intensity of electron Cherenkov radiation per unit of 
solid angle 

nee 2 
Wj = - [EH H; + cos (x F o) (£ oi) |. (36) 


ProB Au 


The argument of the Bessel function is small for electrons having 
a velocity v| which is on the order of the mean thermal velocity of 


plasma electrons (v, ~% vre = Te); 
t m 


n$ 
Riro = a Ski. 


In this case, expression (36) may be simplified 


yet DS oe = 
: 2nv 3 : dn yj tl 27,2 2 PY (37) 
o (1 — u) sin X | cos X int w (nii — ng) (nna) 


where 
U 20: ee 
j = n? sin 9 {cos x ai a T a n +yz] +sinX x 


u(i — 
y ucos (X F 0) 
~ ul rive + u(i — u) x 


v 2 
x Toa UP is . 


x [n in? + n + 


Averaging (37) over the electron distribution, which is assumed to 
be a Maxwell distribution, we may find the contribution made by electron 
Cherenkov radiation to the plasma emissive power at the frequencies 


w v agg? 


n= tan , m wok 7 
(2n)'/*0,, 2 \27 2 2 \2 
o(1—u)sinX cost ol H| (n?,—2n2,) (nir — Re) (38) 
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Comparing the contribution made, in order of magnitude, by ion cyclotron 
radiation (22) and by electron Cherenkov radiation to the plasma emissive 


power, we obtain 
Asi Ons 2s—1 One : 
uy Po Yo (39) 


It can be seen from (39) that the relative contribution of the elec- /204 
tron Cherenkov radiation increases with an increase in the harmonic 
number s, since it was assumed from the beginning that VTi << Vo << Vre’ 


Let us investigate the Cherenkov radiation of a slow magnetosound 
wave by electrons in a low-pressure plasma (n, >> na). We shall assume 


that the argument of the Bessel functions in (35) is small (kj; rọ << 1). 
This can be fulfilled only for very rapid electrons (v| >> VTe)». Employing 


the general expressions for the electromagnetic field components (34), 
(35), the equations (7), (26), and also the condition of electron Cheren- 
kov radiation, we finally obtain 


w, Toke [tet t+ (1-H)? (1 — Bandy? ek (40) 
1°" 220) teh (sec? x — 8? n2)? B? n?nĝ | cos% ` 


. Expression (40) determines the intensity of Cherenkov radiation for 
an electron having the velocity 8} per unit of solid angle in the direc- 


tion x. If By ng < 1, then the radiated wave represents a cylindrically 
converging wave (for x To ). The waves may thus be radiated in any direc- 


tion x. The frequency of a wave radiated at the angle x is 
tg 
v = oot SUES es . 
V | sec? x — 8? n2] 


In another limiting case, when Bj ng > 1, the radiated wave repre- 


sents a cylindrically diverging wave. The vector of the group velocity 
is thus directed at the angle x > Xmin» Where 


(41) 


X 


mince are sec (Pants). 


The phase velocity in this case is directed at the angle 6 < 6, .. 
in the direction of the outer magnetic field, while Goa = Xmin- However, 


we cannot employ the expressions obtained in the case of x > Xmin» since 
in this case w > œ», according to (41). , 
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Averaging (40) over the electron distribution, we obtain the 
contribution produced by Cherenkov radiation, by plasma electrons, of 
a slow magnetosound wave to the plasma emissive power: 


a= = S nen 
(2n)3/205,y nnsin WV | sec? X — BY n? (42) 
We should recall that at any angle x > Xmin’ electromagnetic waves with 205 
two frequencies are radiated -- one with the frequency w > wy; (electrons 


with 68) ng > 1) and the other with the frequency w < wyi (electrons with 
the velocity vj < Vg) Expressions have been obtained (Ref. 6), in which 
it must be assumed that Zp, << 1, for the intensity of Cherenkov radiation 


by electrons of an Alfvén wave and a rapid magnetosound wave in a non- 
isothermic plasma with a low density. 
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BRAKING OF RELATIVISTIC PARTICLES IN LOW ATMOSPHERIC LAYERS 
V. B. Krasovitskiy, V. I. Kurilko 


As is known, when a relativistic charged particle moves in an external 
field, the braking force by radiation increases proportionally to the 
square of the particle energy, and may be greater than the Lorentz force 
(Ref. 1). The influence of the medium on the braking radiation is of 
great interest, since -- when particles move in the lower atmospheric 
layers -- the braking force by radiation may differ significantly from 
the vacuum force. The expression for the braking force by radiation of 
an oscillator in a medium was obtained and studied in the non-relativis- 
tic approximation in (Ref. 2, 3). The spectral radiation density was 
studied in the relativistic case in (Ref. 4 - 7). 


In order to investigate the influence of the medium on the braking 
radiation, it is necessary to know the total energy losses of the particle, 
i.e., it is necessary to specify the medium properties. We shall formulate 
a model for a medium with an isotropic dielectric with a given dispersion /206 


2 
vo 
(eo) lars 
4me*ng 
where 2 are the atom resonance frequencies; w = —— no -- density 
m 


of the medium particles. 


The following expression was obtained in (Ref. 4 — 7) for the spec- 


tral density of particle radiation in the frequency region w >> wn (ug = 


wg Yg _ g2 ) 


: (1) 


P(e) = SPG pe) ft +20) Jy, O1— Font) — 
£ 


<< Jy, (7)] ar) Be > 1, 


where t=2. a 1 — Be)"; J, (8) and K (6) are the Bessel functions. 
ERSA 


Let us employ formula (1) to determine the total losses by radia- 
tion as a function of the particle energy. Since it is difficult to 
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integrate equation (1) in the general case, let us find the frequency 
regions in which the spectral density has the maxima (g v 1), and let 
us determine the radiation intensity in each of these regions. 


1. In the low-energy region, when the following inequality is ful- 
filled 


v wt 
a << BK E=(1— py, (2) 


the frequency On ~ Oy (1 - g2)7! ~~ which makes the maximum contribution 


to the vacuum (Ref. 8, 9) -- is small as compared with the resonance fre- 
quency wn << Q. In this case, the radiation spectral intensity has 
maxima at the frequencies wp and 2. The braking radiation at the fre- 


quencies wp causes losses which coincide with vacuum losses 


2, w2 
WoT e (3) 


The maximum close to the frequency determined from condition B2e (w) X 1 
is also caused by the curvature of the particle trajectory. However, in 207 
this case the influence of the medium is significant, as a result of which 
the radiation is considerably greater than in the preceding case 


E?) s 
~ 3205 (=). (4) 


Thus, the braking radiation is larger (in the energy region under 
consideration) than in a vacuum, and increases with the energy E propor- 


4 
tionally to E /3, 
2. The inequality (2) assumes the following form for given character- 


istics of the medium and unchanged magnetic field strength with an increase 
in the particle energy: 


oy 1 o (5) 
2 <p<aK<i. 
In this case, the braking radiation maximum is located at the frequency 
3, 
we eH 1, The corresponding loss is 
%9 


Tat EF . (6) 
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3. With an increase in the particle energy, the relationship 
between the particle parameters and the medium assumes the following 
form 


1 ao o2 
aKa KaL (7) 


Thus, the radiation maximum in a vacuum occurs in the frequency region 
w >> Q. The force of the braking radiation depends on the parameter 


wW 
u= ae Lin this case. 


OF 


In the case of u >> 1, the braking radiation intensity in the fre- 
quency region w % w, is exponentially small as compared with the corres- 


ponding value in a vacuum. The braking radiation intensity is also small 
WV Qe 
(r ~ ar) In the opposite case (uy << 1), the influence of the medium 


on the braking radiation may be disregarded, so that the particle losses 
coincide with vacuum losses. 


Thus, the medium influences the particle braking radiation only 
when the inequality wi >> Quy is fulfilled. In the region of relatively 


small energies (<3) » the presence of the medium leads to an increase 
6), 
0 


g2 
in the braking radiation, and in the energy region E2 >> —> it leads to 
w 
0 
a decrease in the braking radiation, as compared with a vacuum. Thus, 
the braking radiation is significantly less than the vacuum radiation in 


the energy region F DED Cae 


The change in the braking radiation intensity in a medium may be 
explained in physical terms as follows. In the presence of the medium, 
the frequency determining the maximum of the radiation spectral density, 
according to formulas (1), depends on the medium parameters 


e 


on (9) = Gp thO = on. 


For small particle energies, the presence of the medium leads to an in- 

crease in the frequency WCE Cum) > 1, and the corresponding wavelength 
= “0 

Ay = — decreases): 
€ 


Wm < wn < Q. Therefore, the braking radiation 
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4 
intensity, which is proportional to (S) 13, also increases. For large 
particle energies (uw, >> 2), the frequency wn decreases considerably: 
wt << Q << w, (the characteristic wavelength i, increases, e (wn) < 1). 


Correspondingly, the braking radiation intensity decreases. The 
Cherenkov radiation, which can be computed from the customary formulas 
in the case under consideration, is larger than the braking radiation 
throughout the entire energy region where the influence of the medium 
must be taken into account. 
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EXCITATION OF WAVES IN A CONFINED PLASMA 
BY MODULATED CURRENTS 


A. N. Kondratenko 


The excitation of waves in a plasma by modulated currents has been /209 
studied repeatedly (Ref. 1 - 3). In all of these studies, it was 
assumed that the plasma was an unconfined, isothermic, and linear 
medium. In this article we shall study the excitation of waves in a 
confined plasma, and we shall attempt to take its nonlinearity into 
account. 


When solving the problem of wave excitation in a plasma by a 
modulated current j, we assume that the current is given, and we shall 
disregard the inverse influence of the wave on the motion of current 
particles. This is valid as long as the energy losses of the current 
particles at the wave length are negligibly small as compared with the 
energy itself. The electric field strength of the wave is determined 
by the equation LE = j, where L is the differential operator, which is 
generally speaking nonlinear. A solution of this equation leads to the 
following value for the Fourier field components: E = i. If the fre- 


quency of the current modulation is such that A vanishes (dispersion 
equation), then the field strength becomes large. It is usually assumed 
that A is limited by particle collisions or Cherenkov absorption of wave 
energy by plasma particles. This is the result of the linear theory, 
and it is valid for insignificant amplitudes of the wave field strength. 
However, in resonance, when A ~% 0, the field strength greatly increases, 
and the linear theory may be invalid. Even in the case of slight non- 
linearity, the limitation on the amplitude of the wave field strength 
due to nonlinear interaction may be more substantial than the limitation 
imposed by dissipation. 


Let us study the plasma layer which is infinite in two directions 
(y, z) and bounded by two parallel metallic plates in a third direction. 
The distance between these plates is 2a. A modulated current having 
the form of an infinite layer of thickness 2b, b < a,moves in the plasma 


along the Z axis. The plasma is located in a constant magnetic field 
directed along the current, which is so strong that the plasma particle 
motion across the field may be disregarded. 


Let us investigate two problems. In the first problem, let us 
determine the wave field strength in the linear approximation for an 
arbitrary dispersion law. In the second problem, let us determine the 210 
wave field strength in the nonlinear, but hydrodynamic approximation. 
We shall assume that the nonlinearity is slight. The following quantity 
is a small parameter of the problem 
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RAK, ey 
~ moV 72)? 
z (1 — 45.) 


where e, and m, are,. respectively, the charge and mass of particles of 
type a (a = i -- ions, a = e -- electrons); Eg -- amplitude of wave 
field strength in the linear approximation; w -- modulation frequency; 
Vg -- phase velocity of a propagated wave, equalling the velocity of 


current particles; Uta = YTo (vpy ~- mean thermal velocity of a particles). 
v 


We can write the system of equations describing wave excitation in 
a plasma under the conditions being considered as follows 


PE aE oF ot | v2 a oem, 
ae + ow? BE p ye 2 va z do + owVy X oa” =, * 
asst a=í e(1— ura) 
3 2 Ue 2\ 9 ua 
x | ( —u DEE aa e (1) 
u, = eE l+ u2, ð u? u2 ð u? 
% mo» (1 —- uz.) l= ae a R a: (2) 
; : 1- 8? Vo 
where E - z = component of the wave electric field; o = 79 B= 5 
4ne2ng 
c -- speed of light}; € = wt - k3z; 92 = ——— (np —- equilibrium 
Mg 


density of plasma particles which is equal for ions and electrons); 
w = wt ivy (vq -- frequency of collisions between a particles); foa -- 


equilibrium distribution function of a particles normalized to unity; 


Uy = Ya (v, -- hydrodynamic velocity of a particles). 


Vo 


The derivation of equations (1) and (2) is given in (Ref. 4). In 
these equations, we must set y = 1, and we must supplement them with 
the boundary conditions, which in this case may be reduced to setting 
the field strength on the wave guide walis equal to zero: 


We can present the modulated current in the following form 211 
A : $ 7 CPV or 0< (x) <b 
J=h (x) sin, h= | 0, b < (x) <a, (4) 
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where e -- charge; po -- density; Vg -- velocity of current particles. 


We can obtain the equation of the linear approximation for the 
field from equation (1), setting u, = 0 in it. The solution of the 
linear equation can be written as E = R(x) cos €. We obtain the 
following equation for R(x) 


Sa , 
ae + k?R = — frooji (5) 
where 
e ga ( VIA do 
2 214 y Qa” $ ðv 
meme ra m j) n= (5a) 


Solving equations (5) and (3) together, we find 


R (x)= È Ra COS anX, (6) 


where ana y 
Tja. Vo ci : 
Ra = 7 In = 2 a sin anb; 


a 


an = 55 (20 + 1), An = a2 — Ri. 


If the velocity of the current particles is close to the mean 
thermal velocity of plasma ions or electrons, the integral in formula 


(5a) must be determined numerically, Im kf and Re ki are equal in order 


of magnitude, damping is large, amplitude of the excited wave field 
strength is small, and the energy losses by the current are insignificant. 
In this case, we can confine ourselves to the linear theory. 


If Vo >> Vte,i OF VTe >> Vo >> Vri» then Im ki << Re ki and at a 
certain frequency of current modulation |A| ~ Im k? << dn may hold. 
The corresponding R, (we shall designate it by Eo) increases greatly, 


and only the resonance component can remain in the sum (6): R(x) % 
% Ep cos apx. 


Since the amplitude of the wave field strength is large, the non- 
linearity of the medium may be significant. When determining the non- 
linear dispersion equation, let us regard the nonlinearity and thermal 
scatter as small independent additions to the linear hydrodynamic 
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dispersion equation (A, = 0, a = Re k$). Therefore, the wave damping /212 
P q n L P foe 


caused by kinetic phenomena can only be taken into consideration in 
the third approximation. 


We may find the hydrodynamic velocity of the linear approximation, 
corresponding to A % 0, from equation (2) 


ul!) == £, COS a,x sin Ẹ. (7) 


Let us substitute (7) in equation (1), and let us retain the terms ve, 


Then, representing the field E in the second approximation in the form 


pC?) = EoR2 (x) sin 2%, we obtain the following equation for Rp» (x) 


dR (8) 
Gt + 9°Re = fe cos? anx, 
where 
.e 
ne | So) a oe 
d ( “ChE a h=- Son, oath rae 
=t a=i 
Solving equations (8) and (3) together, we find 
TE h, 2a? C08 gx fy COS2a„x (9) 
27g 4a? g?  cosqa 2” ta? Gt a 
In the second approximation, the hydrodynamic velocity is 
1+ u? e S i 
uf?) = — zeori ; Fe geosthus| cos 2. (10) 
u, 


After substituting ul), ul?) in the right part of equation (1) of 
a 


the first and second approximation, we obtain the following equation for 
the field E in the third approximation 


aE) : i Ed g2 
: 2 — FR ), E a z 
Ox? + oe E ») TERE N uZ) Az cos 3y 


9 — 8u? + 3u4, G 


E2? G 
-oA a 6A) R Hy 


do -AF 


X cos? a,x } cos anx cos $, 


where A3 is a certain function of x which is determined by the third 
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harmonic with respect to &. 


We can write the solution of equation (11) in the form ECG) = /213 
= EoR3(x) cos £. We obtain the following equation for R3(x) 


dR ; 
TE, + k? Ra = Qi cos anx + Qa cos gX COS an% + Qs cos 3an, (12) 


where 


: o e€ fi 8a? — 34? 
Ša ) ELLE N A a a A OA 
Qı - 4 (1 Zy É (3 Uza 4a2 — q? + g Ea X 


9— 8u?, + 3u4, ] 


l= T 
Q ae z aea h 3 u2 2a? 1 
a AlMv—2\2 g@ = -—- e — . 
a 4(1—u?,)? z ( za) 4a? — q? cos qa’ 
° ‘92 
Q= oY Site f fe Smee yea 9 Bu t a, 
saata Gair eee Soe A 
a 4(1—ui,) 4 4a? — q? 8 1—u2, 


Since Qi 2,3 % e*, we can write the solution of equation (12) in 


the following form (Ref. 5) 


Rs = cos p + B, cos (an + q)x + B2 cos (an — q) x + Bs cos 3anx, 
A =k, + eA. 


We thus find 


= _ Qa __ @ cos (an + 9)* cos (an —9)*] _ 
R; cos(ka 5 js DD ee Tor 


_ Qs 
2 
n 


COS 3a,Xx. 


Employing the boundary condition (3), we find that 


2a,, sin ga 


4a? —¢ 


=0. (13) 


\ 
cos( ~ x a—(— Ip Ss : 


We can write the solution of equation (13) as follows 


kı =a, + ô, b~ <a, 


d 4a? $i 
We fin 2and = Q, + ee „Singa 
4a — g -qa 
(14) 
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Thus, we can write the nonlinear hydrodynamic dispersion equation 
as follows 


2 
2. 


j? | — l a ae I 
ow | 1+ Le | = a? + 2anð, (15) 


and v, includes both pair collisions and Cherenkov wave absorption by 


plasma particles. 


In the case of resonance, when Re kj = aĝ, the field strength ampli- 


tude of the excited wave is 


— Anws}j, 


A Fe y (16) 
mar (eDes-3) | 


— Uia 


€ 


Let us first examine high frequency oscillations: uti e << 1. Dis- 
> 


regarding the electron mass as compared with the ion mass, we obtain 


10w? — 2 2 2\3 
24,8 = — o6?Q2.8,; 8, = 3; 1—2 O (20) „ tgga. (17) 
8 4079-93, Qe (40*— 23) ga 


We can show that in the case of w< Qoe, 51 > O always holds. We obtain 


the following value from the equation (15) for the square of the phase 
velocity 


Baar 
Q2, — o? + ca? 


dateng, | (18) 


2 —o? S A 
[oz OE E cadea 
Since 6, > 0, the phase velocity increases with an increase in the field 


strength amplitude of the wave. We may use formula (16) to determine 
the amplitude of the field strength: 


leas BxwepyV o sin a,b 
aif) +(y]" (19) 


Ve 
If e26; > T3 the maximum amplitude is determined by the nonlinearity: 


Snr 2w°V Žo sina,b (20) 
e023, ” a nî * 


_— 
E = 
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The low frequency oscillations: up, >> 1 >> Uqi.’ 


The amplitude of the wave field strength in this case is 
8xwen,V, sine,b 
E = =e = e I 
2 [f 2y \2 ye aa ° (21) 
afer 


m 
If uel ( = x) then ô> and the phase velocity have the values 


of (17) and (18), respectively, for a replacement of the indices e > i. 
2 


. Q 
In the case of sound oscillations, when K, we have 
Bue 
a= zlo 5 226) 0% | tea, cs 
72 39 \Y Ge g qa’ 
2 2 
. on 25, + ant Dze 
V . 

r È < ef8ys the field strength amplitude is limited by the nonlinearity 


and has the same form as in the case of high frequency oscillations (2), 
with a replacement of the indices e >i, 1 > 2. 
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